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Abstract 

This thesis describes the observation and study of Bose-Einstein condensation (BEC) 
of magnetically trapped atomic hydrogen. The sample is cooled by magnetic sad- 
dlepoint and radio frequency evaporation and is studied by laser spectroscopy of the 
1S-2S transition in both the Doppler-free and Doppler-sensitive configuration. A 
cold collision frequency shift is exploited to infer the density of both the condensate 
and the non-condensed fraction of the sample. Condensates containing 10^ atoms 
are observed in trapped samples containing 5 x 10^*^ atoms. The small equilibrium 
condensate fractions are understood to arise from the very small repulsive interac- 
tion energy among the condensate atoms and the low evaporative cooling rate, both 
related to hydrogen's anomalously small ground state s-wavc scattering length. Loss 
from the condensate by dipolar spin-relaxation is counteracted by replenishment from 
the non-condensed portion of the sample, allowing condensates to exist more than 
15 s. A simple computer model of the degenerate system agrees well with the data. 
The large condensates and much larger thermal reservoirs should be very useful for 
the creation of bright coherent atomic beams. Several experiments to improve and 
utilize the condensates are suggested. 

Attainment of BEC in hydrogen required application of the rf evaporation tech- 
nique in order to overcome inefficiencies associated with one-dimensional evaporation 
over the magnetic field saddlepoint which confines the sample axially. The cryogenic 
apparatus (100 mK) had to be redesigned to accomodate the rf fields of many milli- 
gauss strength. This required the removal of good electrical conductors from the cell, 
and the use of a superfluid liquid helium jacket for heat transport. Measurements of 
heat transport and rf field strength are presented. 

The rf fields in the apparatus allow rf ejection spectroscopy to be used to measure 
the trap minimum as well as the temperature of the sample. 
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Title: Lester Wolfe Professor of Physics 

Thesis Supervisor: Thomas J. Greytak 
Title: Professor of Physics 
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Chapter 1 



Introduction 



1.1 The Lure of Bose-Einstein Condensation 

Bose-Einstein condensation (BEC) of a dilute gas has been a very important goal 
since the beginning of experimental research on spin-polarized atomic hydrogen. The 
original intent [|I| was to study quantum phase transition behavior and search for 
superfluidity in this dilute gas, a weakly interacting system that is much more the- 
oretically tractable than strongly interacting degenerate quantum systems such as 
superfluid liquid "^He. This phase transition is a purely quantum mechanical effect, 
and, unlike all other phase transitions, it occurs in the absence of any interparticle 
interactions. 

Although the search for BEC in dilute gases began in spin-polarized hydrogen, 
the first observations of the effect used dilute gases of the alkali metal atoms Rb 
0, Na 0, and Li 0, || in the summer of 1995. The properties of hydrogen that 
were so attractive in the beginning turned out, in the end, to be irrelevant or even 
a hindrance. Nevertheless, we have finally observed BEC in hydrogen [0, 0. In 
several ways our hydrogen condensates are rather different from those of the alkalis, 
and we use different techniques to probe the sample. Hydrogen thus contributes 
significantly to the extraordinary flurry of experimental activity creating and probing 
Bose condensates, and applying them to intriguing new areas. 

The formation of a Bose condensate out of a gas can be studied as the occurrence 
of a quantum mechanical phase transition, but perhaps even more interesting are 
the properties of the condensate itself. Centered in the middle of a trapped gas, 
this collection of atoms exhibits such bizarre effects as a vanishingly small kinetic 
energy, long range coherence across the condensate , a single macroscopic quantum 
state, immiscibilty of two components of a gas |^ |T0|, and constructive/destructive 



interference between atoms that have fallen many millimeters out of the cloud ||11|| . 
These systems hold promise for creating quantum entanglement of huge numbers of 
atoms for use in quantum computing and vastly improved fundamental measurements. 
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1.2 The Significance of the Work in this Thesis 



In this thesis we describe the culmination of a 22 year research effort: the first ob- 
servations of BEC in hydrogenQ. Furthermore, we demonstrate a new technique of 
probing Bose condensates, optical spectroscopy. This tool allows us to measure the 
density and momentum distributions of the sample, and thus infer the temperature, 
size of the condensate, and other properties. We find that the condensates we create 
contain nearly two orders of magnitude more atoms than other BEC systems, and 
the condensates occupy a much larger volume. In addition, we find that nearly 10^° 
atoms are condensed during the 10 s lifetime of the condensate. This high flux of 
condensate atoms implies that this system might be ideal for creation of bright, co- 
herent atomic beams. Indeed, hydrogen has proven to be a rich and promising atom 
for further studies of BEC. 

The thesis is divided as follows. Chapter || and the associated appendices provide 
a detailed description of the behavior of the trapped gas both in the classical and 
quantum mechanical regimes. Results are obtained to guide the reader's intuition 
throughout the remainder of the thesis. Chapter ^ describes changes to the apparatus 
that allowed BEC to be achieved. Chapter ^ details use of the new capabilities of 
the improved apparatus. The most important contributions described in this thesis 
are found in chapter ^ where studies of BEC are presented. Finally, conclusions and 
suggestions for future work are made in chapter ||. 

Much of the work described here was done in close collaboration with a group of 
people, as shown by the several names on our papers. In particular, this thesis should 
be considered a companion thesis to that of Thomas Killian |]13|. Both should be 
read to obtain a complete picture of the whole experiment. The work builds on ear- 
lier work by Claudio Cesar, whose thesis fl^ describes the first 1S'-2S' spectroscopy 
of trapped atomic hydrogen; Albert Yu, whose thesis |T^ on universal quantum re- 
flection includes important insights into non-ergodic atomic motion in the trap; Jon 
Sandberg, whose thesis describes the physics of 1S'-2S' spectroscopy in a trap and 
the development of the uv laser system we use; and John Doyle, whose thesis ||17 



describes trapping and cooling of spin-polarized hydrogen and the magnet system we 
use to trap the gas. 



1.3 The Basics of Trapping and Coohng Hydrogen 

There is a long and interesting history of spin-polarized hydrogen in the laboratory, 
summarized well by Greytak |jl8|| . Here we briefly summarize our method of creating 
spin-polarized hydrogen, loading it into a trap, and cooling it to BEC (quantum 



degeneracy). A broader description is in the thesis of Killian ||13 |. 

Hydrogen is composed of a proton and electron, whose spins can couple in four 
ways: the four hyperfine states are labeled a-d (these states are derived in section 



two dimensional analog of BEC has been observed very recently in atomic hydrogen confined 
to a surface of liquid ^Hc |l2[ 
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Figure 1-1: Hyperfine-Zeeman diagram of ground state atomic hydrogen. The eigen- 
states of the combined hyperfine and Zeeman Hamiltonian are described in section 
3?T|. Atoms in states c and d can be trapped in a magnetic field minimum. 



p.l.l| ). The high- field seeking states (a and h) are pulled toward regions of high 
magnetic field, and the low-field seeking states (c and d) are expelled from high field 
regions. Figure |1-1| shows the energies of the four states as a function of field strength. 
Our experiments use the pure, doubly-polarized (i-state. 

At the beginning of an experimental run molecular hydrogen is loaded into the 
cryogenic apparatus by blowing a mixture of H2 and ^He into a cold (T ~ 1 K) 
can, called the dissociator, where it freezes. Then, for each loading of the trap, H2 
molecules are dissociated by pulsing an rf discharge (the dissociator is in a region 
of high magnetic field, 4 T). The low field seekers are blown into a confinement cell 
(4 cm diameter, 60 cm length), over which the inhomogeneous trapping field is super- 
imposed. The trap fields are created by currents in superconducting coils, described 
by Doyle Jl^ and Sandberg |16|. The coils create a trap with maximum depth 0.82 T, 
which, for hydrogen's magnetic moment of 1 Bohr magneton, corresponds to an en- 
ergy et/ks = 550 mK. (See appendix | for conversions between various manifestations 
of energy in this experiment.) The trap depth is the difference between the field at 
the walls of the containment cell and the minimum field in the middle of the cell. 
There are seventeen independently controlled coils in the apparatus which are used 
to adjust the trap shape. Figure p.-2| is a cutaway diagram of the system. Expanded 



views of the top and bottom of the cell are in figures 3-8 and 3-9. The trap shape 



used for capturing the atoms is basically a long trough. The field increases linearly 
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Figure 1-2: Cutaway view of apparatus. Molecular hydrogen is dissociated in an rf 
discharge in a 4 tesla region. The atoms move into the confinement cell which has 
superimposed on it a trap, created by the quadrupole coils (radial confinement) and 
the pinch coils (axial confinement). The coils arc bathed in 4.2 K liquid helium. The 
cell is centered in the vacuum can with three Kevlar strings. It is thermally anchored 
to the mixing chamber of a dilution refrigerator (not shown). The laser beam passes 
through the upper window, retroreflects on the mirror at the bottom of the cell, and 
exits through the upper window. Atomic fluorescence photons pass through the MgF 
window and Lyman-a filter, and are counted on the microchannel plate detector. 
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away from the z axis of the trap (the potential exhibits near cyhndrical symmetry 
about the z axis); the potential is small and roughly uniform for about 20 cm along 
the z axis. "Pinch" coils at each end provide the axial confinement. The field profile 
indicated in figure |1-2| is along the z axis of the cell. 

The trap depth sets an energy scale: for atoms to be trapped their total energy 
must be less than the trap depth. Two techniques are used to cool the atoms into the 
trap. First, after H2 molecules are dissociated [0, the hot atoms are thermalized 
on the walls of the cell, held at 275 mK. In order to prevent the atoms from sticking 
tightly to the cold surfaces, a superfluid ^He film covers the walls and reduces the 
binding energy of the H atoms to a manageable 1.0 K pO| , ^ WA - When the atoms 



leave the wall their total energy {?>kBT /2 of kinetic energy and et of potential energy) 
is still greater than the trap depth. The second stage of cooling into the trap involves 
collisions among the atoms that are crossing the trapping region. Sometimes these 
collisions result in one atom having low enough energy to become trapped. The 
partner atom in the collision goes to the wall and is thermalized. The gas densities 
expected in the cell correspond to a collision length many times greater than the cell 
diameter, so atoms pass through the trapping region many times before becoming 
trapped. A recent study of the trap loading process will be published elsewhere [ff^ 



After loading the trap, the cell walls are quickly cooled to below 150 mK to ther- 
mally disconnect the sample from the confinement cell. At these lower temperatures 
the residence time of an atom on the surface of the cell becomes much longer than the 
characteristic H+H— > H2 recombination time, and so the surface is "sticky"; all the 
atoms that reach the surface recombine before having a chance to leave the surface. 
Thermal disconnect between the cryogenic cell and the trapped sample is achieved 
because no warm particles can leave the walls and carry energy to the trapped gas. 
Atoms can leave the trapped gas and reach the wall, however, if they have large 
enough total energy to climb the magnetic hill. As these very energetic atoms irre- 
versibly leave the trapped sample, it cools. In fact, this preferential removal of the 
most energetic atoms, called evaporation, is very useful. It is the cooling mechanism 
we exploit to create samples as cold as 20 /iK, as described in chapter ^. A thorough 
description of the heating and cooling processes which set the temperature of the 
system is found in appendix ^ 

We catch both the d and c low field seeking states in the trap. However, inelastic 
collision processes involving two c state atoms quickly deplete the c-state population, 
and the remaining atoms constitute a doubly polarized sample (both electron and 
proton spins are polarized). Because the spin state is pure, the spin-relaxation rate 
is rather small. For peak densities in the normal gas of n ~ 10^^ cm~^ (nearly the 
largest for these experiments) the characteristic decay time is 40 s. 

The trapped gas is probed using various techniques. The simplest method involves 
quickly lowering the magnetic confinement barrier at one end of the trap, and mea- 
suring the flux of escaping atoms as a function of barrier height [1^, In a limited 
regime this technique can be used to infer the gas temperature and density. This tech- 
nique, and its limitations, is described in chapter ^. Another technique is rf ejection 
spectroscopy, described in chapter ^. The third probe technique is laser spectroscopy 
of the 1S-2S transition, envisioned [0, |16[, realized |jl^, |2^, and perfected [|1^ in our 
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lab for a trapped gas. This diagnostic tool can give the density and temperature of 
the gas over a wider range of operating conditions than the other techniques, and has 
proven crucial for studies of BEC. The laser probe will be described in more detail in 
chapter ^, and is a principle subject of the companion thesis by Killian |jl3 . 
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Chapter 2 

Theoretical Considerations 



This chapter addresses two important topics from a theoretical perspective. First, we 
explain why previous experiments could not bring trapped hydrogen into the quantum 
degenerate regime. This understanding allowed us to solve the problem and achieve 
BEC. The second topic is related to the small fraction (only a few percent) of atoms 
we are able to put into the condensate. This situation differs markedly from that 
observed in Rb and Na BEC experiments. An explanation of this difference and 
others is the second topic we explore theoretically. 

The starting point for understanding the behavior of the trapped gas is classical 
statistical mechanics. Good explanations have been developed elsewhere p6| , |17, 
For completeness, we present a development in appendix Since the gas 



exists in a trap of finite depth e^, we use a truncated Maxwell- Boltzmann occupation 
function. Knowledge of the trap shape allows us to calculate the population, total 
energy, density, etc. The effects of truncation of the energy distribution at the trap 
depth are included in these derivations. We see that for deep traps (rj = et/ksT > 7) 
the truncation does not significantly alter the properties of the system from estimates 
based on an untruncated distribution. Appendix |^ also explains the density of states 
functions (see section [A.1.2| ) used in the remainder of the thesis and summarizes 



evaporative cooling (see section |A.2 



2.1 Dimensionality of Evaporation 



The temperature of the gas is set by a balance between heating and cooling processes 
as described in section A. 2. 2. Previous attempts to attain BEC in hydrogen |29 



failed because the cooling process became bottlenecked by the slow rate at which 
energetic atoms could escape, thus reducing the effective cooling rate. To understand 
this bottleneck we must first consider the details of the trap shape. We then study 
the motion of the particles that have enough energy to escape. 

The trap shape used to confine samples at T < 200 fiK is often called the "loffe- 
Pritchard" [0, |2^ type (labeled "IP"). Using axial coordinate z and radial coordinate 
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p, the potential has the form 



Vjp{p,z) = ^{apY + {Pz^ + eY-e (2.1) 

with radial potential energy gradient a (units energy/length), axial potential energy 
curvature 2/? (units of energy/length^), and bias potential energy 6. See section |A.1.2 
for a summary of the density of states functions for this trap. 

In the limit of p <C 9 /a, the loffe-Pritchard potential is harmonic in the radial 
coordinate, as may be seen by expanding the potential in powers of ap/{(3z'^ + 6): 

The trap is harmonic in the radial direction when the third term is much smaller 
than the second term. This is true for radial coordinates p <^ Panharm, where we have 
defined the "anharmonic radius" at which the second term matches the third term: 

Panharm — • V'^''^/ 

a 

The trap appears harmonic in all three directions to short samples for which the radial 
oscillation frequency is essentially uniform along the length of the sample. This occurs 
for temperatures T ^ 46 /ks- In the harmonic regime, the axial oscillation frequency 
is 

^z = \—, (2.4) 
V m 

and the radial oscillation frequency is 



For the evaporation to proceed efficiently, atoms with energy greater than the 
trap depth (called "energetic atoms" ) must leave the trap promptly, before having 
a collision. As shown in appendix |B], in the vast majority of cases an energetic atom 
that collides with a trapped atom will become trapped again. The rare collision that 
produced the energetic atom will be "wasted" . It is essential to understand the details 
of the particle removal process if maximum evaporation efficiency is to be achieved. 

Previous attempts to cool hydrogen to BEC utilized "saddlepoint evaporation", 
in which energetic atoms escape over a saddlepoint in the magnetic field barrier at 
one end of the trap. To escape, the atom must have energy in the axial degree of 
freedom [z) that is greater than e^. This atom removal technique is inherently one 
dimensional. The collisions which drive evaporation produce many atoms with high 
energy in the radial degrees of freedom, and in order for these to escape the energy 
must be transferred to the axial degree of freedom. This energy transfer process was 



analyzed theoretically by Surkov, Walraven, and Shlyapnikov |^ , and we follow their 
analysis. 
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In a harmonic trap the potential is separable, and the particle motion is com- 
pletely regular; no energy exchange occurs. In the loffe-Pritchard trap, however, 
energy exchange can occur because the potential is not separable; the radial oscilla- 
tion frequency depends on the axial coordinate, z, and so radial motion can couple 



to axial motion. (See equation p]5 



This energy mixing can be understood by considering how rapidly the radial oscil- 
lation frequency changes as an atom moves along the z axis. If the frequency changes 
slowly ( "adiabatically" ) , then the energy will not mix among the degrees of freedom. 
The adiabaticity parameter is the fractional change of the radial oscillation frequency 
in one oscillation period as the atom moves axially through the trap. Strong mixing 
occurs when 

% ~ 1. (2.6) 

Here (Op = {dujp/ dz){dz/ dt). For a loffe-Pritchard trap with a bias that is large 

compared to /c^T, uOp = a/ \J [f3z'^ + 6)m. We have used the expansion of Vip from 
equation ^.2| , which is valid if ksT ap <^ 9. Given that ksT <^ 6, the adiabaticity 
parameter is 



Up Pzy/m 



(2.7) 



We see that several factors contribute to good mixing: large axial velocity Vz (which 
occurs at high temperature), small radial gradient a, small bias field 6*, and large 
axial curvature (5. In practice, however, achieving BEG requires low temperatures 
{vz small) and high densities (obtained with large compressions, and thus large a). 
Gonsequently, the degrees of freedom do not mix and evaporation becomes essen- 
tially one dimensional. Typical values for our experiment are a/ks = 16 mK/cm, 
13/kB = 25 /iK/cm^, O/ks = 30 /xK, T = 100 /xK, z ~ 2 cm, and Vz = 140 cm/s, so 
that u)p/ujp ~ 10~^. For these conditions it takes about 10'^ oscillations to transfer 
energy, but there are only about Up/27iTcoi ~ 200 oscillations per collision for a peak 
sample density uq = 10^^ cm~^ (Tcoi is the elastic collision rate, p |121| ). There is not 
enough time to transfer the radial energy to axial energy before the particle has a col- 
lision. The energy mixing is thus very weak and the evaporation is one dimensional. 
Surkov et al. |^ pointed out the consequences of one- dimensional evaporation. They 
estimated that the efficiency is reduced by a factor of Arj compared to that of full 3D 
evaporation. The evaporative cooling power thus drops dramatically. Experiments 
in our laboratory (unpublished) have confirmed that phase space compression ceases 
near 100 yuK. These results were duplicated and studied in more depth by Pinkse et 

al. g. 

In order to maintain the evaporation efficiency a technique is required that quickly 
removes all particles with energy greater than the trap depth. To this end we imple- 
mented rf evaporation, as described in detail in chapter H. 



2.2 Degenerate Bose Gas 
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2.2.1 Bose Distribution 



The statistical mechanical description of a classical gas, presented in appendix ^ 
is not correct when quantum effects are important. A simple way to understand 
the crossover to the quantum regime is to recall that particles are characterized by 
wavepackets whose size is related to their momentum by the Heisenberg momentum- 
position uncertainty relation. As a gas is cooled, the particle momenta decrease, 
and the wavepackets enlarge. The classical (point-particle) description of the system 
breaks down when these wavepackets begin to overlap. The quantum treatment cor- 
rectly deals with the effect of particle indistinguishability. There are many excellent 
treatments of quantum statistical mechanics |3^. Here we review the basic features 
that are important for our experiment. 

A collection of identical particles with integer (half-integer) spins must have a total 
wavefunction that is symmetric (anti-symmetric) when two particles are exchanged. 
The connection between spin and symmetry has been explained using relativistic 
quantum field theory [Q. Particles with integer (half- integer) spin are called bosons 



(fermions). In contrast to fermions, multiple bosons may occupy a single quantum 
state. 

The occupation function for a gas of identical bosons in a box of volume V, and 
in the limit of A^ ^ oo and V —>■ oo but constant N/V, is called the Bose-Einstein 
occupation function [BH] 

1 

^^^(^) = exp(-M)ABT_i' (2-8) 

where /i and T are Lagrange multipliers which constrain the system to exhibit the 
correct population and total energy through the conditions 

^ ^de^fBEie) (2.9) 



V J V 

and 



^ - fdee^fsEie). (2.10) 



V J V 

The physical interpretation of these parameters is that /i is the chemical potential 



and T is the temperature. Figure shows the energy distribution of the population 
in the trap for a classical gas described by the Maxwell-Boltzmann (MB) distribution 
and for a quantum gas described by the Bose-Einstein (BE) distribution. The two 
functions shown describe systems with identical peak densities, a parameter relevant 
to the creation of a condensate. The MB distribution requires many more particles 
than the BE to create the same peak density because the atoms are distributed at 
higher energies, further from the center of the trap. 

The difference between the BE and MB distributions originates in the different 
assumptions about particle distinguishability. To intuitively understand this, we con- 
sider the outcome of these differing assumptions for sparsely and densely occupied 
sets of energy states. In the classical treatment there are many ways to arrange the 
distinguishable particles among sparsely occupied energy states because interchang- 
ing any pair of particles leads to a different arrangement. For very densely occupied 
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Figure 2-1: Comparison of Maxwell-Boltzmann (dashed line) and Bose- Einstein 
(solid line) distributions of particle energies, plotted for various chemical potentials. 

The peak density, a quantity related to the onset of EEC and to evaporation and decay 
rates, is identical for each pair of curves. The vertical axis is the number of atoms with 
the given energy. This is proportional to the occupation function weighted by the 
total energy density of states for a loffe-Pritchard trap, oc €^ + 29e^. Here 9/kBT — 1. 

states, interchanges of distinguishable particles often result in the same arrangement 
because the particles often do not change energy level. For a collection of distinguish- 
able particles with a given total energy, there are "extra" arrangements in which 
the particles are sparsely distributed. In contrast, for indistinguishable particles any 
interchange of two particles results in the same overall arrangement of particles, re- 
gardless of whether the states are sparsely or densely occupied. There are thus no 
"extra" sparse arrangements of particles. 

In a gas near the quantum degenerate regime the occupation of the energy levels 
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ranges from dense (~ 1 at low energies) to sparse (<^ 1 at high energies). Using 
the assumption of equal a priori probabilities , distinguishable particles will most 



likely be found in those arrangements with the more sparse level occupations (higher 
overall energies) because there are so many more of these arrangements. Indistin- 
guishable particles give no rise to such "extra" arrangements, and so the most likely 
arrangements will involve more dense level occupations (lower overall energies) than 
predicted by classical theory. Figure |2-1| demonstrates this effect. For each pair of 



curves the peak density, no, is identical. The peak density is the experimentally ob- 
servable quantity. Since the MB occupation function favors arrangements with higher 
energy (and thus atoms are distributed over a larger volume in the trap), more atoms 
are required in the sample to produce the same peak density. The curves thus do not 
exhibit equal area. 

Bose-Einstein condensation occurs when the chemical potential goes to zero and 
the occupation of the lowest energy state diverges. This occurs at the critical density 



n, = g■,/^{l)/^\T) (2.11) 



where gn{,z) = Y^iZi z'' /l"", (73/2(1) = 2.612. A gas that has undergone the Bose- 
Einstein phase transition is said to be in the quantum degenerate regime because a 
macroscopic fraction of the particles are in an identical quantum state. 

Although a hydrogen atom consists of two fermions, it behaves like a composite 
boson for the studies in this thesis because the collision interaction energies are ex- 
tremely small compared to the electron-proton binding energy . The two fermions 



act as a unit except in high energy collisions where electron exchange is possible. The 
typical interaction energy during low temperature collisions is ~ 1 mK, which corre- 
sponds to 10~^ eV, 10^ times smaller than the binding energy. 

To analyze the behavior of the degenerate Bose gas we shall separately treat the 
condensed and non-condensed portions of the system. The condensate is treated 



in section |2.2.2| . The non-condensed fraction, which we call the "normal gas" or 
"thermal gas", is treated in appendix 0. The results derived here will be used to 
understand how the hydrogen system is different from other gases that have been 
Bose condensed (section p73|) , and to interpret the experimental results in chapter 
^ (note that truncation effects are important because of the shallow traps used for 
those experiments, 4 < < 6). 

2.2.2 Description of the Condensate 
2.2.2.1 Gross-Pitaevskii Equation 

When Bose-Einstein condensation occurs, a macroscopic fraction of the particles oc- 
cupy the lowest energy quantum state of the system, and thus have the same wave- 
function. For a non-interacting Bose gas, that wavefunction is simply the lowest 
harmonic oscillator wavefunction for the trap. Interactions become important when 
many particles occupy this region of space and the local density increases. In this 
case the wavefunction spreads out due to repulsion among the atoms. There is a large 



body of literature on Bose-Einstein condensates (see for a recent review), and so 
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we simply quote here the most important results. 

The Schroedinger equation for the interacting condensate is called the Gross- 



Pitaevskii equation and has the form 

- I^^V^^W + V{r)ij{r) + Uo mr)f ^(r) = /i^(r) (2.12) 

where ip{r) is the condensate wavefunction to be determined. The eigenenergy of the 
wavefunction is /i, which is the total energy of each condensate atom. The quantity 
Uq = Anfi^a/m parameterizes the mean field energy, which is the energy of interaction 
among the atoms per unit density, and is repulsive for s-wave scattering lengths a > 0. 
For hydrogen in the ground state, a = 0.648 A|Q, and Uo/kB = 3.92 x 10~^^ /xK cm^. 



The mean field energy augments the trap potential by an amount proportional to 
the local condensate density ncond(r) = |V'(r)|^- Note that interactions between the 
condensate and non-condensed atoms are neglected here. For the experiments in this 
thesis, this interaction energy is small because the density of non-condensed atoms 
(called the "thermal gas" or "normal gas") is small. Furthermore, the interaction 
energy with the thermal gas is essentially uniform across the condensate because 
the density of non-condensed atoms varies only weakly; the thermal energy of the 
normal gas [T^ ~ 60 /xK) is much larger than the condensate mean field energy 
[fi/ks ~ 2 fiK), and thus the density of the normal gas does not change appreciably 
on the length scale of the condensate. 

We identify the eigenenergy n in equation |2.12| with the chemical potential of the 



system in equilibrium. The chemical potential is the energy required to add a particle 
to the system. When a condensate is present, the normal gas is "saturated", and any 
particles added to the system go into the condensate. The energy required to add the 
last atom to the condensate is /i, the eigenenergy, and so we link /i with the chemical 
potential. In practice we measure fi spectroscopically through the peak density at 
the center of the condensate, which is in turn measured through the cold-collision 



frequency shift. This is explained in section p. 3.1 



2.2.2.2 Thomas- Fermi Approximation 

If the condensate density is large enough so that the mean field energy is much greater 
than the kinetic energy of the wavefunction, then the kinetic energy term (also called 
the quantum pressure term) in equation |2.12| may be neglected. Using this "Thomas- 



Fermi" approximation we obtain the condensate density profile 

ricondir) = np- V{r)/Uo (2.13) 

where Up = fi/Uo is the peak condensate density (the largest density in the condensate, 
found at the bottom of the trap). The Thomas- Fermi approximation is valid over 
most of the volume of the condensate, but near the edges the condensate density 
approaches zero, and the kinetic energy term should be included. We ignore this 
small correction in the calculations which follow because it is minor for interpretation 



of the experiments described in this thesis. See for a good description of various 
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ways to go beyond the Thomas- Fermi approximation. 

The condensate density profile may be obtained without the Gross-Pitaevskii 
equation by assuming the condensate is stationary and its particles are at rest, and 
then balancing hydrodynamic forces P2[. A condensate particle in a region of poten- 



tial energy e has total energy E = e + Ucondi^Wo. Since there must be no net force 
on the particle, E = const = UpUo and Ucondi^) = np — e/Uq. 

2.2.2.3 Population and Loss Rate 

The population of the condensate is easily computed using the Thomas-Fermi wave- 
function in the bottom of the IP trap, which is parabolic for the condensate sizes 
and trap parameters of interest in this thesis. We approximate the potential energy 



density of states (see p |114D as g{e) ~ AIp^/66 and obtain 

cUoUp 



= / de Aipd^/£{n„- s/Uq) 
Jo 



= ±Ajp9U'J\l/\ (2.14) 

As noted above, the total energy of each condensate atom is the chemical potential 
/i = UpUo, so Ec = fiNc. 

The loss rate from the condensate due to two-body dipolar relaxation is found by 
integrating the square of the density over the volume of the condensate: 

N2,c = - deAipe^s{g/2\){np-e/U^f 

•J 

= -}^9a^q^V2.,2 (2.15) 
105 2! ° P ^ ' 



where the 2! accounts for correlation properties of the condensate |^3|, Here the 
dipolar decay rate constant g is slightly different from that given in equation |A.19 
because the rate constant Gdd^cd should be multiplied by 2; the energy 9 liberated 
in this process is large compared to fi. This term is small, so g is not effected much. 
The energy loss rate is £'2,c = /^-^2,c- The characteristic condensate decay rate in a 
parabolic trap is 

2 

- ^ = ijgnp, (2.16) 
which is 1.6 s^^ for a typical peak density = 5 x 10^^ cm~'^. 



2.3 Properties of a Bose-Condensed Gas of Hydro- 
gen 

Hydrogen differs in several ways from alkali metal atoms that have been Bose con- 
densed. The principal differences are its small mass and small s-wave scattering 
length. How do these properties influence the system? 
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First, the small mass implies that BEC occurs at a higher temperature for a given 
peak density: from equation oc 1/m. We have observed the transition at 

temperatures roughly 50 times higher than in the other systems. 

Further, as estimated by Hijmans et al. and will be shown in the following 
sections, the maximum equilibrium condensate fraction is small for hydrogen. This 
will be explained by noting the relatively high density of the condensate, as compared 
to the thermal cloud. This high density leads to high losses through dipolar decay, 
which result in heating of the system. This heating must be balanced by evaporative 
cooling, which proceeds slowly in hydrogen because of the small collision cross section. 
The result of these factors is that only small condensate fractions are possible in 
equilibrium. Possible remedies will be noted. 

Finally, hydrogen's small collision cross-section should allow condensates of H to 
be produced by evaporative cooling that contain many orders of magnitude more 
atoms than those possible in alkali-metal species, as explained in the section |2.3.4. 



2.3.1 Relative Condensate Density 

For trapped hydrogen, the condensate density grows much greater than the density 
of the non-condensed portion of the gas at even very small condensate fractions [|18] 



a distinct difference from other species that have been Bose-condensed. This is note- 
worthy for possible future studies of behavior near the phase transition. As shown in 
section |2.3.2| , it also has implications for the maximum condensate fraction that may 
be achieved in hydrogen. In this section we calculate the ratio of the densities as a 
function of the condensate fraction. We compare hydrogen to Li, Na, and Rb. 
The condensate fraction is 



where the population of the thermal cloud, Nt-, is found using equation |C.7| and the 



population of the condensate is given by equation 2.14| . The fraction can be written 



in terms of the more convenient population ratio, / = N^/Nt, as F = //(I -|- /). For 
a given peak condensate density, sample temperature, and set of trap parameters, the 
population ratio is 

_ / 2^^W^\ / \ / n^^ \ 

^ v v V ' 

/o ftrapi'P) f sample 

where /o = 7.37 x 10^'^'' (cm^)^/^yuK^ for H; ftrap indicates whether the trap shape 
experienced by the thermal cloud is predominantly harmonic (large (p) or predomi- 
nantly linear (small 0) in the radial direction {(j) = O/ksT is a unitless measure of the 
trap bias energy); and f sample carries the details of the thermal cloud and condensate 
occupation. 

The ratio of the peak condensate density to the critical BEC density, R = 
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species 


m 


a 


Tc 


Ratom 





ftrap{4>) 




B 
Li 

Na 
Rb 


[ 

i 

47 
44 


1 


1 
7 

23 
87 


0.648 
-14.4 
27.5 [ 
57.1 [ 


4U| 
46 1 
48|] 


60 
0.30 

2.0 
0.67 


230 
98 
26 
16 


0.6 

2.2 X 10^ 
30 
160 


0.0043 
0.012 
0.013 
0.011 


94 
25 

6.8 
4.4 



Table 2.1: Comparison of the parameters governing the ratio of condensate density 
to thermal gas density for the various BEC experiments. The units are amu, 
The density ratio (for small condensate fractions) is the last column multiplied by the 
condensate fraction to the 2/5 power. This ratio is much larger for H than for Na or 
Rb, indicating that two-body loss from the condensate becomes important at much 
lower /. For Li, hydrodynamic collapse of the condensate occurs before the two-body 
loss process studied here becomes important. 



np/nc{Tc), can be expressed in terms of the occupation ratio as 



3/10 / r A /■ iN\ 2/5 



~ U29/10^3/2(l)7rVl0j [mkBT^a') [ ftraM ) ^^'^^^ 



i?0 R, 



atom 



The prefactor is Rq = 0.0457364. Table p.l| gives the parameters appearing in equa- 
tion p.l9| for the various BEC experiments. For even a small occupation ratio of 
/ = 0.05, the H condensate will be 28 times more dense than the thermal gas, which 
can be compared to a Rb condensate that will be only be 1.3 times more dense. The 
loss rates from the condensate are thus very different for the two systems, as shown 
in the next section. Note that the trap oscillation frequencies play no role in these re- 
sults. The only assumptions are that the trap is of the IP form with bias fieldQ 0, that 
the condensate is in the Thomas-Fermi regime, and that fi <^ ksT so that mean-field 
interaction energy of the condensate with the thermal cloud can be neglected. 

2.3.2 Achievable Condensate Fractions 

The temperature of a trapped sample, and thus its condensate fraction, is set by a 
competition between heating and cooling, as described in section |A.2.1| . Hydrogen's 



anomalously low scattering length means that the elastic collision rate is small, and 
thus evaporation proceeds slowly. The evaporative cooling power is small, and the 
heating-cooling balance favors higher temperatures and lower condensate fractions 
than would occur if hydrogen had a larger collision cross section. Higher temperatures 
are also favored by the large density of the condensate, which places an extra heat 
load on the system due to dipolar relaxation. This extra heat load can easily be much 
larger than the dipolar relaxation heat load of the thermal gas. In this section we 



^When (j> is small, the condensate experiences a harmonic potential while the thermal gas expe- 
riences a predominantly linear potential in the radial direction. When (f) is large, both condensate 
and thermal gas experience the same potential functional form. 
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Figure 2-2: Diagram of system consisting of thermal cloud and condensate. We 
assume the particle and energy exchange rates between the condensate and thermal 
cloud are fast compared to the loss rates. The dynamics of the gas are set by evap- 
oration and dipolar decay from the normal gas when the condensate is small. When 
the condensate is large, however, N2,c 3> -/V2,t, and the dynamics are dominated by 
the dipolar decay rate from the condensate in balance with the evaporation rate. 



examine each of these factors in more detail. 

We model the system as two components, the condensate and thermal gas, that 
are strongly linked. As indicated in figure |2-2| , losses from the system occur through 
dipolar decay in the thermal gas, dipolar decay in the condensate, and evaporation 
from the thermal cloud. We assume that particles and energy are exchanged between 
the condensate and thermal gas quickly compared to all loss rates. The system is in 
a dynamic equilibrium. In section |2.3.3| we examine the validity of the assumption of 
fast feeding of the condensate. 



2.3.2.1 Dipolar Heating Rate 

Here we study the rate of heating the system due to dipolar decay in the condensate 
and in the thermal gas. We define heating as the removal of atoms from the sample 
which carry less than the average amount of energy per particle. We find that for 
hydrogen the condensate heat load exceeds the thermal gas heat load at condensate 
fractions of 0.3%. We then compare this ratio for H to that of the atomic species in 
other BEC experiments. While illustrative, this comparison is not completely appro- 
priate because the loss mechanisms in the other systems are different. Hydrodynamic 
collapse of the Li condensate prevents it from growing larger than about 1000 atoms 



in the experiments of Hulet et al. R 50 1. Three-body loss process are dominant over 



two-body dipolar decay in the Rb and Na experiments. The comparison made in this 
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section therefore serves simply to indicate how different hydrogen is from the other 
species. 

The heating rate for a process (labeled a), with energy and particle loss rates E^j 
and Nrj, is the difference between the average energy per particle in the system and 
the energy per particle that is removed, multiplied by the rate at which particles are 
removed: 

H, ^ iV. (I - I) (2.20) 

where E/N is the average energy per particle for the whole system. In this analysis 
N = Nc + Nt. The characteristic rate for this process is H„/E. 

We wish to find the ratio of the heating rates due to dipolar decay from the 
condensate (He) and from the thermal cloud (Ht) as a function of the occupation 
fraction /. Losses from the condensate play a significant role in the trap dynamics 
when Hc> Ht. 



By inverting equation |2.14| to get Up, and using the definition of /, we obtain A^2,c 
(from equation |2.15| ) in terms of / and Nt. For a small condensate (i.e. fi -C Et), 
straightforward algebra obtains 

I 7r2/5 (vr^ + 6OC(3)0)^/5v4o(r7, 0)^/^ 




221/10 105 (f)2/5Q^ (A „^ (l _ Q2{c^,v)Ao{v,m+f) 



(2.21) 



where B{(j)) was defined in equation |C.12| . The function C{(f),r],f) has only weak 



dependence on / for / <^ 1, so the dominant dependence of the heating ratio on the 
occupation fraction is /^/^. The function C(0, ?7, / = 0) is plotted in figure E^. 



For the experiments described in this thesis ~ 0.6 and ~ 5. Then C ~ 17, 
and for small / 

^ = /7/5 3,9 X lOl (2.22) 

The heating rates are equal when / = 0.27%. This means that heating of the sample 
due to losses from the condensate becomes a very significant problem even when 
the condensate is still quite small. It is difficult to hold the system in equilibrium 
for large condensates because the heat load that must be balanced by evaporation 
quickly grows too large. To obtain a large condensate fraction, the normal gas could, 
of course, be removed, and then E/N would be fi. The heating rate would be zero. 
However, the condensate would decay quickly and there would be no reservoir from 
which to replenish it. 

The experiments on hydrogen may be compared to other BEC experiments using 
Na, Rb, and Li. Relevant values for the quantities in equation p.21| are tabulated in 



table p.2| . These values help us understand why high equilibrium condensate fractions 
can be created in Na and Rb systems. Note that three-body loss processes, while 
insignificant in hydrogen, are the limiting factor which precludes larger condensate 
densities in Na and Rb. 
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Figure 2-3: Behavior of function C{(j),ri, / = 0), plotted for various trap bias fields 
(j) = O/ksT. The asymptotic behavior for large can be seen. 



2.3.2.2 Evaporative Cooling Rate 

In this section we compare the evaporative cooling rates of hydrogen and other species. 
We use expressions derived using the MB energy distribution for the thermal gas 
because analytic expressions exist and should yield rates reasonably close to those 
obtained using the BE distribution and the quantum Boltzmann transport equation. 
We postulate this because evaporation is driven preferentially by atoms with high 
energy, which are in regions where the energy distribution is essentially classical. 
Furthermore, we examine here only the scaling of the evaporation rate with the atomic 
properties m, a, and the experimentally chosen temperature T; exact rates are not 
of interest. The correct particle loss rate may be smaller by a factor of order unity 
from the results presented here, but the energy carried away by evaporating atoms 
should be very similar for the two distributions. The escaping atoms must have energy 
greater than the trap depth whether the gas is classical or quantum. 

As shown in section |A.2.1| , the energy loss rate due to evaporation is 




evap 



(2.23) 
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species 




C(0,r7 = oo,/ = O) 


/= 


H 
Li 

Na 
Rb 


[ 

i 

47 

M 


1000 

130 
70,000 
130,000 


16 
7.2 
7.4 
8.3 


0.28 % 
0.91 % 
2.3 % 
3.1 % 



Table 2.2: Comparison of the parameters entering the heating/coohng balance for 
the various BEC experiments. The units are amu. A, /iK. The second column, 
ma^Tj?, is proportional to the cooling rate. The third column carries the trap shape 
dependence of the dipolar decay process in the thermal gas. The forth column, 
/= = {RatomC{(j), oo, 0))"^/"^, is the occupation fraction at which the heating rate due 
to dipolar decay in the thermal cloud matches that due to the condensate. 



for ?7 > 6. The cooling rate is 

Hevap -^evap ( ~7Z "77 j (2.24) 

\-^^ evap / 

which, for large rj, is NeyapknT times a factor which depends linearly on rj. To 
compare the cooling rates of gases of different species, it is thus sufficient to consider 
Nevap/N ~ noave~^{r] + const). When uq = ndTc), the ratio is, for fixed rj, 



Nevap {mkBT.f^ 2 ^^bT, 

OC — 87ra ■ ' 



oc ma^T'^ (2.25) 



which is listed in table |2.2| . We see that, for a similar r/, the cooling rate for H is much 
smaller than that for Na and Rb. To increase the cooling rate in H, the parameter rj 
must be reduced, which reduces the cooling efficiency. 

We conclude that the equilibrium condensate fraction expected in H is much 
smaller than that in Na and Rb for two reasons. First, loss from the condensate 
significantly infiuences the thermodynamics of the system at much lower condensate 
fractions in H experiments since /= is so small. Furthermore, the cooling rate for a 
given 7] is much less in the H experiment. It is thus easier to achieve large condensate 
fractions in Rb and Na than in H. Note that the condensate fraction for Li experiments 



is limited by hydrodynamic collapse of the condensate |0 



A possible remedy for the low evaporation rate in H experiments has been sug- 



gested by Kleppner et al. [0]. A small admixture of impurity atoms in the trap could 
act as "collision moderators" since the collision cross section with H would most likely 
not be anomalously low; it should be on the order of 10^ times larger than that of H-H 
collisions. The increased evaporation rate should allow more atoms to be condensed, 
and with a much faster experimental cycle time. 
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2.3.3 Condensate Feeding Rate 



As the condensate decays by dipolar relaxation it must be fed from the thermal cloud 
(which is assumed to be at the critical BEC density, ric). This feeding process may 
be bottlenecked by the small collision cross section, cr, of hydrogen. To estimate 
the maximum condensate replenishment rate, Rmax, we consider the event rate for 
collisions between two thermal atoms assuming that the collisions occur inside the 
condensate. Ignoring stimulated scattering for the moment, this rate is the collision 
rate per atom times the number of thermal atoms in the region of the condensate, 

Ncoljeed = {nc(TVBEV2) (Vcond ric) (2.26) 

where 



= ,!/!(eXj) (2.27) 

is the mean particle speed in a Bose gas (ignoring truncation; g2{l)/ g3/2{i) = 0.6297). 
The condensate volume (for samples in the Thomas-Fermi regime) is, using equation 



2.13 



V,ond = 2Aipfi'rU^ + ^-]. (2.28) 




For small chemical potentials we approximate Vcond = 2Aip9n^^'^ /3. We assume that 
some fraction ( of the collisions result in an atom being added to the condensate. 
One might expect ( to be less than unity because only a fraction of collisions involve 
atoms with initial momentum and energy consistent with one atom going into the 
BEC wavefunction, which has nearly zero momentum and energy. On the other 
hand, the Bose enhancement factor very strongly favors population of the condensate. 
Calculations by Jaksch et al.^^ (equation 19a) indicate that for a condensate near 
its equilibrium population, and with /i <^ ksT, the fraction is C ~ 1/5- For the 
moment we leave the parameter free. The maximum condensate replenishment rate 
is (Ncoi ,feed- 

In equilibrium the condensate population is varying slowly compared to the feeding 
and loss rates. We may therefore equate the maximum replenishment rate to the 
condensate dipolar decay loss rate (from equation |2.15D : 

^.liAjpeU^nr = ^-^CnlavBEAjpUrnr (2-29) 

We find that the loss rate matches the maximum feeding rate when the ratio of peak 
condensate density to thermal density is 

1/4 



^Sl, . (2.30) 

nc V 9VT^ \ m J 

For hydrogen at T = 60 /xK and C = 1 the ratio is Up/uc = 16. The ratio is strongly 
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influenced by the scattering length, and is about 30 times larger for the Rb experiment 
listed in table p.2| , assuming a decay rate g = 10~^^ cm^/s [^. In the experiments 
described in chapter |^ the largest observed value of this ratio is about 20. 

A detailed understanding of the maximum condensate feeding rate is clearly de- 
sirable, but is beyond the scope of this thesis. Nevertheless, it is clear that the finite 
rate of replenishing the condensate from the thermal cloud can bottleneck growth of 
the condensate. 



2.3.4 Ultimate Condensate Population 

One figure of merit for a Bose condensate is the total number of condensed atoms, 
which is related to the number of atoms in the trap just prior to condensation. We 
investigate here one limit to this quantity. 

For a trap of effective length L and effective radius r which contains a density n, 
the atomic population is 

~ Tir^Ln. (2.31) 

If evaporative cooling is used for phase space compression, one must consider that for 
efficient cooling to occur, energetic atoms must be able to reach a pumping surface^ 
before having a collision. For a cigar-shaped trap this sets a maximum radius given 
by the collision length 

fmax = ^ = . (2.32) 

na 

For a density-limited sample size the maximum trap population is 

N^ax = ^ (2.33) 

which exhibits a strong dependence on the scattering cross section. Spin-polarized 
atomic hydrogen features a cross section more than 10^ times smaller than that of 
the alkalis, indicating that samples with 10^ times more atoms are possible. Of 
course, the alkali traps can work at lower densities, but this necessitates slower forced 
evaporation, and hence increases the loss due to background gas collisions. Use of a 
"pancake" instead of "cigar" geometry would mean the sample size is only limited 
in one dimension instead of two. However, the dimensionality of evaporation would 
also be reduced from two to one dimensions, and evaporation would thus become 
significantly less efficient. 

Table |2.3| summarizes this fundamental limit on condensate population for thermal 
cloud densities of n ~ 1 x 10^^ cm~^. Perhaps more important that producing huge 
condensates is the sustained production rate of cold, coherent atoms for an atom 
laser. This production rate should scale as Nmax times the evaporation rate, which 
scales as ma?T'^. As shown in table the evaporation rate for Na and Rb is about 
10^ higher than for H. We therefore expect the production rate of coherent atoms to 



pumping surface is basically a one-way valve: atoms which reach the surface leave the trap 
and never return. 
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species 


a 




■^max / L 




2 

cm 


cm 


cm"! 


H 


1.1 X 10-15 


10 


3 X 10^6 


Na 


2 X 10-12 


5 X 10-3 


9 X 10^ 


Rb 


8 X 10-12 


1.2 X 10-3 


5 X 10^ 



Table 2.3: Limits on number of trapped atoms when evaporative cooling is the sole 
coohng mechanism. We assume a peak density no 1 x 10^^ cm"3. 

be 10^-2 — 10^ times bigger for H than for Na and Rb. 
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Chapter 3 

Implementing RF Evaporation 



Previous attempts to realize BEC in hydrogen were thwarted by inefficient evapora- 
tion (see section |2.1| ). This chapter discusses the solution to this problem and the 
hardware required to implement this solution. 

In the improved apparatus an rf magnetic field couples the Zeeman sublevel of the 
trapped atoms to an untrapped level, but only in a thin shell around the trap called 
the resonance region. Resonance occurs where the energy of the rf photons, hurf, 
matches the energy splitting between hyperfine states, proportional to the strength 
of the trapping field. Only those atoms with enough energy to reach the high potential 
energy of the resonance region are ejected. The resonance region thus constitutes a 
"surface of death" which surrounds the sample, and so atoms with high energy in any 
direction are able to quickly escape. 

The rf evaporation scheme needed to be incorporated into an existing cryogenic 
experiment. The basic geometry of the trapping cell is constrained by the supercon- 
ducting magnets which create the trap potential. A vertical bore of 5 cm diameter is 
available for the cell along the cylindrical symmetry axis of the magnets, as shown in 



figure 1-2. Into this bore fits the H2 dissociator and the plastic tube which contains 



the gas as it is being trapped. The length of the trapping region is about 60 cm. 
The cell is thermally and mechanically anchored to the dissociator, which is anchored 
to the mixing chamber of an Oxford Model 2000 dilution refrigerator. The mixing 
chamber is situated above the magnets. 

A complete redesign of the cryogenic trapping cell was required in order to incor- 
porate the rf magnetic field. Any good electrical conductors in the vicinity of the rf 
fields and thermally connected to the cell would absorb power from the rf field via 
rf eddy currents, causing the cell to heat up. Metals, some of which had supplied 
thermal conductivity, thus needed to be eliminated. A superfiuid helium jacket was 
employed to provide the heat transport formerly supplied by these metals. In addi- 
tion, the coils which drive the field were compensated so that the field is very weak 
far away from the atoms, where the large pieces of copper that comprise the disso- 
ciator and mixing chamber are located. In this chapter we discuss these and other 
design considerations, provide construction notes, and describe performance tests of 
the apparatus. 
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3.1 Magnetic Hyperfine Resonance 



In this section we determine the magnitude and frequency of the rf magnetic field 
required to drive the evaporation efficiently. We summarize the ground states of 
hydrogen, beginning with the hyperfine interaction, and then adding the Zeeman 
effect. We obtain the four ground states and their energies. We then derive the 
matrix elements for rf transitions between these states. After making a series of 
simplifying approximations, we obtain the Rabi frequency for these transitions. An 
analytic theory for transition probabilities in iV-level systems is then adapted to our 
situation. We then calculate the field strength required to drive evaporation efficiently 
for the trap and sample parameters of immediate interest in this experiment. 



3.1.1 H in a Static Magnetic Field 

The ground states of hydrogen in a magnetic field are influenced by the hyperfine 
interaction (interaction between electron and proton magnetic moments) and the 
Zeeman interaction (the proton and electron magnetic moments interacting with the 
applied magnetic field). Here we summarize Cohen- Tannoudji et al. |Q (we use SI 
units). 

We start in the zero magnetic field regime. The hyperfine interaction Hamiltonian 

is 

nhf = Al-S (3.1) 

where I is the nuclear spin operator, S is the electron spin operator, and Ah = 27r x Vhf. 
The hydrogen ground state hyperfine frequency isvhf = 1.420 GHz. The total angular 
momentum is F = I + S. The eigenstates \F,mF) can be written in the \ms,mj) 
basis: 

(3.2) 
(3.3) 
(3.4) 
(3.5) 

In an applied magnetic field, B = Bqz, the Hamiltonian includes the Zeeman 
term, 

nz = -B-{Ms + Mi), (3.6) 

where = {qge/2me)S, M/ = {—qgp/2mp)I, q is the electron charge (g < 0), Qe — 
2.0023 is the electron gyromagnetic ratio, and Qp ~ 5.585 is the proton gyromagnetic 
ratio. The eigenstates of the combined Hamiltonian Ti. = Tikf + 'Hz are then 

\d) = |1,1) (3.7) 
|c) = cos^|l,0) + sin^|0,0) (3.8) 
\b) = |1,-1) (3.9) 









|i,o) = 




1-,+)) 


1,-1) = 






|0,0) = 




1+,-)) 
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magnetic field [gauss] 



Figure 3-1: Zeeman diagram of hydrogen F = 1 states in low fields. State d is 
trapped, state c is untrapped, and state b is "antitrapped" , expelled from the trapping 
region. At high fields the c state is trapped, as shown in figure |1-1| . 



|a) = sin^|l,0) -cos^|0,0) (3.10) 



X 



2 



tan^ = — ■ (3.11) 



X 



where x = —Bo/Bhf, Bhj = Ah'^ / ^BQei^ + v) = 506.65 G, /i^ = \q\h/2me is the Bohr 
magneton, and (f = gpTrie/ QefTip ~ 1.5 x 10^'^ is a measure of the relative contribution 
of the proton spin to the energy. The energies are, in order of decreasing energy. 



Ed = fiBBo{l-ip) (3.12) 

= -^{l-VT+^) (3.13) 

E, = -i^j,Boil-ip) (3.14) 

Ea = —{l + VTT^) (3.15) 



Figure p-l| shows these energy levels. Confining our attention to the the low fields 
which are of interest for trapped samples near the BEC transition and to the F = 1 
manifold, we see that atoms in the d state are trapped, those in the c state are 
untrapped, and those in the b state are antitrapped (expelled). 
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3.1.2 H in an RF Magnetic Field 





n, 


1 ^bB\ = 


(1 


Hd,a\ 


Pi, 


1 ^bBi = 


(1 


Hb,c\ 


Pi, 


1 ^bBi = 


(1 



Atoms may leave the trap when their state is switched from d to a, c, or b. We drive 
transitions between Zeeman sublevels by applying an rf magnetic field perpendicular 
to the static fieldQ, B^/ = Bix cos Urft which couples states d and b to states c and a. 
Transitions between the states are induced by the interaction Hrf = — B,./-(Ms'+M/) 
with matrix elements -ffa,/3 cos Ur ft = {a \ Tirf \I3)- The non- vanishing terms are 

p)cos9 + {l + ^)sijie (3.16) 

/?)sin^- (l + (/?)cos^ (3.17) 

p)cose -{l + ^)siYie (3.18) 

Hb,aV2/fiBBi = (l + (^)cos^ + (l-(^)sin^. (3.19) 

In the limit of low trapping field Bq, 9 —Bo/Bhf. We make the approximations 
cos6' = 1, sin6' = 0, and = 0, obtaining 

Hd,c = Hb,c = h^n (3.20) 

where = ^bBi/Ti\/2 is the Rabi frequency. We ignore the a-state since it is far 
separated in energy from the F = 1 manifold. With these approximations we have 
an effective three level system with equal energy separation hbBq. 

As the atoms move through the varying magnetic field, resonant transitions occur 
in a region of thickness Ar = hVLn/a which is given by the resonance's energy width, 
TiVLhi scaled by the local potential energy gradient, a. In the vicinity of the resonance 
region, the Zeeman sub-levels can be assumed to vary linearly in space. For an atom 
whose speed v changes only slightly as it crosses the resonance region, the energy 
levels may then be considered to vary linearly in time. The probability of changing 
Zeeman sub-levels as the atom crosses the resonance region is calculated using a 
variant of the Landau-Zener formalism |Q for avoided level crossings. Vitanov and 



Suominen have analytically calculated the dynamics for avoided level crossings 
in systems with an arbitrary number of levels. They show that, for an atom initially 
in the d state, the probability of emerging from the resonance region in the d, c, or b 
state is 

Pd = {l~P? (3.21) 
= 2(1 -p)p (3.22) 
Ph = (3.23) 

where p = 1 — exp(— (^) is the Landau-Zener transition probability for a two-state 
system, and 

C = 2-Khnyav (3.24) 
is the adiabaticity parameter for a particle of speed f in a magnetic field with mag- 



We consider the polarization of the rf field in more detail below 
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netic field gradient a/^B- For large C, the atom absorbs two rf photons while in the 
resonance region, and emerges in the h state, antitrapped. For small C (the diabatic 
case of interest here) there is only a small probability of leaving the trapped d state, 
and the vast majority of atoms that do leave this state will emerge in the c state, 
untrapped. They are pulled by gravity out of the trapping region along the trap axis. 



3.1.3 RF Field Amplitude Requirement 

In order to realize the improved evaporation efficiency we seek, the characteristic time 
required to eject an energetic atom must be less than the characteristic collision time 
of these atoms |Q. This condition sets a minimum rf field amplitude. 



Here we consider an atom traversing the trap radially, passing through the center. 
Atoms in different trajectories will be effected differently by the rf field and sample 
density. However, in order to have a significantly different orbit (while still having 
enough radial energy to reach the resonance region) , the atom would need high energy 
in both the radial and azimuthal degrees of freedom. Since this is rare (see appendix 
0), we ignore these orbits and focus on the simple case of purely radial motion. 

The characteristic time to eject an energetic atom, r^jecu is the oscillation period T 
divided by the probability of ejecting the atom per period, 4Pc- (We ignore coherence 
between resonance crossings). The atom encounters a resonance region four times per 
period since there is a resonance region on both sides of the trap, and the atom passes 
through each region twice per oscillation (assuming the turning point is outside the 
surface of death). The oscillation period of a particle of mass m moving in a linear 
potential with gradient a and with an initial energy et is T = 4:^J2met/ a. For the 
parameters of interest {et/kB ~ 500 /iK, a/ks ~ 16 mK/cm used near the BEC 
transition), we obtain a period T ~ 1 msec. Thus Teject ~ 1 msec/4Pc- 

To find Pc and the minimum rf field strength, we require that the ejection time be 
less than the characteristic collision time t^oi = 1/ riQav (A rigorous analysis includes 
an average of the quantity n{p)v{p) over the atom's trajectory, which produces a 
correction factor of about 3). For atomic hydrogen near the onset of BEC at Tc ~ 
50 /iK, Tcoi ~ 100 msec. The particles thus sample the resonance region roughly 
400 times before having a collision. The ejection probability Pc must be at least 
Pmin = T/4tcoi — 1/400. The minimum adiabaticity parameter is then (min — 10^^. 
For the gradient a and typical average velocity v ^ 100 cm/s, we obtain a minimum 
Rabi frequency Qmin = \J CivC,min/h = 2rc x 1 kHz. The minimum rf field is -Bi,mm = 

hilminV^/ f^B = 1 VliG. The thickucss of the resonance region is Ar = 30 nm, so 
the approximation of a linearly varying trapping field made above is valid. We note 
that for atoms with more complicated hyperfine structure {F > 1) the rf field must 
be larger because an initially trapped atom must be taken through one or more 
intermediate states before arriving at an untrapped state. For a given ^f-factor and 



hyperfine state the required Rabi frequency and rf field scale as ^ma^T^ independent 
of trap parameters (for linear and harmonic traps). 

A potential problem is trap loss and heating through dipolar decay as atoms in 
the c state pass through the cloud of trapped d atoms |50]. This turns out to not be 
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Figure 3-2: Static magnetic field which confines the gas, viewed in a plane which 
cuts through the center of the trap. The vertical axis of the figure is the distance 
away from the central axis (symmetry axis) of the trap, in cm. The horizontal axis is 
the long axis of the trap. The atoms are confined to the bottom of the trap near p = 
and z = 0. The contours are spaced by 100 /iK. The arrows indicate the direction of 
the field. Note that the axes are very different lengths. This trap is used to confine 
samples below T = 130 fiK. In order to eject energetic atoms, the rf magnetic field 
must be perpendicular to this trapping field. The saddlepoint in the trapping field is 
on the left. The top of the cell and dissociator are to the right. 



a problem for trapped H. The time constant for this process is Tdip = l/Gdnc where 
Gcd is the total event rate constant for dipolar decay in collisions between c and d 
atoms. For H, the total event rate constant is Gcd — 10~^^ cm^/s which is about the 
same as Gdd, the total event rate constant for dipolar decay through d-d collisions 
|57| . The density of c atoms is much less than no since they are free to fall out of the 
trap {Tfaii ~ 100 msec). We conclude that Tdip is long compared to the lifetime of the 
sample as set by d-d dipolar relaxation. Therefore, we need not be concerned with 
dipolar relaxation events between the c and d atoms. 



3.1.4 RF Coil Design 

Here we describe the coil design considerations for creating an rf magnetic field that 
will drive evaporation efficiently in the cryogenic apparatus. 

The rf magnetic field must have a component in a direction perpendicular to the 



trapping field, as discussed above. Figure 3-2 indicates the direction of the trapping 



field in a plane that cuts through the trap. This field is for a typical trap used for 
rf evaporation experiments. Along the sides of the trap the rf field must point in 
the z direction. This field is generated by coils wrapped around the the cylindrical 
confinement cell. These coils are termed the "axial coils" . 

It is important that the rf field be weak enough at the top of the cell to avoid 
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Figure 3-3: Winding pattern of axial rf coil, shown wrapped on a section of the cell 
inner G-10 tube (3.8 cm OD). The coil is fed by the two wires extending out the top 
of the section. The loops are wound as close together as possible, and the feed wires 
are laid tightly next to each other. The point z = is very near the laser beam focus, 
and also the minimum of the trapping potential. 



heating the large copper pieces which make up the cell top, discharge and mixing 
chamber (heating rates are outlined in section |3.2.1|) . The field generated by a loop 
of diameter a decays as {a/zY at large distances from the loop along the loop axis. 
However, the field can be made to decay as (a/z)^ if the loop is compensated ap- 
propriately. In the far-field the arrangement appears as a quadrupole instead of a 
dipole. The heating rate is proportional to the squared magnitude of the field, so the 
heating rate falls as {a/zY^ for a compensated loop. The geometry of the axial coils 
was chosen to provide large fields over the length of the cold samples and to produce 
fields which decay rapidly outside the trapping region. The coils are wrapped around 
the outside of the inner tube of the cell (3.8 cm diameter). They consist of three loops 
in the "forward" direction at 2; = ±4 cm, 1 loop "backward" at 2; = ±3 cm, and two 
loops "backward" at z = ±1.2 cm. The winding pattern is shown in figure |3-3| . 

Figure p-4| shows a calculation]^ of the magnitude of the rf field generated by the 
axial coils along the trap axis. Figure ^-5| shows how the field decays along the trap 
axis. The nearest big pieces of metal that must be kept cold are at z = 40 cm. 
Because the rf field is compensated well, the dominant heating of the top of the cell, 
dissociator, and mixing chamber arises from other mechanisms, described in section 



^Since the wavelength of this radiation (A = 3 m at 100 MHz) is much larger than the length 
scale of interest we disregard propagation effects and calculate the dc field. 
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Figure 3-4: Profiles of tlie rf and dc magnetic fields along the axis of the trap. The 
dashed line is the rf field magnitude, and the solid line is the trapping field for the 
trap used to achieve BEC. 



3.4.4 



The rf fields are primarily used to remove the most energetic atoms, but it is 
also sometimes desirable to remove the atoms at the very bottom of the trap where 
the trapping field is aligned along the z axis. Another set of rf coils, called the 
"transverse coils", generate a field predominantly perpendicular to z. The winding 



pattern is shown in figure 3-6 



3.2 Mechanical Design of the Cell 

The designs of the cells used previously in the experiment are described by Doyle ||r^ 
and Sandberg [|16|. The principle difference between those cells and the cell used for 



this experiment is the presence of the rf magnetic fields. Previous designs relied on 
metals to provide thermal conductance along the length of the cell, but metals must 
be excluded in the new design. They would allow rf eddy currents to flow, which 
would heat the cell and which would screen the atoms from the applied rf fields. In 
section p.2.1| we examine the problem of rf heating. In section |3.2.2| we explain the 
design of a jacket of superfiuid ^He which surrounds the cell and provides the required 
thermal conductance. 
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Figure 3-5: Profile of tlie rf magnetic field along the axis of the trap showing the 
field strength at z = 40 cm where the nearest pieces of cold metal are located. 

3.2.1 Exclusion of Good Electrical Conductors 

RF eddy currents in electrical conductors will lead to heating of the cell. This heating 
during rf evaporation must be limited. The cell must be kept at a temperature below 
~ 150 mK to ensure that walls are sticky for H atoms and to ensure that the vapor 
pressure of the superfluid ^He film, which is crucial for loading the trap, does not rise 
high enough to create a significant background gas density. A failure in either regard 
would allow high energy particles (T ~ 100 mK) to knock cold atoms out of the trap. 
A technical limitation also exists: it is desirable to analyze the trapped sample very 
soon after rf evaporation, and this analysis requires the cell to be below 90 mK. If 
the cell is heated too much during the rf evaporation process it will take too long to 
cool below this temperature. 

The typical thermal conductivity along the length of the cell is ~ 1 /iW/mK so 
that heat loads of only 50 fiW are tolerable for a mixing chamber temperature of 
100 mK. We therefore must beware of power deposited via rf eddy currents in any 
metals thermally connected to either the refrigerator or the cell. In previous designs, 
good thermal conductivity along the length of the cell (a tube about 4 cm in diameter 
and 60 cm long) was supplied by about 100 Cu wires [|l^. These wires, and many 
other materials, needed to be replaced in the new design in order to eliminate rf eddy 
current heating. 

To justify the effort required to remove good electrical conductors from the cell, we 
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3.7 cm 

Figure 3-6: Winding pattern of transverse rf coils. The circular ends wrap around 
the cell, and the straight sections run parallel to the symmetry axis of the cell. The 
coil is centered at z = 0. 



here estimate the heat per unit length deposited in a wire by eddy currents. Consider 
a cylindrical conductor of radius a, length L, and electrical conductivity a immersed 
in a uniform magnetic field S, aligned along the axis of the wire, that is oscillating at 
angular frequency a;, as shown in figure |3-7|. For a good conductor the skin depth is 



roughly 8 = ^j2/afiuj [|58| where /i is the permeability; for the non- magnetic materials 
of interest here, we take fi = fiQ = An x 10"'' H/m, the permeability of free space. 

For frequencies larger than Us = 2/afia'^ the skin depth is smaller than the wire 
radius. For u ^ Ug we may assume that all the eddy current flows in a shell of 
thickness S, and that the magnetic field is extinguished in the wire everywhere except 
in this shell. We thus simplify the problem to a shell carrying current in the azimuthal 
direction through a cross- sectional area 6L and around a length 27ra. The conductance 
around the shell is then S = a6L/27ra. The shell is pierced by the magnetic field in 
an area A = 2Tra6. The power deposited (per unit length of wire) is P/L = SV'^/L 
where V = ABuj is the induced voltage around the shell. We obtain 

L Joji 



Typical values are i? = 10 ^ G, = 27r x 10* Hz. The Cu wire formerly used 
for thermal connection along the length of the cell has radius a = 0.5 mm and 



conductivity a ~ 10^" {n m)-^ |59[, so that P/L = 20 nW/cm. It is clear that 



for 100 wires of length L = 10 cm this heating rate is non-negligible. Since this 
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Figure 3-7: geometry for calculation of rf eddy currents in a wire of radius a and 
length L immersed in a uniform magnetic field B aligned along the axis of the wire. 
The eddy current flows around the cylindrical shell of thickness S, height L, and 
perimeter 27ra. The effective area for magnetic induction is 2TTa6. 



heating rate is so close to the limit, an alternative mechanism for heat transport was 
necessary. A further problem with the Cu wires is that they shield the interior of the 
cell from the rf field. 



3.2.2 Design of the Superfluid Jacket 

In order to obtain good thermal conductance along the length of the cell we exploited 
the extraordinary heat transport properties of superfluid helium. The new cell design 
involves two concentric G-10 tubes with a layer of helium between them, about 
2.2 mm thick. Figure 3^ is a cross-section through the top of the cell. There are 
several issues that were addressed in the design of the cell. These are discussed in 
this section. 



3.2.2.1 Thermal Conductivity 

The superfluid jacket must be thick enough to carry heat sufficiently well from the 
bottom to the top of the cell where it is transferred to the refrigerator. To calculate 
the quantity of helium necessary, we must estimate the required thermal conductance 
and understand heat flow in the superfluid. 

There are three temperature regimes which demand good conductance. During 
loading of the trap, when T = 250 mK, the conductivity must be at least S = 
0.1 mW/20 mK so that S/T^ = 0.3 W/K^. During laser spectroscopy, when T = 
140 mK, the conductivity must be at least S* = 0.2 mW/20 mK so that S/T^ = 
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Figure 3-8: Cross-section through dissociator and upper end of cell. The origin of 
the scale on the right hand side is at the center of the trap. The cell and dissociator 
are inside a brass vacuum can which keeps 4 K liquid helium out. The dissociator 
couples to the cell via a copper mating collar and an indium o-ring. A thin copper fin, 
covered with silver sinter, makes thermal contact with the superfluid helium shell that 
is confined by two concentric G-10 tubes. The window at the top of the dissociator 
provides optical access to the trapped atoms for laser diagnostics. Not shown are 
electrical feedthrough tubes which are soldered into the copper mating collar and lie 
in grooves on the outside of the dissociator. 



Z = 50cm 



4 W/K"^. Finally, during bolometric detection, when T = 80 mK, the conductivity 
must be at least S = 10 fiW/lO mK so that S/T^ = 2 W/K^ 

In metals heat is carried by electrons, which typically scatter off impurities on a 
length scale which is short compared to the dimensions of the object. In contrast to 
this effectively diffusive heat flow behavior, heat transport in superfluid is essentially 
ballistic because there is no impurity scattering. Furthermore, the heat is carried 
by phonons which have a different excitation spectrum than electrons. The notable 
implications are that the effective conductance scales with temperature as in He 
but only as in metals, and that the conductivity of a column of superfluid depends 
on the diameter of the column; phonons scatter when they hit the walls. The empirical 
conductivity relation for a tube of superfluid is n = (20 W/K^ cm^) dT^ where d 



is the diameter of the vessel and T is the temperature. In our design the outer radius 
of the inner tube is a = 1.92 cm and the inner radius of the outer tube is 6 = 2.14 cm. 
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We take d = b — a = 2.2 mm. The conductance along the cell is 



^=AT = X = L (^-2^) 

where L = 66 cm is the length of the cell. For our geometry S/T^ = 0.2 W/K^. The 
thickness of the shell causes scattering on a length scale d in the radial dimension, but 
the scattering length is significantly longer in the azimuthal direction. We thus expect 
the conductivity to be larger by a factor of two or three. Furthermore, the empirical 
formula for k above applies to surfaces from which the phonons scatter diffusely. If 
specular reflection occurs, then the effective scattering length is much longer . 



The characteristic phonon wavelength at 100 mK is 300 nm since the surfaces of 



the coated G-10 tubes exhibit specular reflection of optical wavelengths (A ~ 500 nm, 
observed by shining a lamp on the surface), the superfiuid phonons should also reflect 
specularly. As will be explained in section p. 4. 2| below, we have observed the superfiuid 
conductivity to be 5 times higher in tubes which exhibit specular optical reflection, as 
compared to tubes for which the optical reflection is diffuse. Putting all these factors 
together, we expect a conductance S/T^ = 3 W/K^. There is not much margin for 
error in the design. Measurements of the thermal conductivity are detailed in section 



3.2.2.2 Thermal Link from Superfiuid to Refrigerator 

Once heat has been transferred along the length of the cell, it must be transferred from 
the liquid into metal links which lead to the refrigerator. This process is frustrated 
by the large boundary resistance between the bulk superfiuid and the metal. This 



resistance, called the Kapitza resistance [BUI , arises in two ways from the factor of 25 



mismatch in the speed of sound in the two materials. First, for a phonon to propagate 
into the metal it must approach the surface at very nearly normal incidence; phonons 
with wavevectors outside a small cone of acceptance angles experience total internal 
reflection. Furthermore, the very different speeds of sound give rise to a significant 
impedance mismatch, and so most of the phonons within the acceptance cone will be 
reflected. The empirical Kapitza resistance equation is Rk = (0.02 m^ K'^/W) /AT^ 
5^ for a superfluid-Cu boundary of area A and temperature T near 100 mK. For 



a modest maximum resistance Rk = 10 mK/100 /iW at 140 mK a boundary area 
of at least A = 700 cm^ is required. This is far too large to fit conveniently in the 
apparatus. 

One way to create an effective surface area that is much greater than the simple 
geometric surface area is to apply sintered silver to the surface [^. The network 
of voids between the grains of the powdered metal creates a huge surface area. It is 
possible to create an effective surface area of 10^ cm^ for a coating of silver sinter on an 
area of copper that is straightforward to incorporate into the apparatus, A = 25 cm^. 
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3.2.2.3 Heat Capacity 



After loading tlie trap at T = 250 mK it is important to cool the cell quickly to below 
150 niK to achieve thermal isolation between the trapped sample and the cell walls, 
as explained in section |L^. Thus, it is important to ensure that the heat capacity of 
the superfiuid is low. 

The total heat that must be extracted from the superfiuid is the integral of the heat 
capacity between the initial and final temperatures. The heat capacity is C = CqT^ 
with Co = 8 X 10-2 J/mol p|. The quantity of heat to remove from N moles of 



liquid cooled from T2 to Ti is f/ = CqN{T2 — T^^)/4. From the 6 moles of superfiuid 
in the jacket we must thus remove 0.4 mJ. The refrigerator cooling power is about 
1 mJ/s, and so the time constant is less than a second. This is acceptable. 



3.2.2.4 Filling and Emptying the Jacket 

Since superfiuid is such a good heat conductor, we must be concerned about heat 
links from the coldest part of the refrigerator to warmer regions through the fill 
line which leads to the jacket. Heat may flow through two mechanisms. When 
a tube is filled with superfiuid the thermal conductivity is given by the empirical 
relation quoted above (p|50|). In this case the total heat conductance is proportional 
to d^/L. Smaller diameter tubes are obviously much better. To limit the heat flow 
to less than 10 /xW between the 1 K and 0.1 K stages of the refrigerator requires 
(P/L < 6 X 10"^ cm^; for L = 50 cm we obtain d < 0.3 mm. The other heat 
conduction mechanism occurs in tubes with only a superfiuid film; the film flows to 
a warm region, evaporates, and sends gas back to colder regions. As it recondenses it 
deposits heat. The upper limit on this process is given by the superfiuid film flow rate. 
A saturated film of thickness t ~ 30 nm flowing at the critical velocity Vc = 45 cm/s 
up a tube with limiting perimeter nd will evaporate, and then recondense depositing 
power Pfiux = Lvctud = (0.9 mW/cm)(i where L = 2 J/cm^ is the latent heat. To 
deposit less than 10 /iW at the cold end of the tube we require 0? < 0.1 mm. This 
is an overestimate of the power transport because the warm vapor does not make it 
all the way to the cold end of the tube before condensing. To minimize heat flow 
through the fill line we chose a segmented design; 40 cm sections of 0.2 mm ID tubing 
connect heat exchangers which are thermally anchored to progressively colder parts 
of the refrigerator. In this way heat is deposited as far as possible from the coldest 
part of the refrigerator. 

If the cell is warmed much above 4 K, the vapor pressure of the liquid helium rises 
far above atmospheric pressure. The tiny fill tube is a huge impedance to the gas 
flowing out of the jacket. Effectively we have created a big bomb. A pressure relief 
valve was developed to bleed off excess pressure in an emergency, thus protecting the 



cell from bursting. More details are given in section 3.3.4 
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G-10 tube 

superfluid helium m, 

— Z=-20cm 

mirror for 243 nm light 

MgF window 
G-10 mating ring- 
Figure 3-9: Cross-section through the bottom end of the celL The origin of the scale 
on the right hand side is at the center of the trap. The mating ring at the bottom 
joins the two concentric tubes. The MgF2 window is set against the bottom of the 
inner tube. A uv mirror retroreflects the laser beam which passes along the axis of 
the trap. The bolometer is placed about 1 cm above the window. 

3.3 Construction Details 

In this section we discuss some of the techniques that were crucial for the construction 
of the cell. 



3.3.1 Materials and Sealing Techniques 

It is a significant accomplishment to create joints between different materials that 
are leak-tight to superfluid flow at low temperatures. The cell is a composite unit 
consisting of a MgFg window, two concentric convolute wound G-10 tubes, an OFHC 
copper collar that seals the superfluid jacket at the upper end and mates the cell to 
the dissociator, two rf feedthroughs, and many low frequency electrical leads. The 
sealing techniques described here proved reliable. 

The G-10 tubes came from the factory with rough surfaces. To make them 
leak-tight and to create the specularly reflecting surface crucial for high thermal 
conductivity we coated them with an epoxy mixture of 100 parts by weight Epon- 
828 to 32 parts Jeffamine D-230 catalyst. The tubes were hung vertically 
while the epoxy cured and was then heat set at 60 C for 2 hours. The rf coils 
were wrapped onto the tube before the coating process to allow the epoxy coating 
to produce a smooth surface for specular reflection of the superfluid phonons. The 
tubes were joined at the bottom to a G-10 mating ring with Stycast 2850FT epoxy, 
cured with 24LV catalyst See figure |3-9| . The MgF2 window was also glued in 
place with Stycast 2850FT epoxy. In order to prevent the glue from contaminating 
the inside surface of the window (and blocking Lyman-a transmission) the following 
gluing procedure was used: the window was placed on a pedestal, and the outer 
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edge of the window was coated with glue. The cell was then held above the window 
and slid gently downward. If necessary, addition glue was dabbed to create a fillet 
around the bottom of the window against the G-10 tube. The cell was held vertically 
while the glue cured. Because 122 nm light must pass through this window, the 
transmission of a window glued this way was tested in the ultraviolet, using easily 
accessible 243 nm radiation. No degradation was observed. Windows glued in this 
way have been thermally cycled to 1.4 K many times without ever leaking. Thermal 
shock is a problem, though, so the windows were never brought in direct contact with 
liquid nitrogen during the leak checks. 

Dual thermometers and a heater are placed at each end of the cell. Other wires 
lead into the superfiuid jacket space and connect to pads or meshes used to produce 
the electric field required for rapid quenching of the metastable atomic 25* state. All 
of these leads are brought into the jacket space through a simple epoxy feedthrough; 
the wires pass through a 1.67 mm ID tube which is filled with 2850 epoxy for several 
centimeters. Each rf coil is driven by a twisted pair of superconducting leads (Cu clad 
NbTi wire with a single 33 fim diameter core) The twisted pair for the axial coils is 



connected to a semirigid coax line (2.16 mm OD Cu which passes through the 



Cu mating block at the top of the cell. The other end of the line is soft-soldered to 
a brass cylinder containing a hermetic SMA connector soldered in place. This 
sealing technique works reliably. The twisted pair which feeds the transverse rf coil 
passes through an epoxy feedthrough and is connected to an SMA connector in the 
vacuum space that surrounds the cell. 

The semirigid coax cables which deliver the rf power to the cell are composite units 
that optimize the compromise between low thermal conductivity and high electrical 
conductivity. Cables with high electrical and thermal resistivity (Lakeshore model 
CC-SR-10 ^8|) carry the rf power to the 1.6 K pumped ^He pot, where the cables are 
thermally anchored. The ohmic heating of these lines is too much of a heat load for 
the mixing chamber, however, so lower resistivity lines (Lakeshore type C) carry the 
rf power to the cell (which is thermally connected to the mixing chamber). To reduce 
the thermal conductivity of these lines, a 1 cm section of the high resistivity cable 
is spliced in. The thermal conductance is dramatically reduced, but the electrical 
conductance is high enough to eliminate excessive heating. 



3.3.2 Sintering 

Effective heat transfer from superfiuid to copper is difficult to achieve, as explained 
above. A layer of sintered silver was fused onto each side of the Cu fin, shown 
in figure |3-8|, to increase the surface area for heat transfer and reduce the overall 



thermal resistance. The fin is high purity copper [£9|, 0.5 mm thick, and is brazed 
|7(]| to the Cu mating collar. The sintering process consists of several steps [^, [72| . 
First, the surface of the fin is cleaned by immersion in 1 molar nitric acid for 30 s, 
then rinsed with deionized water. The surface is then plated with silver to increase 
the bonding of the sinter to the Cu fin. The surface is immersed in a plating solution 
(electroless silver cyanide [|7^]) for 15 minutes at 95 C. The surface is then rinsed 
again with deionized water, and the assembly is baked at 120 C for 35 minutes in an 
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air atmosphere. The baking was suggested to "diffusion weld" the plating to the Cu 
substrate. The Ag sintering powder (nominally 700 A particle size) is first cleaned 
by baking at 50 C for 30 min in an H2 atmosphere. The powder is then compressed 
onto the Cu fin in a specially designed curved press which matches the curved surface 
of the Cu fin. The fin was sintered sequentially in eight radial segments at a pressure 
of 5 X 10^ Pa. The press apparatus and mating collar were heated to about 100 C 
during the operation. About 7 g of powder were bonded to the fin, producing an 
estimated 10 m^ surface area. 



3.3.3 Bolometer 



The bolometer is used as a diagnostic of the trapped gas . As atoms are released 
from the trap, their flux is measured by detecting the molecular recombination heat, 
4.6 eV per event. A resistive element serves as both a thermometer and a heater; 
the electrical power required to maintain the bolometer at a constant temperature is 
recorded. The bolometer used in these experiments is similar to previous designs. It 
is constructed from a quartz plate 1 cm x 1 cm x 60 fim. The plate is suspended 
by nylon threads from a G-10 cradle which is glued into the cell. Electrical contact 
is made through superconducting filaments extracted from multifilament supercon- 
ducting wire by dissolution of its Cu matrix in nitric acid. The wires are contacted 



to evaporated Au pads on the quartz by conducting epoxy overlaid with Stycast 
1266 epoxy. The Au pads touch graphite resistive elements which were formed by 
dabbing Aerodag ||T6[ onto the surface. The differences in the present bolometer de- 
sign from previous designs arise from efforts to reduce rf eddy current heating. Large 
gold contact pads have been replaced by pads of a few mm^ surface area. A redundant 
design was used consisting of four contact pads in a row with three resistive elements 
between them. The pads are spaced very closely so that the resistive elements could 
have resistances of a few kQ without large mass. This bolometer has been reliable. 



3.3.4 Pressure Relief Valve 

The vapor pressure of the six moles of liquid helium in the superfiuid jacket of the cell 
can cause the cell to explode if the temperature is raised too high before the liquid is 
removed. This conjecture has been tested. Recall that the fill line for the jacket has 
a very low gas conductance, and so it takes many hours to remove the liquid. In the 
case of an accidental overheating, it is desirable that the liquid in the jacket have an 
escape path with a high conductance. 

We constructed a reusable pressure relief valve which provides this high conduc- 
tance escape path. The reusable valve is at the same temperature as the superfiuid 
jacket so that large diameter tubing can connect the two; there are no heat leak 
concerns which dictate small tubing. The design of the relief valve is conceptually 
divided into a pressure-sensing part and a seal-breaking part. The pressure sensor 



is a BeCu bellows tube ^Ij which lengthens as the pressure inside increases. The 
seal to be broken is a 25 /zm thick piece of brass shim stock which is soft-soldered 
to a brass flange, and is replaceable. When the seal is broken the gas escapes into 
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Figure 3-10: Cross-section through the pressure rehef valve. A VCR fitting 
brazed into the base makes a superfiuid tight vacuum connection to the superfiuid 
jacket through 2 mm ID tubes. The BeCu bellows tube fits around a Cu guide 
cylinder that prevents buckling and keeps the hole in the brass fiange aligned with 
the hypodermic needle. A thin brass sheet is soft-soldered to the brass fiange to 
provide a vacuum seal that is easy for the needle to puncture. 



the large vacuum space which surrounds the cell. The fiange is brazed to the bellows 
and moves toward a stainless steel hypodermic needle as the pressure increases (The 
bellows are pre-stretched during installation. They do not start to move until there 
is a pressure differential of about 0.5 bar). A touch sensor makes electrical contact 
with the brass shim stock before it meets the needle, alerting the operator that seal 
rupture is imminent. The relief valve is shown in figure 3-10 . The relief valve has 
been tested. It is simple and reliable. 



3.4 Measurements of Cell Properties 

Measurements have been made of several key parameters, such as thermal conduc- 
tivity, heat capacity, rf field strength, and rf heating rates. In this section we outline 
these results. 

3.4.1 Measurements of RF Field Strength 

The strength of the rf field generated by the axial and transverse coils was measured 
by inserting a pick-up coil inside the inner G-10 tube during construction of the cell. 
The absolute field strength and field profile along the cell axis were found to agree well 
with calculations for frequencies up to about 5 MHz, where the inductive impedance 
matched the 50 fl driving impedance and the response began to drop. When the cell 



56 



is cooled to 100 mK the response could be different because of the lower resistance of 
the superconducting wire which forms the coils. 

The axial field was measured with a 1.6 cm diameter pickup loop positioned in 
the center of the cell (the test was done before attaching the bottom window onto the 
cell). To match the impedance of the loop to that of coax signal cable, an op amp 
Il79| was connected directly to the loop. The output of the unity gain amplifier drove 
the signal cable through a 50 resistor. The frequency response of the amplifier 
was measured, and so the pickup loop probe has a calibrated response to 200 MHz. 
Knowledge of the geometry of the loop and of the gains in the system allow us to 
measure the field with an accuracy of about 20%. An HP4195A network analyzer 
was used both to calibrate the probe and to measure the fields. The coils were driven 
through the high thermal impedance microminiature coax cable [^] used to convey 
rf power from room temperature through the refrigerator to the cell. The response of 
the coils is shown in figure |3-11| . The structure above 50 MHz is most likely due to 
resonances in the feed line and coil structure. 

The cell resides inside a brass vacuum can in the experiment, and this conductor, 
near the coils, can influence the field; a counter current is generated which tends to 
partially cancel the field. This effect has been modeled, and the reductions shown in 
figure p-ll| are consistent. 

Note that as the apparatus is cooled to operating temperature the resistance of the 
axial (transverse) coil changes from 13 Q (26 Q) to zero. We thus expect a 30% (50%) 
increase in field at low frequencies, and the corner frequency should move higher by 
30% (50%). 



3.4.2 Thermal Conductivity 

The thermal conductivity along the length of the cell has been measured in the 
temperature range 50 mK to 450 mK. As described above, the superfiuid jacket is 
a shell between two concentric cylinders. The outer tube has an inside diameter of 
4.27 cm, and the inner tube has an outside diameter of 3.83 cm. The thickness of 
the shell is thus 2.2 mm. The vertical distance between the thermometers used for 
the following measurements is about 59 cm. The jacket space is cluttered by 5 mil 
diameter manganin wires which connect to thermometers, heaters, the bolometer, 
and plates used to apply electric fields. The thermal conductivity is reduced since 
phonons scatter off these structures. 

Thermal conductivity measurements were made by controlling the temperature 
at the top of the cell, applying heat at the bottom of the cell, and measuring the 
temperatures. The two resistance thermometers at the bottom of the cell allowed 
consistency checks against each other; one of the thermometers at the top of the 
cell verified the operation of the temperature controller, which used the other cell top 
thermometer. The thermometers consist of a RuO paste on a ceramic substrate, man- 
ufactured for surface mounting on printed circuit boards |8^. The four thermometers 
used in the cell were calibrated in a previous cooldown against another RuO resis- 
tance thermometer, which was itself calibrated by Oxford Instruments. The heaters 
at the top and bottom of the cell were each made of two 5 kfl metal film resistors. 
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Figure 3-11: Measured frequency response of axial and transverse rf coil. The solid 
line is the field when the cell is surrounded by the brass vacuum can. The dashed line 
corresponds to no conductors in the vicinity of the cell. The coil is driven through 
the coax lines in the cryostat. 



The two resistors for each heater were positioned on opposite sides of the jacket, 
but at the same vertical [z) position. The heaters were at about the same vertical 
position as the thermometers. The temperature of the top of the cell was controlled 
by a proportional-integral controller (1 s proportional time constant, 3 s integrator 
time constant). The thermal conductance at a given temperature was extracted from 
a series of measurements such as those in figure P-12| . Conductances were measured 
at several cell temperatures. They are compiled in figure |3-13| . All temperature rises 
were less than 10%. The conductances measured here are adequate for the experi- 
ments we wish to perform. 

The thermal conductance across the sinter plates, where heat is transferred from 
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Figure 3-12: Typical measurements used to extract thermal conductance of cell 
at one temperature. The temperature of the bottom (open symbols) of the cell 
increases linearly as heat is applied to the bottom. The temperature at the top (solid 
symbols) stays roughly constant since it is temperature controlled (using a different 
thermometer as the sensor). The lines are linear fits for the range of powers from 
to 2.5 /iW. The slopes of the two upper lines give the thermal resistance of the 
cell at this temperature. The resistances measured by the two thermometers are 
1.62 ± 0.04 mK//iW and 1.56 ± 0.04 mK//iW. The 1 mK offset between the two 
bottom thermometers is due to miscalibration. 



the superfluid to metal and eventually to the refrigerator, may be measured in an 
analogous way. The temperature of the dissociator was controlled with a thermometer 
and heater near the top of the dissociator (see figure |3-8|) . The thermal gradient 
measured thus includes the thermal resistance across the indium o-ring vacuum seal 
(superconductivity quenched by 4 T field) and along the length of the dissociator. 
The data from the two thermometers at the top of the cell were fit to straight lines, 
similarly to the data in figure 3-12 . Two temperatures and conductivities were thus 
obtained for each temperature control point. These two data sets are shown in figure 



3-14 



The conductance of the cell is similar to the conductance across the sinter. Any 
significant improvements to the thermal conductivity of the system would need to 
address both of these segments of the heat transport channel. 

In a separate experiment the effects of specular and diffuse phonon scattering 
were studied. A cell similar to that described above was made. The surfaces of the 
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Figure 3-13: Thermal conductance along the length of the cell. The two data sets are 
from the two thermometers at the bottom of the cell. Approaching low temperatures 
one expects the conductance to fall like T^, as indicated by the line. The inset 
shows the conductance with this dependence factored out. The error bars are the 
statistical fit uncertainties, and do not include systematics. The vertical error bars 
in the inset are dominated by the temperature uncertainty. Data without error bars 
have uncertainties smaller than the symbols. 

bottom 23 cm of both the inner and outer tubes were sanded with 120 grit paper 
to create an optically dull finish. The remaining 42 cm exhibited the shiny surface 
characteristic of the epoxy coating. The thermal conductance of the two sections was 
measured in the range 50 mK to 300 mK. The section with diffuse optical reflection 
had a conductance three times lower than the other section. We thus conclude that 
diffuse reflection reduces the conductivity by a factor of 3 x 42/23 ~ 5. 

3.4.3 Heat Capacity 

The heat capacity of the cell and dissociator have been measured together in a typ- 
ical trap loading sequence. The dissociator was held by a temperature controller at 
280 mK for 7 sec while pulsing an rf discharge inside. The heater and discharge where 
then turned off, and the cell plus dissociator assembly was cooled. The temperature 
of the bottom of the cell is shown in figure P-15| . The cooling rate is expected to be 
dominated by the large heat capacities of the large copper pieces which comprise the 
dissociator and mixing chamber. Note that the molar specific heats of copper and 
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Figure 3-14: Thermal conductance across the Ag sinter, from the superfluid jacket 
to the metal that comprises the dissociator. The two data sets are from the two 
thermometers at the top of the cell. Simple acoustic mismatch theory predicts a 
dependence. This is factored out in the inset. The line, a fit to the data points, is S = 
(8.6 X 10~^ /xW/mK^'^) x T^'^. The error bars reflect the statistical flt uncertainties. 
Data without error bars have uncertainties smaller than the symbols. 

^He are nearly identical at 100 mK, and there is much more Cu than ^He. 



3.4.4 RF Heating 

When rf power is introduced into the cryostat the cell and mixing chamber are 
warmed. Ohmic heating occurs in the coax lines that deliver the power to the cell, 
and rf eddy currents cause heating of the metals in the vicinity of the coils. Mea- 
surements of these heating rates is complicated by imprecise knowledge of where the 
power is deposited and by disruption of the thermometer readout circuits by the rf 
flelds. Nevertheless, some information can be obtained. 

The most reliable piece of information relates to the ejection of atoms from the 
trap. If the cell becomes too warm, ^He atoms evaporate and careen through the 
trapped sample, knocking atoms out of the trap. In many cases the maximum rf 
power is set by this criterion. A test was done of this ejection process by simply 
ramping the rf frequency across the range of interest (23 MHz to 2 MHz) at various 
rf powers, and noting how many atoms were ejected. To prevent resonant ejection a 
large bias field was applied so that the minimum resonance frequency of the trapped 
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Figure 3-15: Temperature of the bottom of the cell while cooling from 280 mK. 
The cooling rate is dominated by the heat capacity of other metal parts, not the 
superfluid. 

atoms was always above the applied rf frequency. No ejection was noted for powers 
up to +27 dBm outside the cryostat (corresponding to about 30 mG fields). The 
cell thermometers indicated 220 mK (bottom) and 180 mK (top). The maximum rf 
powers used for the experiments in this thesis were +15 dBm, so ejection of atoms 
from the trap was not a problem. 

Another piece of reliable data is that the cell will stay below 160 mK when rf fields 
of roughly 5 mG are applied at frequencies of 23 MHz and below. This temperature 
is low enough to allow the cell to cool reasonably quickly to temperatures at which 
atoms may be dumped out of the trap and analyzed. 

We note that there are a few resonances at which significant heating occurs. These 
are above 25 MHz. One possible rf leak might occur where the superconducting 
twisted pair is joined to the coax. Further investigations of where heating occurs 
should be pursued. Methods could include a consideration of how fast objects cool 
when the rf is turned off and the time required for heat to reach various thermometers. 
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Chapter 4 



Manipulating Cold Hydrogen by 
RF Resonance 



Resonant RF ejection of atoms from the trap is a tool for manipulating and studying 
the hydrogen gas. This tool has not previously been applied to trapped hydrogen, and 
so the peculiarities of this system have not been explored. We find — as anticipated — 
that rf ejection can efficiently remove energetic atoms that would otherwise not find 
their way out of the trap over the magnetic saddlepoint barrier. Experiments are 



described in section TA. Furthermore, we find that rf ejection spectroscopy can 
constitute a valuable diagnostic technique near BEC, where the trap dump technique 
used previously fails. RF ejection spectroscopy is also an essential tool for measuring 
the magnetic field at the bottom of the trapping potential. These topics are covered 
in section W72. 



4.1 RF Evaporation 

4.1.1 Need for RF Evaporation: Orbits with Long Escape 
Times 



In section |2.1| we explained that, for low temperature samples in tightly compressed 
traps, the saddlepoint evaporation technique leads to one dimensional evaporation. 
Recall that atoms with high radial energy are unable to escape before having a colli- 
sion. In this section we demonstrate this problem experimentally. 

Three successive samples are identically prepared by cooling into a trap with 
depth et/ks = 1-1 mK, set by a magnetic saddlepoint. The trap radial gradient is 
(y/kB = 4:A mK/cm. After the evaporation cycle is completed, the sample is held 
for many seconds (~ 20), on the order of three characteristic collision times. The 
saddlepoint barrier is then reduced to zero in 5 s while the atom fiux is recorded 
as a function of barrier height by sensing the heat of molecular recombination on 
a sensitive bolometer (section [4.2.3|) . In its simplest interpretation, this fiux profile 
is the energy distribution of the trapped cloud. (Corrections to this picture will be 



described in section [4.2.1|) . Figure ^4-l| shows atom fiux profiles for the three samples. 
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Figure 4-1: Comparison of dump signals from samples with and without rf ejection. 
The sohd sloping line shows the trap depth as it is reduced from 1.1 mK to zero. 
The curves indicate the atom flux as the trap depth is reduced. The dashed curves 
correspond to a samples that were exposed to rf tuned to the trap depth for 9 s and 
23 s. The solid curve corresponds to a sample with no rf exposure. The rf fields eject 
a sizeable fraction of the sample and change the shape of the profile. The response 
time of the detection process is seen at the end of the ramp, at t=5 s. 



The largest profile is for a sample (sample 1) that was held for 10 s, and then dumped 
from the trap. For the smaller two profiles, a weak rf field resonant at the trap depth 
(et = hurf — 6) was applied for 9 s (sample 2) and 23 s (sample 3) at the point in 
the sequence where sample 1 was held. The samples were then dumped out of the 
trap in the same way as sample 1. We see that a significant fraction of the atoms 
have been ejected from samples 2 and 3. To be ejected, these atoms must have had a 
total energy that was above the magnetic saddlepoint energy. Since these atoms did 
not leave the trap in sample 1, they were apparently in orbits with escape times of at 
least ten seconds. This experiment demonstrates the long times required for energetic 
atoms to escape when the energy barrier is set only by a magnetic field saddlepoint 
at one end of the trap, and the reduction in escape time when the atoms can leave 
in any direction. Furthermore, the decrease in trap population between samples 2 
and 3 indicates the long thermalization time of the system as it equilibrates to the 
lower effective trap depth for the radial degrees of freedom. This trap had a radial 
confinement that was a factor of ~ 3 weaker than those used for attaining BEC as 
described in section E^, and so these experiments place a lower bound on the problem. 
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Figure 4-2: Demonstration of slow transfer of energy from axial to transverse mo- 
tional degrees of freedom. The five curves are the flux of atoms escaping the trap as 
the magnetic energy barrier at one end of the trap is slowly ramped to zero (linearly 
sloping line), for each of five samples. The samples labeled 1-4 correspond to vari- 
ous durations of exposure of the sample to 10 MHz rf fields with strengths of a few 
milligauss. Sample 5 was confined in a short trap, and was exposed to an rf field at 
10 MHz for 20 s. The field had ten times more power than for samples 1-4. 



4.1.2 Mixing of Energy 

It is postulated that there is only very weak coupling among the motional degrees 
of freedom of a particle moving in the trap potential. This weak coupling between 
transverse and axial degrees of freedom is the impetus for driving evaporation with rf 
fields. In this section we describe a simple experiment which demonstrates this very 
weak coupling. 

Four identically prepared samples are each cooled into a trap with depth = 
1.1 mK and radial gradient a = 4.5 mK/cm. The trap is significantly longer than 
the one shown in figure |3-4| , with the sample extending from roughly z = —6 cm 
to z = +9 cm. The sample is adiabatically expanded by opening the trap slowly, 
reducing the radial gradient to a = 1.3 mK/cm, but leaving the length constant. The 
peak sample density is approximately tiq = 4 x 10^^ cm"'^. The thermalization time 
is thus on the order of 10 s, and the dipolar decay time is on the order of 10^ s. 



Figure [4-2| shows the flux of atoms escaping the trap as the magnetic confinement 
field at one end of the trap is linearly ramped to zero. As explained above, this 
atom flux can be interpreted as an indication of the distribution of energies of the 
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trapped atoms. The four traces labeled 1-4 correspond to the four samples. The first 
sample was simply held in the trap for 20 s after preparation, and then dumped. The 
second, third, and fourth samples had an rf field (a few milligauss strength) applied 
for various times — 35 s, 100 s, and 200 s. The frequency (10 MHz) corresponds to 
slightly less than half the trap depth set by the magnetic field barrier at the end of 
the trap. Even though the rf field is produced by both the axial and transverse coils, 
the rf fields have significant strength (and appropriate polarity) only along the sides 
of the trap, and not at the ends. We expect the rf fields to only eject atoms with 
energy in the transverse (not axial) directions. 

Comparing sample 1 to sample 2, we see that about two-thirds of the atoms are 
ejected when the rf is applied for 35 s, and that the energy distribution is narrowed, 
a sign of cooling. As the rf exposure is continued (samples 3 and 4), a few more 
atoms are ejected, and the sample continues to slowly cool. Two time scales thus 
appear. The fast time scale (~ 10 s) must correspond to removal of atoms with high 
energy in the transverse degrees of freedom. These atoms interact directly with the 
rf fields. The slow time scale (~ 100 s) is for removal of atoms that originally had 
high energy in the axial degree of freedom, but low energy in the transverse direction. 
It apparently takes on the order of 100 s for enough energy to transfer from axial 
motion to transverse motion to allow the atoms to reach the resonance region. 

This interpretation is supported by the data shown in figure [4-2| for sample 5. 
This sample was prepared identically to samples 1-4, but the trap was short. There 
was significant field strength (and appropriate polarity) at the ends as well as the 
sides of the trap, so that atoms with high energy in either the axial or the transverse 
direction could pass through the resonance region (again, both transverse and axial 
coils are used, as is done for the routine rf evaporation in chapter P). The resonance 
region wraps completely around the sample, in contrast to samples 1-4 for which 
the resonance region is essentially a tube around half the length of the sample. We 
observe that the energy distribution is narrowed very significantly with only short 
application of the rf field. This happens because energy does not need to mix among 
degrees of freedom before ejection can occur. 



4.1.3 RF Field Strength Required for Evaporation 

As discussed in section |3.1.3| , the rf field strength used during evaporation must be 
strong enough to eject energetic atoms from the trap before they have a collision. We 
have experimentally determined this power requirement for a typical sample density 
and trap shape. A series of samples are identically prepared at T ~ 70 /iK, and then 
evaporatively cooled to a given trap depth using rf fields of varying power. After this 
rf cooling cycle the atoms are dumped from the trap, and the total number of trapped 
atoms is recorded. For very low rf powers essentially no rf evaporation occurs, and 
the sample is not depleted. For higher powers the sample is significantly depleted, 
up to an rf power that saturates the process. At this high power, energetic atoms are 
promptly removed before having a collision; no further sample depletion is possible 
with increasing power because the atom ejection rate is no longer bottlenecked by the 
rf ejection rate, but rather by the rate at which collisions in the gas create energetic 



66 



O 



'0,8 



o 
E 

2^ 0.6 

CO 

E 
o 



0.4 



CD 



12 



6 12 
rf field strength [mG] 



Figure 4-3: Determination of the field strength required for rf evaporation. Atoms are 
depleted from the sample by fields of various strengths during an rf evaporation cycle. 
The peak sample density is no ~ 7x 10^^ cm^^ and the trap radial gradient is a/ks = 
4.5 mK/cm. The horizontal scale is the peak rf field strength, as estimated from room 
temperature measurements. There is a factor of two scale uncertainty which arises 
from imprecise knowledge of the temperature dependence of the impedances and of 
the details of the field geometry. 



atoms. 

Figure [4-3| shows one such saturation measurement. The rf frequency is swept 
from 17 MHz to 9.7 MHz at various rf field strengths. Two thirds of the atoms are 
removed when the field strength is above 2 mG. This agrees well with the calculation 
in section p.l.3| . Note, however, the large scale uncertainty on the horizontal axis. 
Ideally we would plot an effective field strength which is the appropriate average of 
the ejection rate over the resonance surface (see section [4.2.2|) ; ignorance of the details 
of the rf field direction and strength makes this average difficult. We have therefore 
plotted on the horizontal axis the strength of the field at low frequencies, as measured 
at one point in the cell at room temperature for a given power applied outside the 
cryostat. A correction was made for the frequency response measured in section [4.2.5 



assuming a driving frequency of 10 MHz. The number of remaining atoms continues 
to slowly drop as a function of field strength even for fields much greater than the 
saturation field of 2 mG. This behavior is not understood. It might indicate the 
presence of atoms in high angular momentum orbits which only very rarely cross the 
resonance region. 
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4.2 RF Ejection Spectroscopy 



The simple technique of rf ejection spectroscopy is useful as a diagnostic of the tem- 
perature and density of the trapped gas near BEC. It has significant advantages over 
the trap dump technique used previously, which becomes completely ambiguous near 
the BEC threshold. For probing the condensate, the technique is inferior to optical 
spectroscopy (described in section p.l.lD , but it is much simpler to use. In this section 
we explain the problems which make the trap dump technique ambiguous. We then 
develop a theory that explains the rf ejection measurements. After an explanation of 
the bolometric detection process employed in the rf ejection spectroscopy, we show 
how the technique can be used to measure the magnetic field at the bottom of the 
trap and to measure the temperature of a cold sample. 



4.2.1 Ambiguities of the Trap Dump Technique 

In previous experiments the sample was studied by dumping it from the trap and 
measuring its energy distribution (see section 2.5 and 2.6 in Jl^); the distribution 



was inferred by measuring the flux of escaping atoms as the axial trap barrier was 



reduced (see figure |4-1| for example). Ambiguities arise, however, because of conflict- 
ing timescales. In order to measure an energy distribution, the flux detector must 
respond in a time Tdet which is short compared to the duration over which the trap 
energy threshold is reduced to zero, Tdump- Furthermore, the characteristic sample 
equilibration time, Tequiu must be long compared to Tdump so that the energy distri- 
bution does not change during the dump process. For the dense samples studied near 
the BEC transition the equilibration time is Tequu ~ 300 ms. The detector response 
time is T^et ~ 200 ms. A further crippling complication is that the characteristic 
escape time, Tesci for an atom with total energy greater than the trap depth must be 
short compared to Tdump so that all of the atoms in a given energy slice leave the trap 
before atoms from the next slice begin to come out; in the tightly compressed traps 
used near the BEC transition this escape time has been observed to be many tens of 



seconds, as demonstrated by figure |4-1| . The timescale hierarchy, 



"^esc '^dump '^equil (^■-^) 

is not achievable, and the trap dump technique fails to providing reliable measure- 
ments of the energy distribution of samples near the BEC regime. There is also a 
practical problem. As the axial confinement magnetic field is lowered, the trap shape 
changes. Often a magnetic field zero forms, further complicating the measurements. 

RF ejection spectroscopy, in which a resonance region is tuned through the trap 
to eject atoms with various potential energies, addresses all these problems. The trap 
shape remains unchanged, so no zeros are introduced. The particle escape time can be 
short because atoms are coupled out in all directions. Finally, the energy distribution 
of the sample does not change significantly during the measurement process because 
only a small fraction of the atoms are removed. 
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4.2.2 Theory of RF Ejection Spectroscopy of a Trapped Gas 



In this section we use the results of the Landau-Zener theory discussed in section 
3.1.21 to explain the rf ejection spectrum expected for a trapped gas, i.e. the flux of 



atoms ejected as a function of the frequency of the rf magnetic field. This expression 
is used to interpret the spectra presented later in the chapter. 

Consider a sample trapped in a magnetic potential V^(r) = fiB |B(rap(i')| described 
by the potential energy density of states function g{e). The bottom of the trap is 
at a magnetic field Bq = An rf magnetic field, oscillating at frequency z/^/ 

with strength B.^/, is superimposed on the trap, and resonantly interacts with atoms 
crossing a surface of constant potential Sres = hurf — For the analysis that follows 
we note that this surface area and its differential thickness is related to Q{e) de, which 
is the differential volume in the trap with potential energy in a range de around e. 
This volume may be written as an integral, over the area of the surface of constant 
e, of the differential "thickness" de/\a{T) \ of the surface: 

/de 
dA-—-- (4.2) 
|a;(r)| 

where a(r) is the gradient of the potential at position r. This result will be used to 
derive equation O . 



We wish to find the rate at which atoms leave the trap as a result of Zeeman sub- 
level transitions in the resonance region. In a patch of resonance surface of differential 
area dA at position r, the differential atom loss rate is dN = FP dA, where F is the 
flux passing through the surface from either direction and P is the probability of the 
atom emerging from the resonance surface in an untrapped Zeeman sub level. Equa- 



tions |3.23| and p.24| indicate that, in the limit of low rf field strength, the vast majority 
of transitions are made from the d to c Zeeman sublevels. The transition probability 
is P = Pc = 2(1 — p)p, where p = 1 — exp{—() and C(r,'y) = 7r/i^P^j(r)/a;(r)t>/i is 
the adiabaticity parameter for an atom with velocity v passing through resonance 
at a position r. Here Brf is the component of the rf field that is perpendicular to 
the dc trapping field. In the weak field approximation, P ~ 2p and p ~ The 
transition probability exhibits a velocity dependence. The flux of atoms with velocity 
V passing through the surface is F{v) = n{r)fy{v)v where fviv) is a normalized one 
dimensional velocity distribution function, such as derived in appendix |^, and n(r) 
is the density at position r (either the Bose-Einstein or Maxwell-Boltzmann density, 
as appropriate). Since both P and F exhibit a velocity dependence, the differential 
atom loss rate must be derived from an integral over velocities 

dN = dA J dv n{r)f{v)v2C{r,v) 

= dA f\ ! ^ \ (4.3) 

The total atom loss rate T{e), through a resonance surface at potential energy e. 
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J^{e) = J dN 



= ^-^n{e)B!^,{e),{e). (4.5) 

Here we have defined i?e// as an effective rf field strength obtained by a weighted 
average, over the resonance surface, of the magnitude of the component of the rf field 
that is orthogonal to the trapping field, ^trap- The weight function is the thickness 



of the surface. Using equation 4.2 



^ f-L / dA '°Y';'\^o?'TI''"1 ■ (4.6) 
" \e(e)M.)„ |a(r)|BLp(r) / 

If the polarization of B^/ is randomized all over the resonance surface, then the cross 
product in equation |4]^ is essentially a scalar product, and Beff{e) is the average, over 
the equipotential surface, of the square of the rf field strength. This approximation 
is similar to a random phase approximation. Furthermore, if the length scale for 
variations in the field strength is significantly different from the size of the trap, then 
the field averaged over the surface will be independent of the equipotential chosen. 
Certainly if the field is uniform over the trap (characteristic length long compared 
to size of trap), then the average of the field over the equipotential surface does not 
depend on which surface is chosen. In the other limit, if the characteristic length 
scale of the fluctuations in the field strength is short compared to the size of the trap, 
then the variations average out over the surface regardless of the surface chosen. The 
result in either case is that B^ff is independent of e. The field generated in the cell 
by a given rf driving power, measured outside the cryostat, can vary as a function 
of frequency; this "transfer function" is an artifact of the reactance of the coils and 



the drive line impedance, and is measured in section [4.2.5| . This dependence is not 
included here. 

In our experiment these two approximations (random orientation and incommen- 
surate length scale) are not so bad; the rf field is a superposition of the fields generated 
by both the transverse and axial coils. The strength of the axial field varies on a length 
scale much shorter than the trap length, but much longer than the trap diameter. 



due to the field compensation scheme employed. Figure |3-4| shows the strength of 
the axial field along the symmetry axis of the trap. The field of the transverse coil 
is essentially uniform all across the trap. Since the relative strengths of the rf field 
components vary all over the trap, the direction of the field also varies. Tests of this 
simplifying assumption are outlined in section |4.2.5| . The relative phase of the axial 
and transverse field is assumed to be zero, independent of frequency and position. 
Violations of this assumption have not been considered. 

If Beff has no strong dependence on e, then the flux of atoms ejected from the 
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trap for a given u^f is proportional to the number of atoms with potential energy e, 
n{e)g{e), as found in equation A distribution of potential energies may thus be 
obtained by sweeping the resonance surface through the trap. This distribution can 
be used to infer the sample temperature. 

In the derivation of equation [4.3| , we made the approximation that the rf power 
was low. We assumed that ( <ti 1. Inverting equation H, this requires that Brj ^ 

^Javh/rriJ^. For typical parameters mv'^/2 = (50 /iK)fc^/10 and a/ks = (16 mK/cm) 
the requirement is that Brf 20 mG. Near the bottom of the trap where the 
potential flattens the small gradient implies that the experimenter must use very low 
rf field strength to remain in the low rf power regime. 

A further approximation used in the interpretation of the spectra is that the 
sample's energy distribution does not change while the spectrum is being measured. 
This condition could be fulfilled by measuring the spectrum quickly compared to the 
rethermalization time, but, as explained above, the response time of the detection 
process is too slow. Another way to realize this condition is to remove only a small 
fraction of the atoms while obtaining the spectrum. The number of atoms ejected 
from the trap in a duration Tspec is N eject = lo^^^" dt T{e{t)). For a linear ramp of the 
resonance from the bottom to the top of the trap, we have e{t) = {et/Tspec)t and 

Neject = '^^^^^"ff'^'P'"' f n{e) g{e) de 
net Jo 

Tist 

where N is the trap population. A spectrum should be obtained over a duration, r^pec, 
much longer than the detector response time; in practice we use Tspec ~ 5 s. For a trap 

depth of 1 mK we obtain the condition i?e// <^ \l'hetllT^Tspecl y^B ~ 200 /iG. This 
low ejection rate condition can be difficult to fulfill while simultaneously enjoying a 
large signal-to-noise ratio in the detector. 



(4.7) 



4.2.3 Bolometric Detection 

Atoms coupled out of the trap by rf resonance are detected by sensing a small fraction 
of the molecular recombination heat they liberate. The detection process is explained 



in detail by Doyle fl^. Atoms (in the c-state) leaving the trap stick to the cell 
wall. A fraction experience a spin fiip to the h- or a-state, opening the way for rapid 
recombination. If the residence time on the wall is long enough, and if the time 
associated with sliding along the wall out of the detection region and into a high 
field region is long enough, and if the surface density is high enough, then the vast 
majority of the atoms will recombine. The heat deposited on the bolometer will be 
proportional to the number of atoms which escaped from the trap. The picture is 
modified slightly for the rf spectroscopic studies presented here because the atoms are 
ejected in the c-state instead of the doubly spin-polarized (i-state treated by Doyle. 
State mixing should occur rapidly, however, induced by magnetic impurities on the 
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Figure 4-4: Impulse response of detection system. An rf pulse is applied for 200 ms, 
and the resulting bolometer signal is recorded. For this pulse the fit delay time 
is Tdeiay = 277 ms, the fall time is Tfaii = 340 ms, and the effective rise time is 
Trise = 251 ms. The heat in this pulse is 2.6 pJ, which corresponds to the heat 
from 3 X 10^ recombination events. The quiescent power applied to the bolometer 
was 100 pW, the bolometer temperature was 180 mK, and the cell temperature was 
90 mK. 



walls, and so the remaining considerations are still valid. They indicate that the 
cell walls should be cold (80 mK) and that the atom flux should be high. Both of 
these conditions impose limitations on the techniques described here. We observe 
well resolved rf spectra with atom fluxes on the order of 10^ atoms per second falling 
out of the trap toward the bolometer. 

Crucial to the interpretation of the spectra presented is an understanding of the 
time response of the bolometric detection process, as suggested above. The impulse 
response function of the detection system was measured directly using a series of 
200 ms rf pulses at various frequencies for samples of two different temperatures in 
the range of interest. Figure shows a typical pulse and an analytic fit of the form 



go 



f^-t' /Tfall _ Q-t' /Trise\ t' > 

rise \ / 



q(t') = I r/all-rr^se \ J - (4.8) 

1 t' < ^ ^ 

where t' = t — Tdeiay This form was chosen to model the time constants for an atom 
to stick to the wall and to find a recombination partner. The fit parameters were 
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consistent to about 10% for all the pulses, and the average values were Tdeiay = 299 ms, 
Trise = 246 ms, and Tjaii = 337 ms. The dependence of these parameters on the film 
thickness has not been studied systematically, but at low film thicknesses the H-He 
binding energy increases and the atoms stick promptly; Tdeiay decreases. Because 
the rf pulse has finite duration, the measured response is not, strictly speaking, the 
impulse response function. Nevertheless, it is a good approximation, and the analytic 
fit parameters quoted above were used in calculating the spectra presented in the 
following sections. 

The response time of the bolometer readout electronics was measured by heating 
the bolometer with short pulses of light from a HeNe laser. The response time was 
approximately 10^'^ s, negligible compared to the time constants for recombination. 



4.2.4 Measuring the Trap Bias Field 

The bias field in the loffe-Pritchard trap, 6/fiB, sets the characteristic sample tem- 
perature at which the trap functional form changes from predominantly linear (T ^ 
9/kB) to predominantly harmonic (T -C O/ks) in the two radial dimensions. It is thus 
an important quantity for analysis of Bose condensation experiments. Furthermore, 
to know the trap depth = hurf — 6, one must know both the trap energy barrier 
(set by the experimenter) and the bias energy. 

The bias is directly measured by sweeping up the frequency of the rf field until 



atoms begin to be ejected from the trap. Figure [4-5| shows such a sweep. The atoms 
begin to be coupled out of the trap at about t = 5 s, corresponding to a frequency 
ug = 700 kHz; in this case we find d/ks = ugh/ks = 33.6 fiK. RF ejection is a 
powerful and sensitive technique to quickly and accurately measure the absolute trap 
bias. 

The bias can also be inferred by deforming the trap potential until a field zero 
is created; atoms begin to experience Majorana flips, leave the bottom of the trap, 
and get detected by the bolometer. The zero is made by changing the current in a 
solenoid that creates a fairly uniform field along the length of the trap. Integration of 
the Biot-Savart law gives the field strength per unit current. This method of finding 
the bias field is not as accurate as the rf ejection method since the trap shape is 
deformed for the measurement. Nevertheless, the two methods agree to within the 
experimental uncertainties. 



4.2.5 Measuring the Effective RF Field Strength 

In order to analyze rf ejection spectra, one must know the effective rf field strength, 
Bf.ff given by equation TX. It may exhibit complicated polarization dependences if 
the random orientation and incommensurate length scale approximations made in 
the discussion of equation ^.6| are invalid. Furthermore, there may be a frequency 
dependence in the field strength generated in the cell by a given rf power applied 
outside the cryostat, as discussed in section ^.4.1| . Clearly a measurement of the 
efficiency is in order. 
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Figure 4-5: Measurement of trap bias using rf ejection. The rf field is switclied on 
at if: = and swept from 200 kHz to 1.2 MHz. The atoms begin to leave the trap at 
about 700 kHz. The sloping baseline signal results from atoms evaporating out of the 
trap as the sample continues to cool after the end of the forced evaporative cooling 
cycle. 



RF spectra were measured to verify the ejection efficiency for frequencies between 
1 and 17 MHz. The ejection rate at a given frequency, given by equation includes 
the density, the potential energy density of states, and the effective rf field strength. 
In order to simplify the measurement of the efficiency, the density was kept constant 
over the energy range of interest by using a warm sample, T ~ 1 mK> hVmax/kB- 
The potential energy density of states functional form is known well. The quantity 
obtained, then, is the ejection efficiency as a function of frequency. The flux ejected 
as the frequency of the rf fleld was swept upward is shown in flgure [4-6| for the axial 
coil alone, the transverse coil alone, and both coils together. Another sweep, taken 
with no atoms in the trap but with both coils driven, indicates rf pickup and heating 
effects. 

The simplest model of the frequency dependence of the ejection efficiency ignores 
polarization effects and accounts only for impedance matching (the voltage divider 
effects of the reactance of the coils and the 50 feed lines); we take 
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Figure 4-6: Ejected atom flux as the rf frequency is swept tlirougli tlie trap. Tlie 
lowest sweep indicates tfie rf pickup and lieating loads; no atoms were in the trap 
and both coils were driven. For the middle two sweeps only a single coil was driven. 
For the top sweep both coils were driven. The point at which atoms start to leave 
the trap, t ~ 2 s, indicates the bias field. Atoms stop leaving the trap when the 
rf field is extinguished, t ~ 20 s. At the lower frequencies the rf field strength is 
~ 3 X 10~^ G. Apparently the transverse coil is responsible for most of the heating 
at high frequencies. 



Figure |4-7| indicates the expected spectra for three values of uq. Evidently ujq ~ 
27r X 8 MHz, which agrees reasonably with the measured response of the coils below 
20 MHz at room temperature (figure |3-11D . This value is used in the calculated 
spectra presented next. 



4.2.6 Determining the Sample Temperature by RF Ejection 
Spectroscopy 



As suggested in equation [4.5| , rf ejection spectroscopy can give the temperature of 
the sample since it gives the potential energy distribution of the trapped atoms. This 
spectrum is different from the total energy distribution obtained by dumping the 
atoms out of the trap, as done in previous hydrogen experiments [1^, or by 
saturating the rf transition for the entire sample [^1], The difference between the 



total and potential energy distributions is indicated in figure ^-8| . As outlined above, 
the rf ejection spectroscopy technique is free of several complications implicit in the 
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Figure 4-7: Determination of the frequency response of the coils by a comparison 
of calculated and measured spectra. The measured spectrum has the background, 
shown in figure [4-6| , subtracted. The only adjustable parameters in the calculated 
spectra are the overall normalization, the baseline, and Uq. 



trap dump technique. RF ejection spectroscopy has been employed to measure the 
temperature of the coldest samples of trapped H to date, a factor of 40 colder than 
the photon recoil cooling limit. 

Figure [4-9| shows the spectrum of a sample that was cooled to T = 30 fiK by 
rf evaporation. This sample should contain a small condensate, although it is not 
visible here. The three calculated spectra in figure [4-9| assume quantum degeneracy. 
Furthermore, they include the measured time response of the detection system, the 
measured trap bias, and the measured effective rf field strength. Possible distortions, 
not considered in the calculations, arise from sample depletion during the rf sweep 
and changing detection sensitivity due to nonlinearities in the detection process at 
low atom fluxes. This second effect should be small, as suggested by the quality 
of the calculated rf ejection spectrum, which ignored nonlinearity. The first effect 
can be shown to be small by taking a series of such spectra at various rf powers, 
thus depleting the sample by various (small) amounts. Of course at lower powers 
the signal/noise ratio becomes smaller; the power used for the spectrum in figure 0^ 
involves a compromise between signal size and distortion. The apparent temperature 
is shifted no more than 10%. 

Determination of sample temperature through rf ejection spectroscopy is limited 
by the signal rate, which is proportional to Neject/Tspec in equation The technique 
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Figure 4-8: Probing the distribution of total energy and potential energy of a trapped 
sample. The total (potential) energy distribution function gives the number of par- 
ticles with a given total (potential) energy. The shaded area is the region of atoms 
probed by the technique. The solid lines indicate the trapping potential, and the 
dashed line is the trap depth. 

could be extended to much lower temperatures if the detection efficiency is improved. 

A practical limitation involves the requirement that the cell walls be cold enough to 
support a long atom residence time, and thus high recombination efficiency. Because 
of rf heating effects, this limitation usually requires the rf power to be reduced for 
several seconds before the rf spectrum is obtained. The sample can warm during this 
time since the trap depth is not as well defined. 
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Figure 4-9: RF spectrum of trapped gas at 30 /xK. The heavy hne is the raw data. 
The thin sohd hne is the inferred background, whose functional form has been verified 
independently. The three dashed lines are calculations of the spectrum for different 
temperatures assuming quantum degeneracy. The calculations include the time re- 
sponse of the detection system, the measured trap depth, and the measured effective 
rf field strength. The only remaining adjustable parameter is the normalization. The 
rf field strength used for this spectrum is ~ 200 /iG. 
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Chapter 5 

Observation of BEC 



In this chapter we describe the first observations of Bose-Einstein condensation of 
hydrogen. 

Our primary diagnostic technique is laser spectroscopy, which we briefiy describe 
here; detailed descriptions are given by Sandberg [jl6l and Killian [jl3l (Killian's thesis 
should be considered a companion to this chapter). With this probe we are able to 
monitor the density and temperature of the gas as the temperature is reduced by rf 
evaporation. When the gas is cooled into the quantum degenerate regime signatures 
of the Bose-Einstein condensate appear in the spectrum; these can be analyzed to 
reveal the size and population of the condensate. As anticipated, the condensate is 
huge, containing a factor of 50 more atoms than condensates in other experiments, 
and having a length on the order of a centimeter. We are able to study the time 
evolution of the condensate, which is governed by feeding of the condensate from the 
thermal gas. We explore some of the implications of these results for the production 
of high-fiux coherent atomic beams. 



5.1 Laser Spectroscopy of Trapped Hydrogen 

Laser spectroscopy of trapped atomic hydrogen has been pursued as both a technique 
to study and manipulate a trapped gas, and as a probe of the collection of individual 



atoms for metrology and fundamental measurements p4| , p!6| , |TJ, We focus here 
on spectroscopy as a tool to measure the density and temperature of the trapped gas. 
It has been suggested that the coherence properties of Bose-condensed atoms could 



be exploited to perform much improved fundamental measurements [B^ 



The principal transition in hydrogen, 1S'-2P, has a large transition rate, but is 
difficult to use for sample diagnostics because it has a large linewidth (100 MHz), 
it suffers from complicated magnetic structure, and it is in the uv (121.6 nm). The 
transition has been used by the Amsterdam hydrogen trapping group as a diagnostic 
for warm samples, and also for limited laser cooling and laser induced evaporation 



87| , [88[] . The Doppler effect can be used to measure the sample temperature, 
but the short 2P lifetime limits the technique to temperatures above a few mK. 
In order to circumvent the short 2P lifetime, the Amsterdam group employed two- 
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photon spectroscopy to the 35" (160 ns) and 3D (16 ns) states [^. To enhance the 
two-photon transition rate, one photon was tuned near the Lyman-a resonance and 
the other near the Balmer-a. In principle the 15'-35' transition could be used for 
temperature measurement all the way into the sub-/iK regime. 



5.1.1 Spectroscopy of the 1S-2S Transition 

We employ the 1S-2S transition and drive it with two photons at 243 nm, twice the 
Lyman-a wavelength. The 25" state is metastable with lifetime tq = 120 ms, and 
the natural linewidth is below 1 Hz. The transition thus constitutes a very sensitive 
probe, through shifts and broadenings, of the local environment of the atom. As 
explained below, the transition can be used to measure the density of the gas, in 
addition to the temperature. 

To excite the transition the laser beam is retroreflected to create a standing wave 
in the atom cloud. In this configuration the atom may absorb two photons from the 
same beam or one photon from each beam. Absorption of co-propagating photons 
involves the Doppler shift, and the excitation is said to be "Doppler-sensitive" (DS). 
In the absorption of counter-propagating photons the Doppler effect cancels to first 
order, and the process is termed "Doppler- free" (DF). 

The energy equation for two-photon excitation of an isolated atom from state i 
to /, with initial momentum pi and final momentum P/ = Pi + /i(ki + k2), where ki 
and k2 are the wave vectors of the laser beams, is 




2hu = yp'jc^ + {mc^ + 2hVoY — \IVi(^ + {j^c^Y- (5-1) 

where the rest mass of the atom in the initial state is m and 2vo = 2.466 x 10^^ Hz 
pO| is the unperturbed transition frequency. Expanding, we obtain 

. = z.. + + • (1 - + '^^^ + il-e)-^ +0{e% (5.2) 

V ' " v ' 

Audi Az/r 

Here Aum and Ai'd2 are the first and second order Doppler shifts, respectively, Aur 
is the recoil shift, and e = 2hvo/mc^ = 1.1 x 10~® is a relativistic correction which 
accounts for the mass change of the atom upon absorbing energy 2hvo- For hydrogen 
in the submillikelvin regime, Avm ^ 1 Hz and can be neglected. In the Doppler- 
sensitive configuration, ki = ^^2 and Av^f = 6.7 MHz. (All frequencies are referenced 
to the 243 nm laser source). At a temperature of 50 /iK, Azz/ji ~ 2.6 MHz, and thus 
the Doppler-sensitive peak is well separated from the Doppler-free. In the Doppler- 
free configuration, ki = — k2, and there is no recoil or first order Doppler broadening. 
A spectrum is shown in figure |5-1| . 

A complete discussion of the shape of the Doppler-free line in the various regimes 



of density and temperature is presented in the thesis of Thomas Killian We 
summarize his discussion by noting that the narrow intrinsic linewidth allows the 
spectrum to be used for measuring the density and temperature of the sample. The 
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density is measured through the cold colhsion frequency shift, the low temperature 
analog to the more familiar pressure shift. In the presence of a cloud of atoms of 
density n in the 15* electronic state, an atom in state a experiences a shift of its 
energy by an amount Snh'^aa-isn/m, where m is the atomic mass and tta-is is the 
s-wave scattering length which parameterizes collisions between atoms in state a and 
atoms in the IS* state. Although both the IS* and 25* energy levels are perturbed 
by surrounding atoms, |a25-_is'| ^ |cti5-is'|, so most of the cold-collision frequency 
shift is due to perturbations of the 25* level. The 1S-2S resonance is shifted by an 
energy AE = hxn, where hx = 8iTh'^{a2s-is ~ o-is-is)/''^ is the energy shift per unit 
density. The energy shift appears in the spectrum as a frequency shift (at 243 nm) 
= X^/'^- A detailed discussion of this mean-field picture of the energy shift has 
been presented by Killian [0]; here we merely emphasize that the process leads to 



a frequency shift proportional to the sample density, and thus the spectroscopy is a 
valuable tool for measuring the sample density in situ. We have measured this shift 
to be 

Xm = -3.8 ± 0.8 X 10"^° Hz cm^ (5.3) 



which is in reasonable agreement with calculations [^. The shift was measured for 
densities in the range 2 — 7 x 10^'^ cm~'^ and for temperatures between 100 and 500 
fiK. (For that study the sample densities were measured using the known dipolar 
decay rate constant, g, from equation |A.19| .) The frequency resolution of our laser 
system (currently 1 kHz) limits our use of this technique to densities above a few 
times 10^^ cm~^. The inhomogeneous sample density leads to a broadened lineshape 
that must be interpreted; the interpretation process is responsible for the majority of 
the ~20% uncertainty quoted for that study. 

At low densities the shape of the Doppler-free spectrum is governed by the finite 
interaction time of the laser beam with the atom as the atom traverses the beam 
(the cloud of atoms is much larger than the beam diameter). The characteristic 



linewidth (at 243 nm) is 6o = u/2TTdo where u = ^2kBT/'m is the most probable 
speed and is the waist diameter of the beam. The expected lineshape for a cloud 
of atoms having a Maxwell-Boltzmann energy distribution, interacting with a laser 
beam with a Gaussian beam profile, is the cusp-shaped S{5) oc exp(— |(5|/(5o) for uv 
laser frequency detuning 5 {5 in Hz). The linewidth is proportional to \/T, making 
this feature a useful probe of the sample temperature. For the beam waist diameter 
and laser frequency resolution achieved in these experiments this technique is reliable 
above 100 /iK. This temperature probe suffers from systematic effects arising from 
misalignment between the counterpropagating beams and divergence of the beam 
over the length of the trap, effects which make the lineshape difficult to interpret. In 
an improved apparatus these effects could be better characterized and controlled. 



5.1.2 Experimental Technique 

The 15'-25' transition is excited by narrowband continuous-wave radiation generated 
by frequency doubling the output of a 486 nm dye laser. This laser, which is frequency 
stabilized to an isolated Fabry-Perot reference cavity, routinely exhibits a linewidth 
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of 1 kHz and output power of 480 mW. After frequency doubling and spatial filtering 
there is typically about 9 mW of uv power in a Gaussian mode available at the window 
into the cryostat. 

The trapped gas is illuminated by the uv beam as it traverses the symmetry axis 
of the trap and is retroreflected by a spherical mirror. The beam focus (50 fim waist 
radius) is positioned at the minimum of the trapping potential so that maximum 
optical intensity overlaps maximum sample density. The atoms are excited during an 
illumination interval Texc ~ 500 /xs. Atoms in the metastable 25* state are held for a 
short delay r^eiay ~ 70 /is while background fluorescence dies away. Then a 10 V/cm 
electric field is applied for a short time Tquench ~ 10 /xs; the 2S state mixes with the 2P, 
which quickly decays emitting a Lyman-a photon which is counted on a microchannel 
plate detector. Only photons which hit the detector during the quenching interval are 
considered signal, and so the background count rate is much suppressed. The total 
detection efficiency should be about 2 x 10~^, limited by detection solid angle (0.13%), 
transmission through a window (60%) and a Lyman-a filter (10%), and detector 
quantum efficiency (estimated at 25%). Actual detection efficiency appears to be 
about 10 times lower than this estimate, possibly due to contaminants on windows 
|T^, The background count rate is between 200 and 700 counts/sec, giving about 8 
background counts per second of laser exposure assuming the excitation and detection 
timing sequence outlined above. The minimum detectable signal thus corresponds to 
4 X 10® atoms excited per second of laser excitation. Straightforward improvements 
to the optical access should easily yield a factor of ten improvement in detection 
efficiency. 

5.1.3 The 15-25' Spectrum of a Non-Degenerate Gas 

The spectrum of the trapped hydrogen gas slightly above the quantum degenerate 
regime is shown in figure |5-1| . As described in section |5.1.1| , there are two components 
of the spectrum, corresponding to absorption of co- or counter-propagating photons. 
The wide, low feature on the right (higher photon energy) is the Doppler-sensitive 
peak, and is offset by the recoil shift from the intense, narrow Doppler-free 
peak on the left. The solid line is the Gaussian lineshape expected for a Maxwell- 
Boltzmann distribution of kinetic energies in a sample at 42 /iK. The DS spectrum 
maps the velocity distribution through the Doppler shift; only atoms with the proper 
velocity v = XAu are Doppler shifted into resonance for a given laser detuning Au = 
u — {1^0 + Aur) iy is the optical frequency of the laser). Thus, only a tiny fraction of 
the sample is resonant, and the excitation rate is correspondingly tiny. The fraction 
of atoms which are resonant is proportional to the ratio of the excitation linewidth, 5q 
from time-of-flight broadening, to the Doppler linewidth = ^2kBT log(2)/mA^; 
in our apparatus this ratio is ~ 5 x 10^^. The very high densities of our trapped 
samples have allowed us to measure the DS 15'-25' spectrum of hydrogen for the first 
time. 

The Doppler-free spectrum gives information about the density distribution and 
temperature of the gas, as described above. Both of these effects are shown by the 
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Figure 5-1: Spectrum of non-degenerate gas. The intense narrow peak on the left 
results from Doppler-free excitation and the low broad peak on the right results from 
Doppler-sensitive excitation which probes the velocity distribution. The solid line 
on the right is the expected Doppler-sensitive spectrum for a nondegenerate sample 
at 42 /iK. At this low temperature the recoil-shifted Doppler-sensitive line is clearly 
separated from the Doppler-free line. The error bars are statistical. The ratio of 
peak heights should be taken as only a rough indication of the transition rate ratios 
because of variations in laser power and beam alignment. 



spectra in figure p-2| , just three of many successive spectra that were measured as the 
density of the sample slowly decayed. The initial spectra exhibit a large shift, arising 
from the cold-collision frequency shift at high density. The inhomogeneous density in 
the trap causes a broadening of the spectrum as well. As the sample density decreases, 
the shift and broadening also decrease. At low density, where the density-dependent 
shift is negligible, the exponential lineshape characteristic of transit time broadening 
is evident. The frequency shift of the spectra of high density samples is used to obtain 
the sample density, and the width of the transit-time broadened line is related to the 
sample temperature. 



5.2 Cooling into the Degenerate Regime 

Armed with an efficient cooling method (rf evaporation) and a sensitive probe of 
sample density and temperature (spectroscopy), we are equipped to push toward the 
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Figure 5-2: Doppler-free spectra at various sample densities. The trap depth is 
570 fiK, corresponding to a sample temperature T ~ 70 /xK. The spectra of the 
higher density samples are red-shifted and broadened by the inhomogeneous density. 
The lines connecting the points are a guide to the eye. The samples were confined by 
trap shape A, described in table |5.1|. The highest density is about 10^^ cm~3. 



quantum degenerate regime. The trap is loaded as described in section |L3|, and 



the sample is cooled by magnetic field saddlepoint evaporation. The currents in 
the magnetic coils that create the trapping potential are ramped to slowly change 
the trap from the deep, long, open shape used for loading the atoms to a short, 
tightly compressed, shallow trap from which rf evaporation is commenced. Figure |5^ 
indicates the process. In this section we discuss the path through density-temperature 
space taken by the sample as it is cooled into the quantum degenerate regime. 



5.2.1 Cooling Paths 

Previous attempts to achieve BEG in hydrogen were thwarted by inefficient evapora- 
tion at the lowest temperatures. In order to verify that we have solved this problem, 
the progress of the sample through phase space (density-temperature space) is mon- 
itored as the sample is cooled. We find that the sample is cooled efficiently. We also 
find and explore the phase boundary where the sample undergoes phase separation 
into a condensate and normal gas. This first evidence for phase transition behavior 
is described below. 
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Figure 5-3: Schematic diagram of trap shape as sample is loaded and evaporatively 
cooled to BEC. Room temperature molecular hydrogen (upper left diagram) is cooled 
to 275 mK by cryogenic techniques. It is dissociated, and the atoms are loaded into an 
open, long trap that is about 500 mK deep (upper right diagram). The magnetic field 
in the p — z plane is shown. Note the very different scales for the radial coordinate p 
and the axial coordinate z. During trap loading, the atoms settle into the region of 
low field at the center of the trap. The sample is then cooled by magnetic saddlepoint 
evaporation. It is compressed radially and axially into a trap that is nearly harmonic 
along the z axis and predominantly linear along the radial axis (lower right diagram) . 
It is this trap shape that is labeled "A" or "B" (see section |5.2.1| ). RF evaporation 
is used to further cool the sample while the trap shape remains fixed (lower left 
diagram). The entire cooling process takes about 5 minutes. 



A series of measurements were made to map the sample's trajectory through phase 
space as it was cooled into the quantum regime. Identical samples were cooled into a 
trap of depth et/ks = 1-1 mK. This process used magnetic saddlepoint evaporation 
exclusively. The trap shape was then held constant while the trap depth was reduced 
by rf evaporation. (Three trap shapes will be discussed in this chapter, A, A', and 



85 



7^ Q T* Q TTl OTCtT* 

JJcLi cLIlicLt:! 


trap A 


trap A' 


trap B 


a/ks (mK/cm) 


15.9 


9.5 


P/kB (/iK/cm^) 


25 


e/kB (/iK) 


27 ±2 


35 ±2 


34 ± 2 


no(T = 120 /iK) (xlO^^ ^ni-S) 


8 


1.7 



Table 5.1: Summary of parameters describing the three trap shapes used in this 
chapter. The parameters a, and 9 describe the loffe-Pritchard potential; a and (3 
are calculated, and 6 is measured. The peak density at T = 120 /iK is noted. 



B. Their characteristics are summarized in table 5.1. Here we refer to trap shape 



A.) The evaporation path (the trap depth as a function of time) was identical for all 
samples, but the different samples were taken various distances along the path. The 
rf frequency f{t) that sets the trap depth = hf{t) — 6 was ramped according to the 
function^ f{t) = f start {fend/ fstartY^^'^^^ where the start frequency is f start = 23 MHz, 
the end frequency is fend = 2 MHz, the maximum ramp duration is c? = 25 s, and 
the shape exponent ^ = 1.5. The ramp was executed for a duration r^. For t > Tr, 
the frequency was maintained at /(r^). Evaporation could still occur, but forced 
evaporation halted. After rf evaporation was completed, the sample density was 
measured using the cold-collision frequency shift of the DF spectrum as described 
above. The densities thus measured are shown in figure |5-4| . The density initially 
increases as the trap depth (and thus temperature) is reduced because the atoms are 
confined to a smaller volume closer to the bottom of the trap. Some atoms are lost 
through evaporation, but the good efficiency of evaporative cooling leads to overall 
phase space compression. The cooling efficiency parameter 7 = dlogD/dlogN (the 
phase space density is D = uqA^) is 7 ~ —1.5 during this stage of the coohng. 

At a trap depth of about et/kB = 280 /iK the phase space trajectory exhibits 
a dramatic kink. Further cooling leads to decreasing density in the normal gas; 
apparently population is being rapidly removed. The location of the kink and the 
slope of the subsequent trajectory agree well with the BEC phase transition line, 
ric = 2. 612 {27imk bTY ^ / assuming a ratio of trap depth to temperature of = 
et/kBT ~ 6, which is reasonable. The "excluded region" above the line has the 
following interpretation: as the sample is cooled, the thermal cloud cannot support a 
peak density above ra^; extra atoms are transferred to the condensate. The condensate 
is not directly observed in these measurements because its density is so large (as 



discussed in section |5.3.1| ) that its spectrum is shifted far away from the spectrum of 
the thermal cloud. 

The magnetic trap shape used for these experiments (trap shape A) was chosen as 
a compromise between high radial compression and adequate magnetic saddlepoint 
evaporation efficiency. High radial compression is desirable because the correspond- 
ingly higher densities allow BEC to be reached at a higher temperature. As discussed 



^The evaporation path was not optimized before most of the data presented in this chapter was 
obtained. For future studies an optimization of the functional form of the path is desirable. 
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Figure 5-4: Path of non-condensed portion of the sample through n — T space as 
the sample is cooled into the degenerate regime. Paths are shown for samples in two 
different traps: open circles — trap A, closed circles — trap B. The trap depth is set 
by the frequency of the rf magnetic field, and above degeneracy the temperature is 
about l/?7 ~ l/8th the trap depth. The density is measured through the cold-collision 
frequency shift. The non-condensed portion of the sample cannot have density greater 
than the BEC critical density ric oc T^/^ and thus the upper left region of the phase 
diagram is not accessible. The condensates created in the two traps contain about 
the same number of atoms. 



in section ^.2.2| , evaporation and the approach to BEC should be more efficient at 
higher temperatures because the ratio of "good" to "bad" collisions is more favor- 
able. On the other hand, high radial compression is undesirable because the magnetic 
field saddlepoint evaporation process (see figure p3) that cools and compresses the 



sample into the tight trap is less efficient. As discussed in section |2.1| , tight radial 
compression leads to more regular motion of the trapped atoms, and thus to reduced 
evaporation dimensionality. 

A more open trap, trap B, has also been used. A typical phase space trajectory is 
shown by the filled circles in figure 5^ . The evaporation was done more slowly, with 
maximum duration d = 55 s, but was otherwise identical to the path described above. 
The evaporation efficiency parameter is 7 ~ —2.7 for cooling above the degenerate 
regime. The higher evaporation efficiency has not been systematically studied, but 
could arise simply because the trajectory duration d was optimized for trap B, but 
not trap A. 
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Computer models of the cooling process have been constructed, and are described 
in appendix 0. The trajectories through phase space calculated by them are in good 
agreement with observations. 

5.2.2 Evaporation Strategies 

There are many paths that can be taken through phase space to the BEC phase 
boundary. In this section we discuss the experimental limitations of the present 
apparatus and their implications for the choice of evaporation path. 

As shown in section |A.2.2| , evaporation is more efficient in hydrogen at higher tem- 
peratures. There the ratio of "good" collisions to "bad" collisions is higher. "Good" 
collisions are those which lead to cooling, and they occur at a rate Tevap oc Tcoi oc 
nav oc n\/T for a sample of density n. "Bad" collisions are those which lead to 
dipolar relaxation and heating. They occur at a rate TiUp oc gn. The collision ratio 
^evap/^dip OC Vt IS higher at higher temperatures. More phase space compression 
can be accomplished per atom lost from the trap. This higher efficiency makes it 
desirable to do most of the phase space compression at higher temperatures. 

The competing consideration is that full evaporation dimensionality, and thus effi- 
ciency, can only be realized using the rf evaporation technique (instead of saddlepoint 
evaporation), and this technique is only usable for trap depths less that about 1 mK. 
Creation of deeper traps requires higher frequency and higher power rf fields. These 
heat the current apparatus too much. 

In an ideal apparatus the rf evaporation would be initiated at much higher tem- 
peratures, phase space compression would lead to higher sample densities, and the 
bulk of the cooling and compression would be accomplished more quickly and more 
efficiently than in the present experiments. The sample could then be cooled at nearly 
constant phase space density to cross the BEC phase transition line at the desired 
density (one may desire smaller densities with correspondingly longer lifetimes). In 
the experiments described here more phase space compression was accomplished at 
lower temperatures than is optimal because of the heating effects of the rf power used 
for evaporation. Improvements to the rf parts of the apparatus should allow larger 
condensates to be created. 

To test these ideas about efficiency we cooled to BEC in two traps, A and B 
described above. (These labels refer to the field configuration at low temperature. 
Identical trap shapes are used for initial trap loading (T = 40 mK, uq = 2x 10^^ cm~^) 
and the first stage of saddlepoint evaporative cooling. The trap shape is converted 
into A or B during the last stage of the saddlepoint evaporation.) The evaporation 
occurred along different paths for the two traps, but at the end the phase space density 
was nearly identical; the condensate population was about 10^ atoms for both traps. 
One expects path B to have lower overall evaporation efficiency (and thus fewer 
condensate atoms) because more of the phase space compression is accomplished 
at lower temperatures (where the evaporation efficiency is lower because the ratio 
of good to bad collisions is less favorable). On the other hand, if rf evaporation 
is more efficient than saddlepoint evaporation, then one expects path B to result 
in better evaporation efficiency, and thus more condensate atoms, since more of the 
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phase space compression is accomplished with the more efficient rf evaporation. Since 
the two paths result in the same final phase space density, we conclude that the two 
effects cancel. This necessitates that rf evaporation be more efficient than saddlepoint 
evaporation. It would be advantageous to use rf evaporation for more of the cooling. 



5.3 Spectroscopic Study of the Degenerate Gas 

When the sample is cooled into the quantum degenerate regime we observe signa- 
tures of the condensate in the 1S-2S spectrum of the gas. In this section we discuss 
two of these signatures from which we extract the chemical potential and tempera- 
ture. Knowledge of these thermodynamic parameters allows us to estimate several 
properties of the system. 

Bose condensation involves the macroscopic occupation of the lowest energy quan- 
tum state of the system. In a trap this state is concentrated at the minimum of the 
trapping potential, and it has very small kinetic energy, given approximately by the 
Heisenberg momentum-position uncertainty principle. In the Doppler-sensitive spec- 
trum, which maps the momentum distribution, one would expect an intense line at 
zero detuning from the recoil shift Az/^j, rising above the background spectrum WE 



Figure |5-5| shows the DS spectrum of the normal gas and the condensate together. 

Instead of appearing as a very narrow, intense peak centered at zero detuning, 
however, the line is red-shifted and broadened by the very high and inhomogeneous 
density of the condensate. This same effect is observed in the DF line, displayed in 
figure |5-6| . The condensate spectrum has a very different shape and shift than the 
normal gas spectrum, implying very high densities and a density distribution that is 
different from the normal gas. 

The spectroscopic technique of probing the condensate differs significantly from 
the spatial imaging techniques used by other groups, providing somewhat comple- 
mentary information. The probe is not limited to length scales larger than the spa- 
tial resolution of the imaging optics, and so detection of very small features, such 
as condensate "droplets" would be possible. Also, we are able to measure the 
momentum distribution of the normal gas in situ. Finally, the coherence length of 
the condensate can, in principle, be obtained directly from the broadening of the 



spectrum, as recently demonstrated in Na by Stenger et al. |9^. Because of the 
large condensate densities, we are only able to put a very small lower bound on the 
coherence length for our H condensates. 



5.3.1 Peak Condensate Density 

The DF spectrum of the condensate, shown in figure 5^, can be understood using a 
local density approximation; atoms excited from different regions of the condensate 
exhibit correspondingly different frequency shifts. The signal size S" at a given de- 
tuning A is proportional to the number of condensate atoms on a surface of constant 
density ricond = 2A/xc, where Xc is the (two photon sum) frequency shift per unit 
density for excitation out of the condensate. (The relation of Xc to Xm, the measured 
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Figure 5-5: Doppler-sensitive spectrum of degenerate gas. Zero detuning is taken 
at center of the recoil-shifted spectrum, which is detuned 6.70 MHz blue of the 
Doppler-free resonance, open circles — normal gas, filled circles — spectrum focusing 
on condensate. The dashed line is a fit to the normal fraction data which assumes 
a Maxwell-Boltzmann velocity distribution. The solid line is a calculation which as- 
sumes a Bose-Einstein distribution. The dot-dashed line is the condensate spectrum 
expected for a Thomas-Fermi wavefunction in a harmonic trap, when the dominant 
spectral broadening is the cold-collision frequency shift. The temperature of the Rose 
calculation was chosen to fit the observed spectrum. The condensate spectrum is 
a composite of data from the first sweep of each of three loadings of the trap; the 
spectral sweep took 0.94 s. The normal spectrum is a composite of data from the 
first three sweeps of each of two loadings of the trap; the spectral sweep took 0.94 s. 
This data is for trap shape B. 



frequency shift per unit density for excitation out of the thermal gas quoted in 



equation |5.3| , will be addressed in section |5^. Although there is significant ambiguity 



it is clear that either Xc = Xm oi Xc = Xm/2.) The spectrum calculated using this 
local density approximation in a harmonic potential is parameterized only by the 
peak shift, Ap, and the overall scaling factor Sq. The normalized spectrum is 



5o-7^¥#2Ta/A- ; A„<A<0 



5(A) = <( -Ai-^)^^ ' -f---" _ (5_4) 

otherwise 
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Figure 5-6: Doppler-free spectrum of degenerate gas. The intense, narrow feature on 
the right is due to the normal gas, and is shrunk vertically by a factor of 40 to fit in 
the plot. The wide, low, red-shifted feature is the condensate. The lineshape is that 
expected for a Thomas-Fermi wavefunction in a harmonic trap. The spectrum is a 
composite of the first two spectral sweeps from each of 10 loadings of the trap; the 
sweep took 0.67 s. The normal gas spectrum is a single sweep from a single loading 
of the trap. This data is for trap shape A'. 



(The negative signs are appropriate for red density shifts: Ap < 0.) A fit of this 
function to the measured spectrum is the dashed line in figure |5-6| . The peak shift 
is 920 ± 70 kHz, which indicates a peak condensate density Hp = 2Ap/xc = 4.8 ± 
0.4 X 10^^ cm~^ (assuming Xc = Xm] the uncertainty reflects the statistical fitting 
uncertainty, and ignores the 20% uncertainty in Xm)- If we assume Xc = Xm/2 the 
density is twice as big. The quality of the fit indicates the utility of the Thomas- Fermi 
wavefunction approximation. 

A quantum mechanical approach may also be used to understand this spectrum, 
as developed in the thesis of Killian [|13]. Before excitation, an atom is in the IS* 
electronic state, and its motional quantum state is the condensate wavefunction. 
After excitation the atom is in the 25* electronic state; its motional state has been 
projected onto the motional eigenstates of the modified potential it experiences. This 
new potential, a sum of the trap potential and the mean-field interaction energy 
for a 25" atom in a IS" cloud, is very different because the new mean-field energy 
has a different sign and is much stronger. The parabolic condensate density profile. 
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with density much larger than the background thermal cloud density, gives rise to 
an attractive parabolic potential for the 25" atoms. The motional quantum states 
the 25* atom finds available are simple harmonic oscillator states. The excitation 
spectrum is the probability of exciting atoms from the condensate wavefunction to 
each of these 25* motional states. This probability is the wavefunction overlap, which 
is straightforward to calculate. Calculations by Killian show very good agreement 
between the quantum and local density approximation treatments of the spectrum. 

The peak condensate density indicates the number of atoms in the condensate, 
given a knowledge of the trap shape. (This data is for trap shape A'). Using equation 



2.14 we obtain Nc = 1.2 ± 0.2 x 10 . This is larger by a factor of 50 than the biggest 



condensate created so far in other labs Using equation p.l5| , the characteristic 



dipolar decay rate is T2^c = A^2,c/A"c = 1-7 s~^. The chemical potential for this 
peak density is /i = UpUo = 1.9 fiK. The extent of the Thomas-Fermi wavefunction 
is determined from equation |2.13| by setting the density to zero at the edge of the 
condensate: ricondUo = — Vip{p, z) = 0. We obtain 



and 



Pmax = -aZ/x^ + 2p9 = 7.3 fim (5.5) 
a ^ 



'^ = 2.8 mm. (5.6) 



For trap B, similar data has been obtained. The peak shift of Ap = 620 ± 20 kHz 
corresponds to a peak density 3.26 ± 0.10 x 10^^ cm~'^ (assuming Xc = Xm] the 
uncertainty refiects the statistical fitting uncertainty, and ignores the 20% uncertainty 
in Xm)- In this more open trap, the condensate population is Nc = 1.19 ±0.10 x 10^, 
the chemical potential is n/ks = 1.28 ± 0.04 /iK, the diameter is 2pmax = 20 fim, 
and the length is 2zmax = 4.5 mm. If we assume Xc = Xm/2, then the peak density is 
twice as high, and the derived condensate parameters are correspondingly different. 
These results are summarized in table 



5.3.2 Sample Temperature 

The temperature and chemical potential of the trapped gas allows a full thermo- 
dynamical description of the system, assuming a near equilibrium situation. The 
chemical potential was obtained above. The temperature is extracted here using two 
methods, both of which suffer from uncomfortably large uncertainties. Nevertheless, 
in this section we characterize the temperature with uncertainty on the 20% level. 

The primary thermometry tool is the DS spectrum of the normal gas. As described 
in section f>.l.l\ above, the spectrum is a map of the velocity distribution of the sample. 
Far from quantum degeneracy this (classical) velocity distribution is Gaussian, and 
the normalized spectrum is 
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parameter 


trap A' 


trap B 


a/ks (mK/cm) 


15.9 


9.5 


P/kB (/iK/cm2) 


25 


25 


e/kB (/xK) 


35 ±2 


34 ±2 




28 





220 


et/A;ij (/iK) 


272 


214 


Te (/iK) 


~ 65 




50 


(kHz) 


920 ± 70 


620 ± 20 




Xc = Xm 


Xc = Xm/2 


Xc = Xm 


Xc Xm/2 


Up (xlO-*^^ cm~^) 


4.8 ± 0.4 ± 1 


9.7 ±0.7 ±2 


3.3 ±0.1 ±0.7 


6.5 ± 0.2 ± 1.3 


^i/kB (/iK) 


1.9 


3.8 


1.3 


2.6 


iVc (xlO^) 


1.2±0.2ni) 


6.6 ±1.3(1^) 


i.2±o.i(l[5:5^) 


6.7 ±0.5(1^) 




1.7 


3.3 


1.1 


2.2 


2pmax (/im) 


15 


21 


20 


28 


2zmaa; (mm) 


5.5 


7.8 


4.5 


6.4 



Table 5.2: Summary of parameters describing the two trap shapes used for achieving 
BEC, and summary of the properties of the condensates. The parameters a, (3, and 6 
describe the loffe-Pritchard potential; a and /3 are calculated, and 9 is measured. The 
trap depth at which BEC is first observed is et^c, and the trap depth used for the obser- 
vations listed here is e^. The approximate temperature in the presence of the biggest 
condensates is Tc. The peak cold-collision frequency shift in the condensate is Ap. 
The remainder of the table is divided to show the implications of assuming Xc = Xm 
or Xc = Xm/2 (recall the sum frequency shift is Xm = — 3.8±0.8 x 10^^° Hz cm^). The 
peak condensate density is Up, and the chemical potential is /i. The number of conden- 
sate atoms is Nc. The characteristic condensate dipolar decay rate is T2,c = N2^c/Nc. 
Finally, the length and diameter of the condensates are given. The uncertainties are 
divided into a component which depends on the present experiment (first number), 
and a component reflecting the 20% uncertainty on Xm (second number). 



where a = ^JkBT/mX^ = 373.6yT/l /iK kHz is the linewidth parameter and the 
detuning A is measured in Hz. We have not measured the DS spectrum of the gas 
far from quantum degeneracy because the low sample density reduces the signal rate 
below the background noiseQ. There is sufficient signal near the degenerate regime, 
however. For a degenerate Bose gas at a given temperature, the velocity distribution 
is significantly narrower than its classical counterpart because the particles are dis- 
tributed toward lower energies. In a homogeneous potential the DS spectrum of an 



■^The signal rate at zero detuning is S{0) oc no/\/T, whereas the phase space density is 13 cx 
uq/T'^/'^. The signal rate S oc DT falls less rapidly than the phase space density as the temperature 
increases, indicating the potential feasibility of obtaining DS spectra of samples far from quantum 
degeneracy. 
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ideal gas is cusp-shaped: 



S{A) oc log 

In a trap, however, only the atoms at the very bottom of the trap are in the quantum 
degenerate regime. Atoms in regions of potential energy e > ksT experience a reduced 
phase space density, parameterized by an effective chemical potential /i' = (/i — < 
— l/c^T; their velocity distribution is nearly classical. The measured spectrum is 
an average over the various velocity distributions in the gas. Appendix ^ presents 
a detailed treatment of Doppler-sensitive optical excitation in a trap, showing that 
if the laser beam is small, then atoms in regions of high phase space density are 
preferentially excited and the spectrum reflects the Bose velocity distribution more 
dramatically. The shot noise on our measurements is too big to allow a precise 
comparison between the shapes of the classical and quantum spectra, but evidence of 
the effect is noted in the spectrum below. 

Although it is interesting to probe the Bose-Einstein velocity distribution and 
compare it to the Maxwell-Boltzmann, the disadvantage presented by the effect is 
the dependence of the spectrum on the effective chemical potential, which leads to 
measurement uncertainties. The laser beam geometry and alignment must enter 
the interpretation of the spectrum, and these factors are difficult to measure directly. 
The alignment of the beam with the sample is verified by translating the trap radially 
relative to the laser beam while exciting the condensate. The condensate is used as a 
precise alignment marker in the center of the trap. With this technique we can center 
the trap radially over the laser beam to within about a half a beam radius (20 /xm), 
small compared to the thermal radius of the normal gas (100 /im). Systematic axial 
adjustment is not possible in the current apparatus. The uncertainty in the axial 
position of the beam focus is about 2 cm, about half the length of the normal gas 
cloud. 

To demonstrate the technique we study here the DS spectrum of a sample in 
trap B, in the quantum degenerate regime. Plotted in figure |5-7| is the DS spectrum, 
obtained immediately after forced rf evaporation is completed; the trap depth is 
et = 214 fiK. The first step in the analysis is to fit the wings of the spectrum with 
a Gaussian, an analytic form that is simple to fit. This fit is shown in the upper 
panel along with the spectra expected for temperatures higher and lower than the fit 
result by one standard deviation. Next, spectra are calculated assuming a truncated 
Bose-Einstein velocity distribution, as described in appendix |^. These calculations 
take into account the trap shape, trap depth, laser beam waist, and alignment of laser 
beam relative to the center of the trap (we choose ideal alignment here: beam focus 
at center of trap, beam aligned with the z axis). The calculations are done for a range 
of temperatures. The calculated spectra are fit in the wings by a Gaussian, and the 
spectrum with fit temperature closest to the fit temperature of the data is chosen. For 
the data in figure |5-7| , the fit temperature for the data is 36 ± 6 fiK. The calculated 
Bose spectrum with temperature T = 50 /iK also had a fit temperature of 36 fiK (see 
below). This spectrum, and those for T = 40 /iK and 55 /xK are plotted in the lower 
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MB velocity profile: 
- - - T=30 




laser detuning [MHz at 243 nm] 

Figure 5-7: Doppler-sensitive spectrum of degenerate gas. The condensate signal 
(point at -0.5 MHz detuning) is off scale. In the upper panel the wings of the line 
(all points except center three) are fit to a Gaussian, the expected lineshape for a 
Maxwell-Boltzmann (MB) velocity profile. The fit temperature is T = 36 ± 6 /zK. 
In the lower panel the same data is compared to calculations of the spectrum which 
assume a Bose-Einstein (BE) velocity distribution for a degenerate gas in trap B. 
Spectra calculated for three temperatures are displayed. We conclude that the sample 
temperature is T = 50 ± 10 fiK, corresponding to r] = 4.3 ±0.9. The BE velocity 
distribution fits the low velocity portion of the data better than the MB distribution, 
indicating the presence of quantum effects in the momentum distribution. The data 
is compiled from the first three sweeps of two samples; each sweep required 0.94 s. 



panel of figure 5^. We conclude that the sample temperature is T = 50 ± 10 yuK. 

When a Gaussian is fit to a Bose velocity distribution, calculated for temperature 
Tbe, the fitted temperature Tmb is significantly lower than Tbe because of two effects. 
First, the Bose distribution concentrates population at the low energies corresponding 
to the low velocity area near the center of the spectrum. This effect is tuned by 
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changing the chemical potential from the quantum degenerate regime <^ ksT) to 
the classical regime (|/i| ^ kBT, fi < 0). The fit temperature is also smaller than Tbe 
because the true energy distribution is truncated, and the Gaussian profile disregards 
truncation. This effect is tuned by changing the trap depth parameter t] = et/kBTBE 
(in practice et is fixed and Tbe is varied). For large r] truncation is unimportant, 
while for the small rj values used here truncation effects cause significant deviations. 
The calculation procedure has been checked in these limits. 

A systematic error can arise if the laser beam focus is not aligned in the trap 
as expected. The beam would be probing a region with effective chemical potential 
different from expected. Three beam parameters affect the alignment: the beam 
waist, the radial position of the beam, and the axial position of the beam. We know 
the waist radius to within 30% because of the aperture constraints encountered while 
transporting the beam into the cryostat and back out. The radial position of the 
beam is measured by translating the trap, and is verified by the presence of the 
condensate spike in the center of the spectrum. The axial position of the beam is 
not as important because the divergence length is on the order of the length of the 
sample. We conclude that uncertainties in the beam parameters should not cause a 
major systematic effect compared to the shot noise of the current measurements. 

The other indication of sample temperature comes from the phase space diagram, 
figure |5-4| . For trap B, the phase boundary appears at a ratio of trap depth to 
temperature of ~ 6.5. This would indicate that the temperature of the sample in 
figure |5-7| should be 33 fiK. A systematic error arises, however, in the extraction of 
the peak density of the normal gas (plotted in figure from the DF spectra. The 
lineshape of the spectrum changes significantly near the onset of quantum degeneracy, 
as detailed in section f).5[ This systematic error could shift the apparent density higher 
by up to roughly 30% (measured density higher than actual density), translating to 
a decrease in the indicated temperature by 20% (measured temperature lower than 
actual temperature). This adjustment shifts the expected sample temperature to 
40 /iK, within the uncertainty of the DS measurement. 

We note that the DS spectrum does not constitute an absolute measure of the 
sample temperature in the regimes of quantum degeneracy and/or shallow traps (low 
7]). On the contrary, the details of the trap shape and laser beam geometry are 
folded into the spectrum. In spite of these difficulties, reliable temperatures can be 
extracted. In addition, we observe in figure p-7| hints of quantum perturbation of 
the momentum distribution of the normal gas, a subject which could be studied with 
improved signal rate. 

Our temperature derivation procedure could be complicated by a non-uniform 
temperature in the sample. Loss of atoms from the trap occurs mostly in the con- 
densate when the condensate density is high. As explained in section |2.3.2.1| , the 
effective heating rate of the system is thus higher in the center 10% of the length of 
the trap where the condensate lives. Near the ends of the thermal gas the heat load 
is lower, and conceivably the temperature could be lower. 

In order for a temperature gradient to exist, the characteristic mean free path 
must be less than the length of the sample. At the highest densities achieved, the 
thermal gas has peak density uq = 2 x 10^^ cm~^ corresponding to a mean free path 
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Figure 5-8: Condensate fraction F = //(I + /) as a function of sample temperature, 
plotted for various bias energies. The heavy lines correspond to peak condensate 
density np = 3.3 x 10^^ cm~^, and the light lines correspond to a peak density twice 
as big. The trap depth corresponds to trap B, as do the solid lines {6 /ks = 35 yuK). A 
temperature of 50 /iK indicates a condensate fraction around 2% (12%) for Xc = Xm 
iXc = Xm/2). The calculations, based on equation |2.18| , properly treat the finite trap 
depth. 



I = l/riQcr ~ 4 cm. This is about the characteristic length of the sample. While 
significant thermal gradients are not expected, there is a danger of such effects. 



5.3.3 Condensate Fraction 

The fraction of atoms in the condensate can be obtained using two methods. First, 
given knowledge of the trap shape and sample temperature, we compute the pop- 
ulation of the normal gas using equation p.7| (which assumes thermal equilibrium). 
We then compare to the condensate population computed above using equation p.l4| . 
The trap shape parameters a and (3 cancel in the comparison, and the only quantities 
that enter are measured: the temperature, peak condensate density Up, and trap bias 
energy 9. The population ratio / = N^/Nt is given by equation |2.18| . We use / to plot 
the condensate fraction, F = //(I + /), in figure |5-8| for the trap depth, bias energy, 
and range of temperatures applicable to trap B. The condensate fraction for trap B 
corresponding to T = 50 ± 10 /iK is F = //(I + /) = 2A{t\){t\^)% (uncertainties 
arising from T and Xm listed separately). This is in agreement with estimates by 



97 



Hijmans et al. |55 . 

We have found the condensate fraction just after crossing the BEC phase bound- 
ary; condensate fractions could be determined at other points on the phase diagram, 
but the high loss from the condensate depletes the system as lower temperatures are 
attained, reducing the achievable fraction. Also, we are considering equilibrium con- 
densate fractions. As explained in section p.4| , the condensate is fed from the thermal 
cloud. An extremely high condensate fraction could be created by simply removing 
all the thermal atoms by rf ejection, but the resulting condensate would decay rather 
quickly. 

The other method of finding the condensate fraction exploits the spectrum more 
directly; the population of the normal gas is proportional to the integrated signal of 
the normal component of the DF (or DS) resonance. This integral is compared to 
the integral of the condensate portion of the spectrum, and adjustments are made 
to account for the partial illumination of the normal gas by the laser beam. (The 
condensate diameter is significantly smaller than the beam diameter, but the normal 
gas diameter is significantly larger.) The fraction computed with this technique ||T3[ 
is 2.2tJ.5%. The two methods give consistent results; an increased signal rate would 
reduce the large temperature uncertainties and make the measurements useful. 



5.4 Time Evolution of the Degenerate Gas 

Because the condensate fraction is small, there are enough atoms in the normal gas 
to continuously replenish the condensate for many seconds. As a first look at this 
feeding process we study the time evolution of the condensate. 

Simple estimates of time scales reveal that the condensate lives much longer than 
one might expect in the absence of replenishment. The large condensate density leads 
to fast condensate decay. The characteristic time an atom spends in the condensate 
before undergoing an inelastic collision is r2,c = Nc/N2^c, less than a second for the 



condensates described in table p74 However, the condensate is observed to live for 
more than 15r2,c, indicating that feeding of the condensate by the thermal cloud 
occurs quickly. In this section we examine the time evolution of the condensate. We 
find that a simple equilibrium model reproduces the data. 

To measure the condensate dynamics many sequential condensate spectra have 
been measured. Figure |^ shows part of a series of such spectra, obtained using 
DS absorption, for trap B. The solid curve is a fit to all the data in the series, with 
adjustable parameters being the peak condensate density in each scan, the overall 
amplitude scaling factors for the condensate and normal gas, the temperature of the 



sample, and a time constant for the temperature to change ||9^. The peak condensate 
density decreases only slowly in time, implying feeding of the condensate from the 
normal gas. For this feeding to occur, the normal gas must be losing atoms, and 
thus density. In order for the normal gas to remain on the phase boundary, the 
temperature must therefore be decreasing. 

The time evolution of the condensate can also be probed using the DF spectrum, 
as shown in figure |5-1CI| . Similar results are obtained. Note the systematic misfit 
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Figure 5-9: Time evolution of the Doppler-sensitive spectrum of the condensate in 
trap B. Spectra are obtained sequentially after forced rf evaporation is finished. This 
data is compiled from five identical loadings of the trap. The solid line is a model 
which is fit to the entire sequence of spectra, only four of which are shown here. p2[ 



of points just to the red of the peak shift, Ap in figures |5-10| and The spectra 



suggest that the density distribution of the condensate is more complicated than that 
predicted by the Thomas-Fermi wavefunction alone. Unfortunately, the low signal 
rate precludes further study at this time. 

Figure |5-11| shows the peak condensate density as a function of time, obtained in 
two ways. In the lower panel sequential spectra are taken of the same condensate, and 
so losses induced by the laser beam cause the sample to decay more quickly than it 
would without the extraneous one body loss rate. The upper panel reveals the latter 
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Figure 5-10: time evolution of Doppler-free spectrum of condensate in trap B. The 
noted time is measured from the end of the forced evaporative coohng cycle to the 
middle of the sweep. Each spectrum is fit independently with a two parameter model: 
the peak density Ap and an overall scaling factor Sq; the background is estimated 
independently. The right side is not fit because of signal spillover from the very 
intense DF normal line (not shown). This data is a compilation of sweeps from each 
of five identical loadings of the trap. 



situation. Here the sample is held for various intervals after forced rf evaporation is 
finished, before laser probing is commenced. The sequence of spectra obtained is fit, 
and a peak condensate density at the beginning of the laser scan is obtained. This 
density is plotted in the figure. We see that the laser probe shortens the condensate 
life by only a small factor. We also observe that in the absence of losses induced by 
the laser, the condensate density decays by only a factor of 2 over 15 s (population 
decays by factor of 2^/^ ~ 6). 

The number of atoms lost from the condensate through dipolar decay during this 
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Figure 5-11: Time evolution of peak density of the condensate in trap B. The lower 
panel shows the peak condensate density extracted from a sequence of DS spectra, 



some of which are shown in figure |5-9| . The upper panel is the peak condensate mea- 
sured after waiting various times between the end of the forced evaporation cycle and 
the beginning of the laser scans. This data demonstrates the lifetime of the conden- 
sate when the only loss mechanism is dipolar relaxation. The solid lines are from a 
model of the heating and cooling processes in the normal gas and condensate. For the 
lower panel a one-body loss rate was introduced. The inset shows the calculated sam- 
ple temperature evolution in the absence of extraneous atom loss due to the presence 
of the probe laser. 



interval (i = 15 s can be easily estimated by integrating equation p^.l5 



cAost 



'dtN2c = — -,QoOUl'^ dt nl'^it) 
' 105 2! ° io ^ ^ ' 



(5.9) 
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We estimate N^^iost = ^i-o'l)i-3) x 10^ (uncertainties from present experiment and 
20% uncertainty in Xm, respectively). If we had used Xc = Xm/2 this loss would have 
been 2''/^ ~ 11 times larger, or Nc^iost = 7{1lQ^){tl) x 10-^°. For a temperature of 
50 ± 10 fiK we estimate a normal gas population of Nt = 5^2 x 10^°. The spectra 
show that even after losing Ncjost atoms the trap still contains a sizeable condensate 
and normal gas. It therefore appears unphysical for Xc = Xnif^] however, more data 
should be collected to reduce the large uncertainties in this measurement. 

A computer model has been developed to study the time evolution of the degen- 
erate system. As described in appendix ^ this model incorporates dipolar decay and 
evaporative cooling, and assumes thermal equilibrium and a Thomas-Fermi conden- 
sate wavefunction. The good agreement between the model (solid line in figure [5-1 Ij ) 
and the data indicates that the system is reasonably well understood. An important 
result of the model is that the temperature of the sample changes by 20-30% during 
the 15 s condensate lifetime. This behavior is shown in the inset to figure [5-11 . 



The low condensate fraction and long condensate lifetime in this system suggest 
that the H condensate could be an excellent source for a bright atom laser. On the 
order of 10^ atoms/s could be continuously coupled out of the condensate for 10-20 s 
per loading of the trap (continuous loading schemes could be imagined to circumvent 
this limitation). The outcoupling could be accomplished by photon recoil as atoms 
are excited to the 25* state |^ or by Raman output coupling, as demonstrated in Na 



by Kozuma et al. The beam thus formed should be diffraction limited; the large 
condensate dimensions and short de Broglie wavelengths achievable in this system 
should result in a very low divergence beam. The atoms in the beam could be in the 
metastable 25" state, which could be useful because of the large energy available for 
driving chemistry on a surface the atoms hits. 



5.5 Signature of Quantum Degeneracy in Spec- 
trum of Normal Gas 

The direct condensate signatures presented so far are expected from simple theories 
of BEG. One does not expect the DF spectrum of the normal gas to be much affected 
by the condensate. Nevertheless, a dramatic signature appears. 

Simple arguments suggest that the for this experiment (/x <^ ksT) the normal 
gas should be only very weakly perturbed by quantum degeneracy. The condensate 
volume is about a factor of 10^ smaller than the volume of the normal gas, so the 
high condensate density should not affect the normal gas much. 

In fact, measurements of the DF spectrum of the normal gas in the degenerate 
regime reveal a surprise. The lineshape is markedly asymmetric below the transition, 
but largely symmetric above, as shown in figure |5-12| . Below the transition much 



of the spectral support is red-shifted to large detunings. These detunings, if simply 
interpreted as cold-collision frequency shift with the measured shift-to-density ratio 
Xmi correspond to sample densities a factor of two higher than expected for a degen- 
erate Bose gas in equilibrium at a temperature which is reasonable for the trap. Note 
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Figure 5-12: Doppler-free spectra of normal gas showing asymmetry in quantum 

degenerate regime. The spectra in the upper panel correspond to samples above 
quantum degeneracy; those in the lower panel are for samples that are degenerate. 
The trap shapes are identical, and the trap depths are indicated. Note the difference 
in vertical scale. The hues connecting the points are a guide to the eye. The error 
bars are statistical, and do not reflect systematic errors such as fluctuations in beam 
alignment and laser power. 



that Xm was measured in samples far from quantum degeneracy. If x were larger by 
a factor of two for the part of the cloud in the bottom of the trap where the effective 
chemical potential is small, then the asymmetric lineshape would make qualitative 
sense. The atoms in the highest density (red-shifted) part of the spectrum would tend 
to be shifted even further (because their x is larger) than their low density counter- 
parts. However, no quantitative theory has been formulated predicting that x should 
vary for the normal gas as quantum degeneracy is approached. 
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5.6 Density Shift for Excitation out of the Con- 
densate and the Normal Gas 



A thorough treatment of the relation between 
]T3[. He concludes that x 



Killian 



Xr. 



and Xm has been undertaken by 
/2 because exchange effects should be present 



in the normal gas, but absent in a condensate. Here we discuss experimental evidence 
which suggests that perhaps Xc = Xm- 



Time Evolution As explained in section 5^, the condensate is observed to persist 
for many seconds. Integration, over the lifetime of the condensate, of the rate of 
atom loss through dipolar decay in the condensate is very sensitive to Xc because 
the loss rate is proportional to n^^; a factor of two difference in Xc translates to a 
factor of 11 change in the loss rate. If we take Xc = Xruf^ then the number of atoms 
lost, Nc^iost from equation is as large as the total original population of the trap! 
The estimate here would be flawed if the trap parameters {a,f3,6) applicable to the 
condensate were somehow different from those for the thermal gas. This is unlikely, 
as noted in discussion of the trap shape in appendix ^j. 



DS-DF Comparison Doppler-sensitive excitation results in excited atoms moving 
at the recoil velocity Vrecou ~ 3 m/s. These motional states are very different from the 
low energy motional states inhabited during Doppler-free excitation. Furthermore, 
the high velocity of the recoiling atoms implies that motional coherence between 
the ground and excited states of the atom is very quickly lost during excitation, 
a situation conceivably different for DF excitation. Nevertheless, the condensate 
spectra obtained using the two different techniques are identical to within the (small) 
statistical uncertainty. 



Spectroscopic Measurement of Condensate Fraction The condensate frac- 
tion measured spectroscopically agrees well with the fraction one computes from the 
population in the condensate, assuming Xc = Xm- If Xc = Xraf^ then the calcu- 
lated condensate fraction is 2^/^ ~ 6 times larger than indicated by the spectroscopic 
technique. 



The possibility should be explored that x = '^T^{0'2S~is — o-is-is)/^ (instead of 
4^(^25-15'— c^is-is) / ^) for excitation out of the thermal gas. Then it would be easy to 
understand why Xc = Xm- Furthermore, the asymmetric spectra presented in section 
5.5 could be understood if x changed from the first to the second value in the normal 
gas near quantum degeneracy. If this were true all ambiguity would be relieved. 
Experiments designed to guarantee that the 15* and 25* atoms are distinguishable 
could conceivably provide insight into this issue. This possibility is addressed in 
section |6.3.5| . 

The Bose condensates reported in this chapter open possibilities for many new 
experiments. Some of these ideas are presented in the following chapter. 
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Chapter 6 
Conclusion 



6.1 Summary and Significance of this Work 

We have created Bose-Einstein condensates of atomic hydrogen, culminating a long 
research effort. In the process we have confirmed that, as explained by Surkov et al. 



3^ , the reduced dimensionality of the evaporation process was the crippling factor 



in previous efforts to reach this goal in our lab |^ and in Amsterdam [^]. Our 
condensates, while huge compared to other condensates, are only a small fraction 
of the total trap population, a prediction by Hijmans et al. for spin-polarized 



hydrogen in equilibrium. We have this huge number of atoms in the trap because of 
the cryogenic loading technique we use, which can be scaled to arbitrarily large trap 



sizes (consistent with the scaling arguments in section |2.3.4|) . This cryogenic loading 
method has been generalized for arbitrary paramagnetic atoms and molecules by 
Doyle and coworkers ||9^, |9^, |100|| , and therefore it should be possible to create 



big condensates of arbitrary bosons. The high condensation flux we observe (10^ 
atoms condensed per second for 10 s) should be useful in creating a bright atom laser. 
In addition, we have added another atom to the list of Bose-condensed species. In 
section [2.3| we explored many differences between the hydrogen system and the alkali 
metal gases that have been condensed so far. Finally, we have developed and applied 
spectroscopy as a new probe of condensates. 



6.2 Improvements to the Apparatus 

The experiments described in this thesis were frustrated by several factors. Changes 
in the apparatus could allow us to investigate several new areas. 

One of the most straightforward improvements is to reduce the stray electric fields 
in the cell. Large fields reduce the lifetime of the 25* atomic state by Stark quenching; 
the lifetime in a static electric field E (neglecting the natural lifetime) is ||101| 



/475 V/cm\^ , , 

TSQ = ^—j 1.6 ns. (6.1) 
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In the current cell, with its non-conducting walls, static charges reduce the 25* lifetime 
to 0.1-1 ms, but in previous experiments in a metal cell we have observed lifetimes 
on the order of the natural lifetime, Tq = 120 ms. It should be possible to reduce stray 
fields by more than an order of magnitude by coating the inside of the plastic cell with 
a thin metallic film. Longer 25* lifetimes would allow us to address questions such 
as how interactions with the condensate affect the 25* lifetime. Also, low stray fields 
would allow atoms ejected from the condensate in an atomic beam to travel many 
centimeters toward the detector before decay (see section |6.3.3|) . Reduced stray fields 
would allow optical excitation to remain coherent for many milliseconds, allowing 
interesting studies of the excitation process. The large stray fields increase indirectly 
the background counts in our data: because of the short lifetime, fewer atoms can 
be accumulated in the 25' level before the detection pulse, and so more detection 
pulses must be used to accumulate the same signal. Each detection pulse increases 
the number of background counts. Reduction of stray fields would thus increase our 
sensitivity. 

Another straightforward improvement is to reduce the heat load on the refrigerator 
caused by high power in the rf coils. The rf evaporation could then begin at higher 
sample temperatures, and more efficient cooling would undoubtedly result. More 
atoms in the condensate would increase the signal rates, and more detailed studies of 
the condensate would be possible. The heat load created by the rf coils appears to 
not be a consequence of the fields themselves, but rather of the system that conveys 
the power to the coils. If the coils were thermally anchored to the 4.2 K ^He bath the 
heating of the refrigerator should be dramatically reduced. 

Better knowledge of the trap shape and laser alignment could allow more detailed 
studies of the infiuence of quantum degeneracy on the momentum distribution of the 
sample, as developed in section |6.3.2| . As noted by Killian |jl3|, better knowledge of 



trap and laser geometry would also allow more accurate measurements of x- We could 
better understand the sizes of the condensates we produce (recall how in section ^.4| 
the 20% error on Xm translates into nearly a factor of two uncertainty in the number 
of condensate atoms lost during the lifetime of the condensate). The evidence for 
Xc = Xm presented in section could then be either bolstered or refuted, resulting 
in better understanding of the system. 

The experiments probing the condensate, described in chapter ^, suffer from low 
signal rate. Straightforward improvements to the apparatus should allow at least a 
factor of 10 increase in detection solid angle, thus allowing studies at higher time 
resolution to be feasible. Condensate formation and feeding could be studied in more 
detail. The detection solid angle could be improved by moving the cold atoms into a 
short trap formed by small currents in thin coils. This small trap could be surrounded 
by Lyman-a detectors. 



6.3 Suggestions for Experiments 

The success of the experiments described in this thesis has fertilized thinking about 
many new experiments that could be done. In this section we discuss several of these 
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ideas. 



6.3.1 Creating Bigger Bose Condensates 

Nearly any experiment is easier with a larger condensate, and so here we discuss 
strategies for putting more atoms into a hydrogen condensate. 

First, more atoms could be initially loaded into the trap. We do not currently 
understand the loading process in detail. Investigations are underway to explore 
possibilities for loading more efficiently from the discharge. In addition to increasing 
the loading efficiency, one could also simply increase the physical volume of the trap 
by making it longer and of larger diameter. One should be able to load a factor of 
four more atoms by scaling up the apparatus using available technology. 

The overall phase space compression efficiency exponent for cooling from loading 
conditions to BEG is 7 = —1.5, which is small (see section A.2.3| ). Improvements to 



the cooling efficiency could be made using two strategies. The first is to commence 
rf evaporation earlier in the cooling cycle. The second is to increase the ratio of good 
collisions to bad collisions. It is probable that the admixture of an alkali "moderator" 
atom "M" in the sample could dramatically increase the effective collision rate without 



significantly increasing the dipolar decay rate I^Tj. Presumably the cross-section for 
H-M collisions would be much greater than the anomalously low H-H cross-section. 

Finally, improvements could be made to the cooling procedure that do not require 
any changes to the apparatus. Instead of using the rf cooling paths shown in figure 
^-4| , the sample could first be compressed at constant trap depth, the depth being 
set by the highest rf frequency usable in the apparatus. Phase space compression 
would be accomplished at higher temperature, and thus more efficiently. After this 
compression the sample could be cooled to BEG at high density. 



6.3.2 Probing the Velocity Distribution of the Normal Com- 
ponent of a Degenerate Gas 

In section |5.3.2| we described the infiuence of quantum degeneracy on the velocity dis- 
tribution of the normal gas. The expected spectrum is narrower and more cusp-shaped 
than for a classical gas at the same temperature. Here we consider an experiment 
designed to compare these spectra. 

In order to ensure that the two samples (classical and quantum) being probed 
have the same temperature, we propose to use a single sample, but probe two parts 
of it. The inhomogeneous trapping potential serves to reduce the sample density 
(and thus the local chemical potential) in regions far from the deepest point in the 
trap. Glassical behavior can be observed in a gas containing a condensate when the 
probe laser is translated a distance pp ^ ksT/a from the trap axis (for an IP trap 
with T > O/ks). The beam diameter should be small compared to the translation 
distance, a condition that is easy to fulfill in an open trap (small a). 

The signal rate is dramatically reduced by the small densities far from the middle 
of the trap. In order to obtain the same signal rates as the present experiments 
(in which the high density central portion of the sample is probed), the detection 
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efficiency must be increased by the ratio of the densities, n(p = 0)/n{pp). For trap B 
in table |5.2| , and for pp = SksT/a = 160 pm, the required improvement in detection 
efficiency is about 



n{p = 0) ^3/2(1) ^^1/2(^,0) ^^^2 



n{pp) ^3/2(e-3)Ti/2(r7-3, -3) 



10^ (6.2) 



Here we have used equation |C.2| and figure |C-1| . A redesign of the apparatus for better 
optical access is required. 



6.3.3 Creation of an Atom Laser 

Atoms excited to the 25* state by co-propagating photons receive a momentum kick 
Pr = 2h/X where A = 243 nm is the excitation wavelength. The atoms leave the exci- 
tation region with the recoil velocity Vr = 3.28 m/s added vectorially to their initial 
velocity. In the normal gas there is a wide distribution of velocities, so only a tiny 
fraction are Doppler-shifted into resonance for a given detuning. In the condensate, 
however, the velocity distribution is very narrow, and one might expect to interact 
optically with all the condensate atoms at the same laser tuning. Atoms excited 
out of the (coherent) condensate by a laser beam with a long coherence time should 
constitute a coherent atomic beam. The divergence of this beam of atoms should 
be diffraction limited by the size of the excitation region. Such a scheme has been 
demonstrated by Kozuma et al. |Q using Raman pulses. 



Unfortunately, several complications arise. First, the condensate spectrum is 
broadened by the cold-collision frequency shifts arising from the high and inhomoge- 
neous condensate density. Since these shifts are much greater than the laser linewidth, 
the laser only interacts with a small fraction of the condensate for a given detuning. 
(One can imagine efficient methods of broadening the laser linewidth to match the 



condensate linewidth ||102|| , resulting in excitation of all the condensate atoms, but 
this broadening would write decoherence onto the ejected atoms.) Since the optical 
radiation only excites atoms from regions of a given density (corresponding to the 
detuning), the excitation region can be quite small, and the divergence of the beam 
is enlarged. For example, atoms excited by a laser beam of linewidth 2 kHz that 
is tuned to the peak condensate density will come from a region of diameter on the 
order of Uexc = 1 A^m- The resulting divergence angle is h/priiuJexc = 4 x 10~^ ||103 



The situation is more complicated when the laser is detuned to interact with atoms 
in regions of intermediate density. Here the excitation occurs on a shell; interesting 
diffiaction effects should occur. 

As the 25* atoms propagate with recoil velocity Vr, the potential surface they ex- 
perience is strongly perturbed by the strong attractive mean field potential created 
by the high condensate density. The depth of this potential for the peak condensate 
densities listed in table is on the order of 50 pK. The depth is small compared 
to the 650 pK of recoil kinetic energy, so the atom has enough energy to leave the 
condensate. However, the depth is much larger than the energy in the radial coor- 
dinates, p/ks = 2 pK. One expects focusing effects as the recoiling atoms leave the 
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condensate. 

The 1S-2S excitation process does not effect the hyperfine state of the atom; the 
recoihng atoms are still in the d-state, and are trapped. As the atoms move toward 
the magnetic barrier that forms the end of the trap they slow down. Dephasing also 
occurs because the integrated action for slightly different atomic trajectories quickly 
becomes quite different in the mountainous potential surface the atoms explore as 
they leave the trapping region. In order to circumvent these problems of slowing and 
dephasing we propose to flatten the potential surface: we put the atoms in a magnetic 
field insensitive state by driving the d-a transition for the 2S atoms. The transition 
can be driven selectively for only the 25* atoms because the hyperfine splitting of the 
25* state is 2'^ = 8 times smaller than the IS" hyperfine splitting of 1420 MHz. One 
might be concerned that the rf field would drive d-c state transitions in the cold cloud 
of 15" atoms, and thus change the effective trap depth. This does not occur because 
the frequency is too large (178 MHz corresponds to a trap depth of 8.5 mK). The 
a-state 25 atoms are free to fall out of the trapping region. In the current apparatus 
the atoms can fall a height h = 20 cm before hitting a window. This should take a 



time Tfaii = {—Vr/g){l — yl + 2gh/v^) = 56 ms for atoms with downward directed 

recoil velocity, and Tfaii = {vr/g){l + yl + 2gh/v^) = 790 ms for those with upward 
directed recoil velocity. Here g is the acceleration of gravity. We do not expect to 
detect the atoms launched upward because the natural lifetime (120 ms) is too short, 
and because, in the current apparatus, the atoms will hit a surface before reaching 
the peak of their 55 cm high trajectory. 

Detection of the atomic beam could be accomplished by quenching the atoms 
with an applied electric field pulse when the atoms are at the bottom of the cell 
near the detector. To estimate the signal size, we compare to the signal obtained 
when the atoms are not allowed to move very far before detection, the experimental 
procedure for figure [5-5|. We expect a reduction in the number of fluorescing atoms 



by e~'^f'^"^'^° ~ 0.63 because of decay during the transit time Tfaii- We expect a factor 
of two further reduction because only 1/2 the atoms have downward-directed recoil 
velocity. We expect an increase in signal because the detection solid angle is greater 
near the bottom of the cell. With the current cell design we should be able to achieve 
a factor of seven increase in detection efficiency, but with straightforward extensions 
of the cell closer to the detector this factor could be another factor of ten larger. 
Neglecting this improvement, we expect the atoms in the atomic beam to yield about 
the same signal size as observed in figure The experiment is thus feasible in a cell 
identical to the current one, but a conductive coating to reduce stray electric fields is 
essential. 

In addition to observing an atomic beam ejected from the sample, another in- 
teresting experiment would be a study of coUisional quenching of the 25* atoms as 
they move through the high density condensate. One might expect coUisional de- 
excitation to occur on a time scale comparable to the characteristic collision time: 
Tcoi = ^/npa2s-isVrecoii — 700 /xs. A brief study indicated no coUisional quenching 
on the sub-millisecond time scale, but the study should be repeated for longer times. 
Reduction of the stray fields would allow this. 
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6.3.4 Measuring the Condensate Density through Dipolar 
Decay 

In the experiments described in chapter |^ the condensate density was measured spec- 
troscopically, but confusion arises around a factor of two in the conversion constant Xc 
that relates frequency shift to sample density. We propose an independent measure 
of condensate density which should resolve the ambiguity. The proposed measure- 
ment exploits two-body dipolar decay, which proceeds at a well known rate in the 
condensate [Q. If the condensate density distribution is well approximated by the 
Thomas-Fermi profile in equation |2.13| , then the characteristic loss rate is given by 
equation p.l5| . Using equation |2.14| we find that the peak shift (frequency shift of 



maximally shifted part of the spectrum) should decay in time as 

35xc 

where Ap q is the initial peak shift. The only free parameters used to fit this equation 
to measured data would be Ap q and Xc- The factor of two ambiguity in Xc should be 
easy to resolve. If the Thomas-Fermi density distribution is a poor approximation, 
then the factor 4/35 should be different. The quality of the fit of the Thomas- 
Fermi lineshape to the data in figures |5-5|, p-9|, and |5-10| indicates that this 



approximation is reasonable. 

To measure the condensate decay time cleanly, we require that the condensate is 
not being replenished by the normal gas. The most straightforward way to accomplish 
this is a complete removal of the normal gas by rf ejection. The trap depth would be 
quickly swept down to the edge of the condensate, using high rf power to efficiently 
outcouple all the normal atoms. DF or DS spectroscopy of the condensate would then 
probe the time evolution of the peak shift, Ap(t). 

Since the lifetime of the condensate would be short, there would not be much 
signal accumulation time. The detection efficiency in the current apparatus is too 
low, but a factor of six improvement should make this experiment feasible (a factor of 
six can be obtained by replacing the Lyman-a filter over the MCP with a window). 
Currently, sufficient signal can be obtained by summing 0.8 s sweeps from several 
identically prepared samples. A factor of six improvement would allow the sweeps to 
be shortened to 0.13 s, short compared to the decay time ?ibXc/ ^Q^pfl — Is expected 
for Xc = Xm/2 and Ap o ~ 600 kHz. 



6.3.5 Exchange Effects in Cold-ColUsion Frequency Shift 

In section ^]B|we proposed that perhaps Xc = Xm, which would make sense if Xm did 
not include exchange effects for some reason. This could be tested by measuring the 
cold- collision frequency shift in the classical regime for atoms that are known to be 
distinguishable. To ensure distinguishability, the atoms could be put into different 
hyperfine states, c and d. A c-state atom interacting with a cZ-state atom should 
follow an interaction potential that is very similar to the identical particle case, but 
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the mean-field interaction energy would be smaller by a factor of two because the 
atoms are distinguishable and exchange effects are absent. One expects two peaks 
in the DF spectrum: one from collisions between identical particles, which is shifted 
twice as much as the other, which is from collisions between distinguishable particles. 
If the observed spectrum contains only a single peak, then we conclude that the atoms 
are distinguishable in all collision processes. In that case we expect Xc = Xm- 

We propose to do this experiment in a very deep trap {et/ks ^ (/iB//cB)506 G = 
34 mK) so that the c-state atoms will be confined (see figure We propose to 

use a sample typical of what exists immediately after loading the trap (T = 40 mK, 
no = 2 X 10^^ cm~^). If half the atoms are transferred from the d to the c state, the 
expected decay time of the c-state atoms is r ~ 2 s. The dominant terms in the loss 
rate are Gcc-*aa = 1-4 x 10~^^, Gcc-*bd = 7 x 10"^'*, and Gcc^ac = 6 x 10"^'' cm^/s for 
fields around 600 G [|104|] , which corresponds to the thermal energy. We must measure 
the DF spectrum quickly compared to this time r to observe how the spectrum 
changes. 

Time-of-fiight broadening would be significant in this spectrum, on the 40 kHz 
level for the current laser beam waist (100 /im diameter). This is large compared to the 
expected cold-collision frequency shifts of 2 and 4 kHz. To reduce the broadening, the 
beam diameter could be expanded by a factor of ten. The excitation rate (proportional 
to the square of the radiation intensity) would drop by 10^, but the interaction volume 
would increase by 10^ (the sample diameter is larger than the beam diameter). The 
linewidth would narrow by a factor of ten because of the increased interaction time, 
and so the overall signal rate would drop by only a factor of ten. Samples can 
be created every two minutes and systematics are easy to control, so it would be 
straightforward to accumulate signal. 



6.3.6 Measuring ais-is 

The ground state s-wave scattering length has been calculated to the level of 10^^ 
||40|| , but has not been measured nearly as well. One could use the interaction energy 
UoUp = Anh'^anp/m in the condensate to measure a, given the peak density Up. We 
propose to measure the chemical potential /x directly in an experiment similar to the 
measurement of Xc in section 6.3.4 . Using /i = UoUp we obtain a = fim/Anh'^np. To 
get fi we measure the difference between the potential energies at the center (6) and 
edge (£e) of the condensate. We measure 9 by sweeping the rf frequency upward until 
atoms begin to fall out of the trap, as described in section |4.2.4| . We measure Ee in 
a series of measurements of the sample density; the normal atoms are removed by 
sweeping the rf frequency down to some ending frequency fend, which is varied. For 
fend slightly larger than Ee/h the density measurements should give identical results. 
But for fend < Se/h the condensate is depleted, and the measured density decreases. 

A measurement of a to the 10% level seems feasible. Consider trap B and a peak 
condensate density np = 3 x 10^^ cm~^. The bias field 9 must not drift more than 
1 mG between measurement of 9 and Ee-, this stability should be possible. To measure 
Ee we vary fend and look for changes in the maximum frequency shift, as described 
above. Since the maximum frequency shift scales linearly with the chemical potential. 
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we must measure Ap o to better than 10%. This should be possible with the factor of 
six detection efficiency improvement outlined above. Many spectra would be obtained 
as the condensate decays, and so a group of many data points would be fit to the two 
parameters in equation |6.3| . Precision of 10% seems feasible. 



6.3.7 Dynamics of 2S atoms in IS condensate 



As described in section f).3.1\ and in the thesis of Killian [|T3] , condensate atoms ex- 



cited to the 25* state are confined radially by the deep interaction potential created 
by the large condensate density. The motional eigenstates in the 25" potential are 
rather different from the IS* eigenfunction, the condensate wavefunction. We there- 
fore expect the 2S atoms to exist in superpositions of several of these excited state 
eigenfunctions, and to thus oscillate across the condensate. 

The initial position of the 25* atoms is tuned by the detuning of the excitation ra- 
diation field. The position after some delay time r could be measured using Balmer-a 
absorption imaging. One expects breathing modes. Coherence effects in the excita- 
tion process could be explored with this technique. 



6.4 Outlook 



The studies of Bose condensed hydrogen presented here constitute simply the appe- 
tizer in the feast of experiments to be carried out. Once ambiguities about Xc are 
straightened out and the fascinating asymmetric fines in section ^]5| are understood, 
the creation of the high-flux coherent atomic seems to be the next step. The beautiful 
atom laser experiments of Mewes et al. [ 105|| , Kozuma et al. |9^, and Anderson and 
Kasevich should provide solid stepping stones. 
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Appendix A 

Theory of Trapped Classical Gas 



In this appendix we present the statistical mechanical treatment of a trapped ideal 
gas, consider decay of the sample through relaxation of the spin-polarization, and 
calculate the effects of evaporation. The appendix exists to familiarize the reader 
with concepts and scaling laws useful for hydrogen trapping and cooling experiments. 
Simple relations describe samples trapped in deep potentials. Here we address the 
effects of the finite trap depth, following the work of Luiten et al. |26|. We find that 
the truncation corrections are only important for shallow traps. 



A.l Ideal Gas 

We begin our discussion by describing a trapped ideal gas in the classical regime. 
"Trapped" implies confinement by an inhomogeneous potential with a depth signifi- 
cantly greater than the characteristic energy of the gas particles. "Ideal" means that 
collisions among the atoms, although crucial for maintaining thermal equilibrium, are 
not included in the Hamiltonian. The interaction energy for the densities studied in 
this thesis will turn out to be very small compared to the thermal energy, except in 



the case of when a condensate, as described in section p.2.2| . "Classical" means that 
the gas is far from quantum degeneracy, and the Maxwell-Boltzmann distribution 
adequately describes the correct quantum distribution. These approximations are 
appropriate for the trapped hydrogen gas studied in this thesis as it is cooled from 
the initial temperature of about 40 mK to near the onset of quantum degeneracy 
around 50 /xK. 



In sections [A.l.l| and [A.l.^ we describe the occupation function for the gas and the 



density of states functions for traps of experimental interest. In section |A.1.3| these 



concepts are employed to calculate various thermodynamic properties. In section 
A. 1.4] we consider trap losses arising from decay of the spin-polarization of the trapped 



gas. Collisions and evaporative cooling are the topics of section [A.2| . 
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A. 1.1 Occupation Function 



We assume that the gas is well described by the Maxwell-Boltzmann distribution, 
except that the energy distribution is truncated at some trap depth energy ej. The 
gas is said to be quasi-thermal. This is a reasonable assumption if two conditions 
hold First, the elastic collision rate must be high enough to redistribute particle 



energies quickly compared to the rate at which processes occur which would take the 
gas out of equilibrium, such as changes in the trap shape or depth. This is called 
the "sufficient ergodicity" condition. (Elastic collisions are discussed in section [A. 2 



below.) Second, particles with energies greater than the trap depth must be removed 
promptly compared to the rate at which they are promoted to the energetic states. 



Situations in which this condition is not satisfied have been discussed in section 2.1 



The phase space occupation function, /(e), is the number of particles per unit 
volume of phase space in a region of phase space with total energy e = 7i(p, r). {Ti 
is the Hamiltonian of the particle). We postulate an occupation function /(e) that is 
the Boltzmann factor below the trap depth, and zero above: 

/(.) = I T : : ; :: (^-d 

where A = h/y/27nnkBT is the thermal de Broglie wavelength and hq is a reference 
density; the normalization is the total number of trapped atoms. Evidence support- 
ing this postulate was obtained by numerical integration of the Boltzmann transport 
equation [^. The occupation function of a simulated gas, described initially by 



/(e) =const. below threshold e^ and zero above, was observed to transform to the 
Boltzmann form in equation [A.l| in a few characteristic collision times. The trapped 



gas is described by two parameters, T and uq, which correspond to sample temper- 
ature and maximum sample density in an infinitely deep trap [^. For finite trap 
depths, temperature is poorly defined because the distribution is truncated, and uq 
is larger than the actual maximum sample density, as discussed in section ^. 



A. 1.2 Density of States 

A useful concept is the density of states, the volume of parameter space corresponding 
to the fulfillment of some condition. We use two density of states functions. The 
"total energy density of states" function, p(e), is the differential volume of phase 
space corresponding to a total particle energy e. (The number of particles N{e)de 
within an energy de of e is the occupation per unit volume multiplied by the volume, 
N{e)de = f{e)p{e)de.) The "potential energy density of states" function, g{e), is the 
differential volume of real space with potential energy e. In this thesis e will always 
denote a potential energy, and e will always denote a total energy. Likewise, g{e) will 
always be a potential energy density of states, and p(e) will always be a total energy 
density of states. 



114 



The total energy density of states function is 



P(e) = Y,jd'v dh 6{n{p, r) - e) (A.2) 



where the Hamiltonian is ?^ = /2m + V^(r), and V^(r) describes the trap poten- 
tialQ. Using p(e), integrals over phase space of a function F(H{r, p)), which depends 
only on the total particle energy, can be converted from the six dimensional integral 
/ d^rd^pF{T-C{r,p)) to a one dimensional integral over energy / dep{e)F{e). 
The potential energy density of states function is 



g[e) 



jdh5{V{v)-E) (A.3) 



The trap shapes, ^(r), in our experiments are well approximated by either of 
two model potentials. Both potentials are cylindrically symmetric, and we use z 
for the axial coordinate and p for the radial coordinate. The aspect ratio of the 
confined clouds (length/diameter) is typically between 10 and 400. The "quadrupole 
trap with hard wall ends" [[106| , [1^ approximates traps used at higher temperatures 
(T > 1 mK), and was treated by Doyle |T^. The second model trap shape is often 



called the "loffe-Pritchard" (labeled "IP"). It corresponds to the potential we 

use at the lowest temperatures and has the form 



Vip{p, z) = ^{apf + + ef - e (A.4) 

with radial potential energy gradient a, axial potential energy curvature 2/? (units 
of energy/distance^), and bias potential energy 6. The potential energy density of 
states for this trap is 

Att 

g{e) = ^^^e{e + e). (A.5) 

We define the prefactor Aip = An ja^^f^, with units (volume) (energy) The total 
energy density of states for this trap is 

We define the prefactor Aip = n'^m^^'^Aip / h^2^^'^ , which has units of (energy)"^. 



A. 1.3 Properties of the Gas 

Properties of the gas, such as population, energy, and density, may be calculated 
using the phase space occupation function and density of states functions |^. The 

^For a particle with magnetic moment fj,, the potential energy is V{r) = fj, ■ B(r). If the particle 
is moving slowly enough so that the Larmor precession frequency changes only slowly (lu/lo <C 1), 
then the magnetic moment will adiabatically follow the direction of B(r), and V{r) = fj,B{r). 



115 



number of trapped atoms, A^, is 



TV = rdef{e)p{e) 
Jo 

= uqA^Z (A.7) 
where the single particle partition function is 

Z= r de p{e)exp{-e/kBT) (A.8) 
Jo 

The A^-particle partition function for a non- interacting classical gas is Z]\[ = Z^ /N\. 
For a power law density of states given as p(e) = po e^^^^"^ we obtain 

Z = Z^P{5 + ^/2,7]) (A.9) 

where rj = et/ksT is a dimensionless measure of the trap depth (typically 4 < rj < 10) 
and 

2oo = Po{kBTY^'''^T{5 + 5/2) (A.IO) 

is the partition function for an infinitely deep trap. The truncation correction factor 
P{q,ri) is defined in equation |A.12[ 

The partition function for the IP potential is[| 

where (f) = O/ksT is a dimensionless measure of the trap bias energy and rj = et/kBT is 
a dimensionless measure of the trap depth. We have introduced the function P{q, rj), 
which is a correction factor for the finite trap depth |2^. It is related to the incomplete 
gamma function T{q,7]) = f^'^e'^dt by 



m 

T{q) - T{q,r]) 



(A.12) 



m 

For large rj the correction is negligible, and P{q,ri) = 1. The behavior of P{q,r]) is 
shown in figure |A-1| for various values of q. The parameter q describes the energy 
exponents appearing in the trapping potential. As stated in equation |A.11| , the values 
g = 3 and g = 4 are of interest for the IP trap. Note that for rj > 7 there is only 
a small deviation from the infinitely deep trap (corresponding to r] oo) for the 
exponents 6 of interest in this thesis. 



■^For this and many other non-trivial integrals we use the computer program for symbolic math- 
ematics called Mathematica. 
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Figure A-1: The behavior of the truncation correction factor P{q,ri), plotted for 
various exponents q. The values g = 3 and g = 4 are used in finding the population 
of the IP trap. The value q = 3 (heavy solid line) is for a harmonic trap. 



The effective volume of the trap is 

Veff = N/no = ZA^. (A. 13) 

The energy of the trapped gas is 

E = rdeep{e)f{e) (A.14) 
Jo 

For the IP trap 

p P(5,r^) + |P(4,r^) 

Ejp = 4kBTNip————^^——-. (A.15) 

This result will be used in section |A.2.2| . The average energy per particle is E = E/N. 

The density at a position r (with potential energy e = ^(r)) may be obtained as 
follows. Using the phase space occupation function 

/(r.p) = n„A»exp(-^^M±£!Z!!!!) (A,16) 
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Figure A-2: The effect of truncation of tlie energy distribution on tlie density in tlie 
trap in a region of potential energy e, plotted for various trap depths rj = et/ksT. 



we integrate over trapped momentum states to obtain |g6| 

n(r) = noe-^('-)/*^^'^ (erf(t) - 2te-*7v^) (A- 17) 



where t = ^ {et — V (r)) / ksT and erf()f:) = 2 /q e * ds/y/n is the error function. Note 
that for an infinitely deep trap 

n(r) = noo(r) = noe-^^"-)/^^^. (A.18) 

The ratio of the density n{r) in a region of potential energy e = V{r) to the density 
'^oo(O) at the bottom of an infinitely deep trap is shown in figure |A-2| . We see that, 
for rj > 6, truncation does not significantly effect the density near the bottom of the 
trap where e is less than a few ksT. 



A. 1.4 Dipolar Decay 

In this section we calculate the effect of spin relaxation on the population and energy 
of the trapped gas. Spin relaxation is the dominant loss process in our trap. In this 
analysis we properly treat truncation of the energy distribution, following Luiten et 

Dipolar spin relaxation is a process in which the magnetic dipoles of two atoms 
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interact, leading to a spin flip of one or botli atoms [W\. If one of the interacting 
atoms makes a transition from the trapped (i-state to the a state (see figure 1^1) > 



hyperfine energy, corresponding to 68 mK, is released and both atoms are ejected 
from the trap. An atom that makes a transition to the c-state is weakly trapped in 
the radial direction, but is pulled by gravity out of the trap through the weak fields 
along the trap axis. An atom that makes a transition to the high field seeking b state 
is expelled from the trap. The rate constant for loss from the trapped d state has 
been calculated ||104]| and measured to the 20% level ||107| , p3| . The calculated rate 



constant for loss from the d state through all channels at low temperatures is 



'^{Gdd-^aa + Gdd-^ad + Gdd^ac + Gdd^cc) 

1.2 X 10"^^ cm^ s~\ 



+ G. 



dd-^cd 



(A.19) 



where G^x^^^p is the event rate for the process in which two atoms in the a and 
A states make transitions to the (p and ip states. The rate depends very weakly on 
magnetic field, but we ignore the dependence for the low temperature experiments 
which are the focus of this thesis. 



The dipolar decay rate from a region of density n is h = —gn^. Consequently, loss 
occurs most rapidly from regions of high density near the bottom of the trap. The 
rate at which atoms are lost from the trap through dipolar decay is 

iVs = -g r de Q{e)n^{e) (A.20) 
Jo 

where n{e) is the gas density in regions with potential energy e. An effective volume 
for dipolar decay, V2, is defined by 

N2 = -gnl V2 = -gnoN (^-^^ (A.21) 



The ratio V2/K// for the IP trap is plotted in figure |A-3| . The volume for dipolar 



decay is smaller than the effective volume because dipolar decay occurs preferentially 
in regions of highest density, which are restricted to the bottom of the trap. 



Since relaxation occurs predominantly near the bottom of the trap, the escaping 
atoms typically have less energy than the sample average energy per particle, E. 
Dipolar decay is thus a heating mechanism. The energy loss rate is the integral over 
the trap of the local loss rate multiplied by the local average of the total energy per 
particle: 

E2 = -g f^' de Q{e)n\e) (e + K{e)) . (A.22) 

The average kinetic energy per particle is K = 3kBT/2 in an untruncated Maxwell- 
Boltzmann distribution of velocities, but here we include the finite trap depth in K. 
We calculate the average kinetic energy per particle for particles with potential energy 
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0.40 




Figure A-3: The volume for dipolar decay compared to the effective volume, plotted 
for various bias fields = 9/kBT. The difference between large (harmonic) and 
small (linear radial confinement) arises because in the harmonic trap the density 
distribution is flatter in the middle. For a box potential the density distribution is 
completely fiat, and the two volumes are identical. 



£ in a trap of depth e^: 



K{e) 



2m 



kBT 



Qt + At^ - 3s 
At -2s 



(A.23) 



where s = A/vrerf (t) exp()f:^) and t = y {et — e)/kBT. Figure |A-4| shows the behavior 

of K. The average energy carried away per atom is C2kBT = E2/N2- This is plotted 
in figure A- 5 for the IP trap, for various values of 0. We see that for rj > 6 we can 
reasonably ignore truncation effects. 
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Figure A-4: The mean kinetic energy of atoms in a region of given potential energy, 
plotted for samples in traps of various depths. The trap appears shallower for atoms 
at higher potential energies. 



A. 2 Collisions and Evaporative Cooling 

So far we have not addressed elastic collisions, but they play a crucial role in the 
dynamics of the trapped gas. In this section we examine we examine the role of 
collisions in maintaining thermal equilibrium and in evaporatively cooling the sample. 
We calculate the equilibrium temperature of the gas. And we conclude the section 
with a discussion of an impediment to efficient evaporative cooling, the impediment 
that prevented the creation of BEC in hydrogen before this work. 

Only low energy collisions occur in our cold samples, so the elastic scattering 
process is described by quantum mechanical s-wave scattering theory^ |108|, The 



result of this theory is that elastic interactions are parameterized by a single quantity, 
the s-wave scattering length a. The cross section for collisions between identical 
particles is a = Svra^. In a gas of density n and mean particle speed v = JskBT/irm, 



the collision rate is Tcoi{n) = nav\/2. The collision rate determines the time scale for 
energy redistribution in the gas. 

A sample in a trap of finite depth is never in complete equilibrium since collisions 
constantly tend to populate that portion of the equihbrium energy distribution which 
lies above the trap depth. The removal of these energetic atoms is called evaporation. 



^s-wave scattering is allowed for identical bosons, but is prohibited for identical fermions 
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Figure A-5: The average total energy per atom that is lost from a loffe-Pr it chard 
trap through dipolar relaxation, in units of /c^T. The parameter cf) = d/ksT is the 
strength of the bias field, and rj = et/ksT is the trap depth. We see that we can 
ignore truncation effects for rj > 6. For the harmonic trap (large 0) the particles have 
less potential energy because the potential is quadratic; in the linear radial potential 
(small 0) there is more potential energy per particle. 



Its use as a cooling mechanism for a trapped atomic gas was first proposed by Hess 
and was first demonstrated in the MIT spin-polarized hydrogen group ||109 



Evaporation occurs when a collision of two atoms transfers enough energy to one 
atom that it may escape the trap. This atom leaves the trap with an energy e > e^, 
much greater than the average energy per atom in the sample, E. The average 
energy per remaining particle is lower, and after rethermalization of the trapped gas, 
the temperature is lower. Cooling has occurred. Note that rethermalization is given 
a headstart by the preferential evaporation of atoms from the upper end of the energy 
distribution, closer in energy to the trap depth, et, as described in appendix M. 



A. 2.1 Evaporative Cooling Rate 

The evaporative cooling rate for a sample described by a truncated Boltzmann oc- 



cupation function has been calculated by Luiten et al. using the Boltzmann 
transport equation. We quote their results here. 
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Figure A-6: The quantity Ki;ap/K// for an IP trap as a function of r/, plotted for 
various bias fields (j) = 6/kBT. The inset shows KmpC'', which becomes linear in rj for 
large rj. This linear relation is useful when estimating evaporation rates. 



The particle loss rate due to evaporation is 



evap 



-TIq (7 V Vevap 



-N 



-L col (no) K 
V2 V, 



evap 



(A.24) 



a product of the total population, a characteristic collision rate, and the fraction of 
collisions which produce an atom with enough energy to escape. The ratio Kiiap/K// 
for an IP trap is plotted in figure |A-6| for various = O/ksT (we use equation 43a 



from 26 



The energy loss rate due to evaporation is given by equation 40 in ||26|| : 



E. 



evap 



-^evap^BT {tJ -\- 1 evap / ^evap} 



(A.25) 



where X^mp is a truncation correction factor. Figure |A-7| shows the behavior of 
Xevap/Vevap for th© IP trap (we use equation 43b from [^). For large rj, the quantity 
in curly braces in equation |A.25| is nearly rj + 1. 



A. 2. 2 Equilibrium Temperature 

The temperature of the trapped gas is set by a competition between heating due 
to dipolar relaxation and cooling due to evaporation. The sample is in dynamic 
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Figure A-7: The quantity Xevap/Vevap for an IP trap as a function of r] = et/kBT, 
plotted for various bias fields = OlksT. The volume Xevap is a small contribution 
to the evaporation volumes for ?7 > 6 



equilibrium when the temperature is no longer changing, meaning that the average 
energy per particle is stationary. Setting d{E/N)/dt = dictates that 

- = ^ = ^--"^ + (A.26) 

A^ N N,,ap + N2 

If the temperature drops below equilibrium, the evaporation rate drops exponentially. 
In addition, the effective volume decreases, the density increases, and the heating rate 
due to dipolar decay increases. On the other hand, if the temperature increases above 
equilibrium, the heating rate decreases and the evaporative cooling rate increases. A 
dynamic equilibrium is thus maintained. 



The equilibrium temperature is found by numerically solving equation |A.26| . For 



notational simplicity we define the factor B which indicates whether the IP trap is 
primarily linear or harmonic in the radial direction, 

_ P(5,.,) + ,»-P(4.^)/2 
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Figure A-8: Graphical solution of the heating/cooling balance equation, plotted for 
traps of various depths ej. The bias energy is 6'//;;^ = 34 fiK. The heavy lines are the 
right hand side (rhs) of equation |A.28| . The light lines are the left hand side. The 



solid dots indicate the equilibrium temperature for the various trap depths. At lower 
temperatures the equilibrium r] is smaller, indicating less efficient evaporation. 



Then equation |A.26| may be written 

N,,ap _ 4B{r],(f)) -C2 



■ - X -■ (A.28) 

N2 r/ + l + ^-45(r/,0) 



The left and right hand sides are plotted in figure |A-8| for typical trap depths and a 



bias field O/ks = 34 /iK. The equilibrium i] is larger at higher temperatures because 
the mean particle speed, and thus collision rate, is higher. The higher 77 at higher 
temperatures indicates that evaporative cooling is more efficient there; fewer atoms 
are expended in the battle against heating caused by dipolar decay. 



A. 2. 3 Forced Evaporative Cooling 

The analysis above indicates that, for a given trap depth, a minimum temperature 
will be achieved which is set by the competition between heating and cooling. In 
order to cool the sample further the trap depth must be reduced. A comprehensive 
review by Ketterle and van Druten ||5^ of forced evaporation analyses the process 
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in detail. For the experiments described in this thesis the trap depth is reduced by 
either lowering the magnetic confinement field at one end of the trap, or by reducing 
the frequency of an RF magnetic field that selectively removes the most energetic 
atoms. This second method is the focus of chapter |^. 

A computer model (described in appendix 0) has been used to evolve in time the 
total energy and total number of particles in the sample as it is cooled. Temperature 
and density are calculated along the way. The effects of changing the trap depth 
and changing the functional form of the trapping potential have been included. This 
model agrees well with experiments. 

The rate at which the temperature of the gas can be changed is an important pa- 
rameter. Luiten et al. have shown that a classical gas with an energy distribution 



far from equilibrium will reasonably approximate a Maxwell-Boltzmann distribution 
after about three collision times. One therefore expects fluctuations of the energy 
distribution to dissipate on this time scale. To maintain quasi-equilibrium the trap 
depth and the trapping potential should thus be changed slowly compared to this 
collision time. 

Since trapped atoms are a scarce commodity, a flgure of merit for a cooling process 
is the amount of phase space compression achieved per atom lost. Far from quantum 
degeneracy the phase space density is the number of atoms per cubic thermal de 
Broglie wavelength, D = nA^. A useful measure of the cooling efficiency is the cooling 



exponent ||56[, which is the differential fractional change in phase space density per 



differential fractional change in population, 

rflogA^ ^ ^ 

A large, negative value indicates that the phase space density increases by a large 
factor for a decrease in population by only a small factor. Cooling exponents as large 
as 7 = — 3 have been observed ||110| . 



In cooling clouds of alkali metal atoms, a process called "runaway evaporation" is 
observed, in which cooling of the sample leads to an increase in sample density because 
the atoms are more localized in the bottom of the trap. In turn, the increased density 
leads to an increased collision rate and faster evaporation. The cooling accelerates, 
a desirable phenomena because the number of atoms lost through background gas 
collisions during the time used to cool the sample is reduced when the cooling time 
is shorter. The conditions required for runaway evaporation are carefully explained 
by Ketterle et al. ||5^. Runaway evaporation does not occur in clouds of trapped 
hydrogen because the dominant loss rate, due to dipolar relaxation, increases with 
increasing density with the same exponent as the evaporation rate. The ratio of the 
evaporation loss rate ("good" losses) to the dipolar decay loss rate ("bad" losses). 



Ne.ap _ CTvVe^ap . ^q^.,^^ _ 4.5) (A.30) 



is independent of density, depending only on a/T for a given rj. The effective volumes 
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for dipolar decay and evaporation used in equation A. 30 are for ~ 1 and rj > 7. 
From this equation we see that evaporation is more efficient at higher temperatures, 
as noted above. 
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Appendix B 



Fraction of Collisions Which 
Produce an Energetic Atom 

In a collision between two atoms with a medium energy, what is the probability of 
creating an "energetic atom", i.e. an atom with enough energy to escape the trap? 
Consider a collision between two distinguishable particles of identical mass m and 
initial momenta pi and p2. After the collision the momenta are labeled qi and q2. 
The total and relative momenta are P = (pi + P2) and p = Pi — P2, and similarly 
for Q and q. 

Conservation of momentum leads to the equation P = Q. We orient the axes of 
the coordinate system so that Q = Qz and we work in spherical coordinates. The 
total energy is 

Eo = {P^ +p'^)/4:m. (B.l) 

Conservation of energy demands that = q^. We thus have four equations and six 
unknowns, leaving two undetermined coordinates, which will be labeled 9 and in 
the description of 

ci — q{{x cos + y sin 0) sin ^ + i cos 9} . 

After aligning the coordinate system we see that the there are two conserved scalar 
quantites, P — Q and p = q. We can compute the fraction of this allowable phase 
space corresponding to various conditions by assuming that, in a collision, 9 and (p 
may take any value with equal probability. Wc define the function f{9, (f)) dfl as the 
volume of phase space within a differential solid angle dQ of the given angles. Since 
the collisions arc postulated to be independent of angle, wc take f{9, (f)) = fo = l/47r. 

In the calculations to follow we scale the energies to the trap depth, Ef. We use 
the dimensionless variables a = q/\/AEtm and h — Q/\/AEtm. The total energy is 
iheuEQ^Etia^ + h^). 

Wc first find the fraction of phase space corresponding to one atom with energy 
greater than the trap depth. We call this an "energetic atom". The condition is 
expressed as Ei = ql/2m > Et or E2 = q'^/'^m > Et. Using q ■ Q = qQ cos 9 and 
defining the "energy distribution" parameter ^ = ab cos 9 we translate the conditon 



129 



above to 

\^\>l-Eo/2Et. 

The fraction of phase space (consistent with the given parameters a and h) that 
corresponds to the existence of at least one energetic atom is then 

G{Q,q) = [ de sin ed(j) f {6, (j)) 

J\^\>l-Eo/2Et 

du (B.2) 

where u has been substituted for cos^ and the lower limit on the integral is a = if 
Eq > 2Et or, if Eq < 2Et, the smaller of = (2 - - 6^)/2a6 and 1. We obtain 

uo>l 

G{a,b)^{ 1-Uo 0<'Uo<l (B.3) 

1 uo<0 



This quantity may now be averaged over all the initial conditions of interest to 
obtain the average fraction of collisions which produce an energetic atom. The al- 
lowable initial conditions are that Ei < E-t, E2 < Et, and Eq > E^ since we assume 
both particles are initially trapped but there is enough total energy for one particle 
to escape. As we sum over all the possible initial states with the given a and b, we 
must weight the average by the fraction of initial states which have both particles 
trapped. This factor is simply the volume of phase space left over from the above 
computation applied to the initial conditions, 1 — G. Since P — Q and p — q, the 
average fraction is then 



Je 



d-'Qd-'ciGil-G) 

Et<Eo<2Et 

where Fq is a normalization factor to be computed below. Exploiting the spherical 
symmetry and switching to scaled variables we obtain 

fV2 ramax „ 

F = Fo (47r)2 4:Etm / b^db / a^da G{a, b){l - G{a, b)) (B.4) 

JO d'min 

where a„j,j,„ = — 6 + a/2 and qmax = ^2 — 6^ are obtained from the condition that 
< Mo < 1- The limits on the b integral are chosen to keep the limits on the a 
integral reasonable. The result is 



F = Fo 167r^ Etm 



20 - 67r 



The normalization constant Fq is simply the weighted integral over the allowable 
initial states, 

i^o= il d^'Qld^il-G)] \ 

\JEt<Eo<2Et J 
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Figure B-1: The fraction of collisions which produce an atom with enough energy 
to escape the trap. A truncated Maxwell-Boltzmann distribution of velocities is 
assumed. The parameter rj is the truncation energy divided by fc^T. 



We follow the same scheme, but note that when mq > 1 the integrand, 1 — G, remains 
at unity. We obtain 

We conclude that the fraction of the collisions which will produce an energetic atom 
is 

F = 5-^^ 0.288. 



5-1 



The above calculation assumed that the probability of having any of the initial 
states was identical. In fact, the higher energy states are less likely than the lower 
energy states, as dictated by the Boltzmann factor. Into the above integrals, then, 
must be inserted the factor exp(— Fq/A^bT) = exp(— //(a^ -|- 6^)) where r] = Et/kBT. 
The b integral must be computed numerically. We obtain behavior as indicated in 



Fig |B-1| . We see that for typical values of the parameter r], a collision which produces 
an energetic atom is an extremely rare event. 

Another quantity of interest is the probability that an energetic atom which col- 
lides with another atom will becoming trapped. This is essentially the reverse of the 
process calculated above. We conclude that in the vast majority of collisions between 
a thermal atom and an energetic atom, neither atom will be energetic afterward. 
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Appendix C 



Behavior of Non-Condensed 
Fraction of a Degenerate Gas 



In this appendix we consider a weakly interacting Bose gas of temperature T and 
chemical potential /x < confined in a trap of depth e*. A truncated Bose-Einstein 
energy distribution is assumed. Here we derive results that will be used in section p.3| . 
Truncation effects are important because of the shallow traps used for the experiments 
in chapter |. 

As discussed in section 2.2.2| , the interactions in the condensate increase the chem- 



ical potential slightly above zero. The mean-field repulsion energy of the condensate 
deforms the effective trap shape slightly, moving the energy at the bottom of the trap 
from e = to e = fi and flattening the potential there. However, we ignore these 
effects on the thermal gas because for our experiments ksT ^ fi; the magnitude of 
the deformation is much smaller than the characteristic energy of the atoms. 



C.l Ideal Degenerate Bose Gas 
C.1.1 Density 

Here we calculate the density of the gas in a region of the trap with potential energy 
e. We include truncation effects because, as will be shown in the experimental data in 
chapter ^, the typical rj is between 4 and 6 when the condensate is present. Integrating 
the occupation function over all the trapped momenta we obtain a density 

n(e;T,^,e,) = -y^ ^^^y,^^,_^y,^^ _ ^ dp (C.l) 

where Pmax/'^''^ — ^t—£ is the largest allowable kinetic energy for atoms with potential 
energy e. We separate the result of the integration into an expression for the density 
in an inflnitely deep trap and a unitless correction factor T1/2 which accounts for 



133 




Figure C-1: Behavior of the truncation correction factor T1/2, plotted for various 
effective chemical potentials. The largest deviations from unity occur for small 77', 
corresponding to potential energies within a few ksT of the energy truncation point e^. 
The range of /i' shows how the correction factor changes between the classical (/x' ^ 
— 1) and the degenerate quantum (|yu'| -C 1) regimes. The asymptotic functional form 
for low and high /i' can be seen. 



truncation of the distribution. We obtain 

,{^l-e)/kBT\ 



n{e;T,ij,,et) = 
The correction factor is 



A3(T) 



e^~^^ 



dx 



X 



gx^M' _ 1 



dx 







e /i — £ 
knT ' knT 



(C.2) 



(C.3) 



where r]' is an effective unitless trap depth and fi' is an effective chemical potential. 
In the limit of large t]', Tfc(?7',/i') = 1. The behavior of T1/2 is displayed in figure 
C-l|. As expected, the truncation correction factor is most pronounced at potential 



energies within a few /c^T of the energy threshold e^, for which r]' is small. Also, for a 
nearly degenerate gas {fi' small) the population is concentrated more at low energies 
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so the effects of truncation are less important. For trap depths greater than AksT 
the truncations effects can be neglected. 



C.1.2 Kinetic Energy 

The mean kinetic energy of the particles in a region of potential energy e is 
1 rpniax Anp"^ 

K{e- T, fi, et) = ^ipy2m+e->.)/k,T _ i ^ ^p. (C.4) 

Again we break the result into an expression for K in an infinitely deep trap and a 
correction factor. We obtain 

where rj' = (e^ — e)/kBT and fi' = {fi — e)/kBT. The first term is the classical term, 
obtained for a Maxwell-Boltzmann distribution. The middle term is the difference 
between classical and Bose statistics. The third term is the truncation correction, 
which is unity for deep traps. For fi' = the middle term is 0.513. This result is 



used in the calculation of the energy loss rate due to dipolar decay in section C.2 



This result tells us that the velocity distribution of a degenerate Bose gas should be 



narrower than its classical counterpart, an effect observed in section 5.3.2. 



C.1.3 Population 

The number of atoms in the thermal gas is obtained by integrating the Bose occupa- 
tion function, weighted by the density of states, over the trapped energy states: 

The lower bound on the integral is the bottom of the trap. For the analysis which 
follows the thermal gas will exist in the presence of a condensate, and so we set fi = 0. 
Specializing to the loffe-Pritchard trap, we obtain 

N,{T, /i, e,) = AipiksT^^^^^^^^ Ao{v, <P) (C.7) 

15 

where C, is the Riemann zeta function (C(3) = 1.202) and y4fc(r7,0) is a correction 
factor for the truncation of the energy distribution: 

ax 

ax 

- 1 
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Figure C-2: Behavior of the truncation correction factor Aq as a function of scaled 
trap depth r] = et/ksT, plotted for various bias field energies = 9/kBT. The 
asymptotic functional form for low and high can be seen. 



The behavior of Aq is shown in figure |(J-2 . 



C.1.4 Total Energy 



The total energy of the thermal gas is derived in a manner analogous to the occupa- 
tion: 



Et{T,fi,et) = / 
Jo 



e de. 



We obtain, for a degenerate gas (/x = 0) in a lofFe-Pritchard trap, 



(C.9) 



Et{T,tt)=Ajp{kBT) 



J 2^ 
\ 15 



+ 24C(5) Ai(r/,, 



(C.IO) 



where C(5) = 1.037. The behavior of A i is shown in figure G-3 . The average energy 
per particle in the thermal gas, Et, is then 



(c.ii; 
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Figure C-3: Behavior of the truncation correction factor Ai as a function of scaled 
trap depth r] = et/ksT, plotted for various bias field energies = 9/kBT. The 
asymptotic functional form for low and high can be seen. 



where 

_ ^ 07r^ + 18OC(5) 
^^'^^ = ^^ + 6OC(3)0 

. 2.701^±i^ (C.12) 
0+1.351 ^ ' 

is a factor that indicates whether the trap is predominantly linear or harmonic in 
the radial direction. This average energy is shown in figure |C-4| . For a classical gas 
the average energy per quadratic degree of freedom is kBT/2, but for a degenerate 
Bose gas the energy is less. For large the potential is harmonic and there are six 
quadratic degrees of freedom, implying a mean energy of 3k bT for a classical gas; for a 
degenerate Bose gas the mean energy is 2.7 fc^T. The difference between the classical 
and degenerate Bose results (3-2.7) is not as great as the factor ~ 0.5 mentioned for 
/i' = in the discussion of equation |U.5| . This is because in calculating the mean 
kinetic energy an average is made over all the atoms in the trap with their many 
effective chemical potentials. Atoms far from the center of the gas have fi' <^ — 1, 
and are essentially in the classical regime. 
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Figure C-4: The average energy per particle in the thermal portion of a degenerate 
Bose g function of scaled trap depth rj = et/kBT, plotted for various bias field 

energies = O/ksT. For low the potential is linear in the radial direction. More 
energy can be carried as radial potential energy for low (p than for large 0, which 
corresponds to harmonic potentials. 



C.2 Dipolar Decay 

The rate at which atoms leave the thermal gas through dipolar decay is calculated 
by integrating h = —gn"^, over the volume of the gas. For a degenerate gas in the IP 
trap we obtain a two-body loss rate 

N2AT,et) = r g{e){-g)n\e;T,fi = 0,et) de 
Jo 

= Qi{<P,V) (C.13) 

where Qi is the (unitless) integral of the density squared, and includes truncation 
effects: 

Qi{4>, r,) = v^(x + 0) ^?|/2(e~")T?/2(r7 - x, -x) dx. (C.14) 
The behavior of Qi is shown in figure |C-5| . 
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Figure C-5: Behavior of the squared density integral Qi as a function of scaled trap 
depth T] = et/kBT, plotted for various bias field energies (p = O/ksT. In the figure Qi 
is scaled by 0.337 + 0.5890, the high 77 asymptote. 



The rate at which energy leaves the thermal gas through dipolar decay is 



g{e){-g)n\e;T,fi = 0,et) {e + -f{e)) de 



A6 



(C.15) 



where 



rv 



2 5'3/2(e-^) Ti/2(77 -x, -x) 



dx 



(C.16) 



The behavior of Q2 is shown in figure P-6| . The average energy carried away by each 
particle that is lost through dipolar relaxation is then C2kBT = kBTQ2/Qi, and is 



shown in figure |C-7| . As expected, the average energy lost per particle through dipolar 
relaxation is less than E, the average energy per particle throughout the sample 
(compare figure |C-7| with figure |C-4|) , and thus dipolar relaxation leads to heating. 
Compared to the Maxwell-Boltzmann distribution, the Bose-Einstein distribution 
has more particles at lower energy. We thus expect C2 to be smaller, as is seen by 
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Figure C-6: Behavior of the dipolar decay energy loss integral Q2 as a function of 
scaled trap depth 77 = et/kBT, plotted for various bias field energies = 9/kBT. In 
the figure Q2 is scaled by 0.819 + 1.0640, the high i] asymptote. The asymptotic 
functional form can be seen for high 0. 



comparing figure |(J-7| to figure |A-5| . 
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Figure C-7: Average energy of atoms lost from thermal portion of a degenerate gas 
through dipolar decay, in units of ksT, plotted for = and various bias energies 
4> = 0/kBT. The asymptotic functional forms for small and large (j) are apparent. 
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Appendix D 

Model of the Dynamics of the 
Non-Degenerate Gas 



This appendix describes a computer model of the dynamics of the trapped gas, in- 
cluding evaporation, dipolar decay, and changes of the trap shape and depth. The C 
source code is available. 



D.l Overview 



We desire a general purpose computer model of the time evolution of our trapped gas 
to assist in choosing paths through parameter space toward BEC. Previous models 
jlTl employed approximations that are not valid at the low temperatures we wish to 
access and are not easily adaptable to the more complicated trap shapes we wish to 
employ. 



D.1.1 Approximations 

The primary approximation is that the phase space distribution function for the 



gas is always a Maxwell-Boltzmann truncated at the trap depth. Luiten et al. 
have confirmed this approximation for slow evaporation. This approximation is valid 
as long as two additional assumptions are valid: (a) instantaneous ergodicity and 
(b) instantaneous removal of energetic atoms. The validity of these approximations 
depends on some time scales. For approximation (a) we need the time required for the 
gas to re-equilibrate to be short compared to the speed at which the trap shape and 
depth are changed. For approximation (b) we need the energetic atoms to be removed 
quickly compared to the energy redistribution time. These are the approximations 



described in section A.. 1.1 



D.1.2 Describing the Trap 

We make no approximations about the trap shape. Instead the trap is described by 
energy and potential energy density-of-states functions (p(e) and Q{e)) which may 



143 



involve any number of parameters, all (possibly) time dependent. The energy trun- 
cation point is also given as a time-dependent parameter. 



D.2 Statistical Mechanical Description of the Gas 
D.2.1 Variables 

We describe the trap by an evaporation energy threshold, et(t), and the density of 
states functions and their set of (time dependent) parameters, ao-(^)- These are our 
independent variables. Some function is given by the user to relate these trap shape 
parameters to time, the ultimate independent variable. 

The primary dependent variables are E and A^, the total energy and total pop- 
ulation of the sample. At each time step during the evolution of the trapped gas 
we calculate changes to the total energy, dE, and to the total number of atoms in 
the trap, dN. The new value for T may then be calculated from E. If the form of 
p(e) is known, it is often possible to analytically relate T to E. Otherwise T must 
be determined numerically, as is done in this model. The density parameter no is 
trivially related to by equation |A.13. 



D.2. 2 Computing Properties of the Gas 

The population, A^, total energy, E, and density n{r) are computed as described in 
section |A.1| . 

For notational simplicity, we introduce the functions 



V^„= r dep{e)e''e-'/^^^ (D.l) 
Jo 



and 

ipn= I depie)e". (D.2) 



The Helmholtz free energy is given hj A = — ksT log Z^-, where is the A^- 
particle partition function defined below equation |A.(j| . We can define a "single par- 
ticle" free energy 

Ai = —ksT log Z 



so that A = NAi + ksTlogNl. The entropy is 

(D.3) 



OA 

^~ df 



V 

To keep Veff constant, we replace Z with Veff/A^ (using equation |A. 13|) to obtain 

3N 



S = kB(^NlogZ + — -logW}j. (D.4) 
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D.3 Treatment of Physical Effects 



The model accounts for cooling due to evaporation, heating due to two-body spin 
relaxation, adiabatic trap shape changes, and the skimming of energetic atoms as the 
energy threshold is lowered. The finite trap depth is fully integrated into the model. 
At each time step we calculate 

E = Eevap + E2 + Etrap + Egkim (D-5) 

and 

N = Ne,ap + N2 + NsHm. (D.6) 

Each of these terms is explained below. 



D.3.1 Evaporation 

Luiten et al. p6| have described the evaporation process, starting from the kinetic 



equation. We summarize their results here. Following section |A.2.1| , the particle loss 
rate due to evaporation is 

Nevap = a V Vevap (D-7) 



where a is the collision cross section, v = \ SkBT /-Km is the mean speed 



evap 



^1 



is the effective volume for evaporation, and 77 = tt/kBT. The energy loss rate due to 
evaporation is 



E, 



evap 



n^av {{et + kBT)Ve^ap - ksTX^^ap} 



(D.9) 



where 



X. 



evap 



A^e-" / * rfep(e) [e'^/'^s^ - (r^ + 1 - e/kBT)e 



3 —rj 



A^e 



[i)o - (r/ + 1)^06-'^ + 



^1 ^ 
kBT 



(D.IO) 



D.3. 2 Dipolar Spin Relaxation 

The particle and energy loss rates N2 and E2 due to dipolar relaxation are given by 
equations |A.2CI| and |A.22| in section |A.1.4 . 



D.3. 3 Trap Shape Changes 

We assume that the trap shape is changing slowly compared to the transit time of an 
atom across the trap. Trap shape changes are then adiabatic. We can use dS = Q to 
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compute dT from changes in the parameters to p(e), and eventually obtain dE. 
We break dE into its constituents 



,r. dE 
dE= — 
dT 



dT+Y: 



dE 



dan 



(D.ll) 



where the first partial derivative is done at constant trap shape and the sum in the 
second term is over all the trap shape parameters. Since the process is adiabatic, 
dS = 0. From equation p.4| we see that if dS = and is fixed we must have 
dZ = 0. Then we may find dT by writing out dZ. We obtain 



dT 





darj 

T 


dZ 
dT 





(D.12) 



Note that 



Remembering that E 



n+1 



dT 

Nipi/ipo we get 



E, 



trap 



E 



E 



Otr. 



(D.13) 



(D.14) 



In practice the partial derivatives are evaluated numerically. 



D.3.4 Skimming Energetic Atoms 



As the trap evaporation threshold is lowered, atoms are skimmed from the top of the 
energy distribution. Using equation |A.7| we calculate 



dNgf^j^jji 



det 



dN 
d^ 



T,a„ 



and 



E 



skim 



(D.15) 
(D.16) 



D.4 Software Implementation 

The evaporation simulator is a set of routines which numerically integrate the differ- 
ential equations for E and A^ through time. A user of these routines must write C 
routines that describe the trap to the integration routines. The description takes the 
form of an energy density of states, a potential energy density of states, and an energy 
threshold above which the energy distribution is truncated. Other code (described 
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below) is also expected to parse the command line and do initialization. Ideally, all 
of the information that specializes the general evaporation routines to a specific trap 
and ramp function is contained in a single file in the directory specialized. This 
file is mentioned in the top level Makefile. Henceforth this file will be referred to as 
userfile. 

Information about the trap is passed to the evaporation simulation routines in a 
data structure, mytrap, defined in prog.h. The components of this data structure 
are pointers to functions and arrays. The data structure is a global variable, and is 
initialized by the routine initial ize_mytrap in userfile. This routine is called at 
the very beginning of the program. 



D.4.1 Program Flow 

Execution of the program flows as follows. First, initialize_mytrap() is called to 
set up the data structure. Then a routine, pointed to by mytrap, is invoked to parse 
the command line and stash user-provided information into the appropriate arrays. 
The differential equation integrator is then started. It steps E and N through time 



using the Runge-Kutta routines from Numerical Recipes ||111|| . At each time step the 
routines from userfile are used as follows. First, the threshold and parameters used 
for calculating the density of states are calculated. These might correspond to magnet 
currents or rf field frequencies, as computed for controlling the experiment. These 
values are printed to a log file. Since these parameters may often be manipulated into 
a form that economizes computation when calculating the density of states, another 
routine, pointed to by mytrap, is invoked to do this. Finally, we are ready to calculate 
ipn and (fn for this instant. The new temperature and density are calculated. The 
effects of evaporation, relaxation, skimming, and adiabatic trap shape changes are 
determined, and we are ready for the next step. 



D.4.2 Data Structure 

The data structure mytrap is referenced often is these calculations, and is the link 
between user-supplied routines and the general simulation machinery. The entries in 
the data structure are: 

parse Pointer to a function which parses the command line, obtaining the 
initial density and temperature of the simulated gas, and the duration of the 
simulation. This information is returned to the simulator. Further information 
might be the beginning and ending currents in certain magnets, shape factors 
for ramps, etc., which would be used to specify the trap shape at each instant 
of time. 

threshold Pointer to a function that gives the energy threshold (atoms with 
more than this energy are assumed to escape from the trap), as a function of 
time. 
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calc_dos_parains Pointer to a function which calculates, at a given time, the 
parameters used to calculate the density of states. The values calculated here 
will be printed in the log file for each time step. 

n_dos_params An integer, giving the number of density of states parameters 
calculated by calc_dos_parains. These are printed in the log file. 

precalculate Pointer to a function which manipulates the density of states 
parameters to speed up actual calculation of the density of states. For example, 
in a lofTc trap the density of states is proportional to the square root of the radial 
confinement field. This square root may be calculated here in this routine. 

edos A pointer to a routine which returns the energy density of states at a 
given energy. 

pedos A pointer to a routine which returns the potential energy density of states 
at a given energy. 

trapspecs An array of information used (possibly) by calc_dos_params and 
threshold. These data are not touched by the simulator. The only C file that 
looks at these is userfile. 

ntrapspecs Size of the array trapspecs. Only used in file userfile. 
Several examples of userfile exist, and may be copied and hacked. 



148 



Appendix E 

Bose-Einstein Velocity Distribution 
Function 



To find the one-dimensional velocity distribution function which is used in calculating 
the Doppler-sensitive lineshape observed in section |5.3.2| one should use the Bose- 
Einstein (BE) occupation function instead of the more convenient Maxwell-Boltzmann 
(MB) occupation function. For a gas at temperature T and chemical potential /i the 
occupation function is 

In this appendix we derive the BE one-dimensional velocity distribution function and 
use it to calculate characteristic spectral lineshapes. Truncation effects for a trap 
of depth et are included, but shown to be small when the laser beam intersects the 
middle of the trap and preferentially excites atoms from the deepest part of the trap. 
In this case, quantum degeneracy is shown to narrow the spectrum. 



E.l Velocity Distribution 

The one-dimensional momentum distribution function is obtained in a way analogous 
to the density, in equation |C.2| , but the momentum integral is over only two dimen- 



sions. We define fpipz) dpz as the probability of an atom having a momentum along 
the z axis within a range dpz of Pz'-, the normalization is 1. The distribution is 

P±,max 27rp± 

. .u e{pi/2m+py2m+e-^.)/ksT _ ^ 

fp^P'^ = 7v J:^^ (^-2) 

g(p2/2m+£-M)/fcsr _ I 







where p± is the momentum perpendicular to z and P±max/'^''^ = et — e — p1/2m 
is the maximum allowable kinetic energy perpendicular to z. The denominator is 
h^n{e;T, ii,et), calculated above. The numerator can be broken into an expression 
which neglects truncation and a correction factor Tq defined in equation |U.3|. We 
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Figure E-1: Behavior of the truncation correction factor Tq as a function of effective 
trap depth rj', plotted for various effective chemical potentials ii'. 



obtain 



fp{Pz;T,fx,et,€) 



—27imkBT 
h^n{e;T,ix,et 



log 1 



xTfi 



et - e-pl/2m /i 



knT 




(E.3) 



For large rj' truncation is unimportant 



The behavior of Tq is shown in figure |E-1 . 
and To = 1. 

The velocity distribution function is fvivz) = mfp{mvz). Figure pi)-2| shows this 
distribution function for various potential energies in the trap. At higher potentials 
the effective chemical potential becomes more negative, and the distributions look 
more classical. Only very near the trap minimum is there significant clumping of the 
velocities around zero. 

Figure [E-3| directly compares the quantum velocity distribution (small /i) with 
the classical distribution (large /i). As expected, the classical velocity distribution 
at £ = O.lfcsT" and e = IksT are nearly identical, but differ substantially from 
their quantum counterparts. For e = 5k bT the classical and quantum distributions 
are nearly identical because even the degenerate gas becomes dilute at this high 
potential energy. Truncation effects are important, as evidenced by the sharp cutoff 
of the velocity distribution. 
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Figure E-2: One dimensional velocity distribution function, plotted for regions with 
various potential energies e (listed in units of ksT). The sample temperature is 
T = 50 jjK, the chemical potential is = — 2 x 10~^A;bT, and the trap depth is 6kBT. 
The inset reveals part of the £ = distribution that is not visible in the main figure. 



E.2 Spectral Signature of Velocity Distribution 



The Doppler-sensitive spectrum is a direct measure of the one-dimensional velocity 
distribution function fv{vz) through the (radian) frequency shift A = k • v, where 
k = 2-112/ X is the wavevector of the radiation and A = 243.2 nm. The spectrum, 
S'(A), is given by considering the spectrum of the atoms at each potential energy in 
the trap: 

S{A)=So deg{e)n{6) dvJ,{v,)S{A - kv,) (E.4) 

■'0 J —Vmax 

where = et — e \s the maximum kinetic energy allowable for atoms with 

potential energy e. The inner integral vanishes if |A/A;| > Vmax- The upper integration 
limit on the outer integral thus becomes Smax = — mA^/2/c^. The spectrum is 

5(A) = ^ 1'™^ deg{e) log (l - e('^— ^V2fe^)ABT>) 

^ [ et-e-mAy2k^ - e - mAy2k^ \ 

k^ ' k^ ^^-^^ 
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Figure E-3: Comparison of one dimensional velocity distributions for samples near 
to (far from) quantum degeneracy, having chemical potential fi/ksT = —2 x 10^^ 
{fi/kBT = —20), plotted with heavy (light) lines. The sample temperature is T = 
50 /iK. The trap depth is G/c^T. The velocity distribution for each chemical potential 
is plotted for three different potential energies, listed in units of fcsT. For the e = 5 
case the two lines overlap; also, truncation effects are clear because of the shallow 
effective depth. 



This calculation has assumed uniform laser illumination of the sample, has ignored 
coherence effects for atoms that pass repeatedly through the laser beam during ex- 
citation, has assumed all atoms are driven for identical durations, and has neglected 
the cold-collision frequency shift. 

The non-uniform laser illumination can be inserted into equation |1l).5| as a factor 



R{e), an average of the square of the laser intensity over the surface of constant 
potential energy e in a trap potential V{r), 

R{s)ll = -^Jdh /2(p, z) 5is - Vir)) (E.6) 

where I{p,z) = Ioexp{—2p^/uj'^{z)) describes the spatial dependence of the intensity 
of a Gaussian beam [ p.03|| ; the beam radius a distance z from the focus is uj{z) = 



cjoyl -|- {z/dy where ujq is the waist radius and d = nul/X is the divergence length. 
We assume the beam focus is at the center of the trap and the beam is oriented 
parallel to the cylindrical symmetry axis of the trap. For a loffe-Pritchard trap 
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potential energy, e [jiK] 



Figure E-4: The average of the square of the laser beam intensity on an equipotential 
surface of energy e, plotted for various beam waist radii Uq. The trap has radial 
energy gradient a/ks = 16 mK/cm, axial curvature jS/ks = 25 /iK/cm^, and bias 
d/ks = 34 fiK. At large e the average vanishes very slowly because the beam diverges 
and illuminates more of the surface. 



potential V{p, z) = y(«p)M-7^^^^~T^)^ — 9, we obtain 




4 (g + 9Y - {(3z^ + 
Ml 1 + {z/df 



This average has been calculated for the trap shape of immediate interest for the 
experiments described in this thesis. See figure [E-4| . The laser beam primarily excites 
the atoms at the bottom of the trap where the beam intensity is largest. The energy 
distribution truncation effects are thus not very important, and the shape of the veloc- 
ity distribution is well approximated by — \og{l— exp{{p—e—mv'^ / 2)/ ksT)). Because 
the the effective chemical potential is small for these atoms, however, quantum effects 
can be important. 

Figure [li)-5| shows calculations of the lineshape, using equation [E.5| and including 



R{e), for samples in the quantum degenerate regime at various temperatures for 
a set of trap parameters of immediate interest to this experiment. We see that 
the distributions are significantly narrower than those arising from samples far from 
quantum degeneracy. If the spectrum for a degenerate sample is fit to the distribution 
expected for a non-degenerate sample, the quoted temperature will be about 40% too 
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detuning [MHz] 

Figure E-5: Expected Doppler-sensitive lineshape for a gas near to (far from) the 
quantum degenerate regime, plotted for samples at various temperatures. The thick 
(thin) hnes are for samples with chemical potential ///fee = — 1 x 10~^(— 5 x 10^) //K, 
corresponding to the degenerate (non-degenerate) regime. For these calculations = 
300 /xK, a = 16 mK/cm, /? = 25 /xK/cm^, ^ = 34 /xK, and uq = 48 /xm. 

low; the actual temperature is about 1.6 times higher for the conditions of immediate 
interest in our experiment. 
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Appendix F 

Model of the Dynamics of the 
Degenerate Gas 



F.l Overview 



This appendix describes a computer model of the dynamics of the degenerate trapped 
gas, including evaporation and dipolar decay losses from the condensate and thermal 
gas. As depicted in figure |2-2| , the trapped gas is modeled as a system composed of 
two tightly coupled components, the thermal gas and the condensate. We assume 
the loss rates from the system are slow compared to the rates of transferring particles 
and energy between the two components. A dynamical equilibrium is thus postulated. 
The results of this model are shown in figure |5-11| . 

The trap shape is assumed to be static and of the loffe-Pritchard form. The depth 
is fixed. The Bose- Einstein occupation function is used for calculation of dipolar decay 
losses from the thermal gas, and the Thomas-Fermi condensate wavefunction is the 
basis for calculations of loss from the condensate. 

The particle and energy loss rates due to evaporation should be calculated using 
a quantum kinetic theory derived from the quantum Boltzmann transport equation. 
This approach would correctly include the Bose statistics. As a very simple approxi- 
mation to the desired, but complicated, formulation, we simply apply the evaporation 
rates calculated for a Maxwell-Boltzmann distribution by Luiten et al. |^ (equations 
in section |A.2.1|) . One expects a lower evaporation rate near quantum degeneracy be- 
cause atom energies are distributed to lower energies relative to the thermal energy 
fc^T than is typical in a classical gas. Roughly speaking, the effective t] is larger. 
Another way to think about the situation is to realize that collisions which lead to 



evaporation occur preferentially among higher energy atoms. As shown in figure 2-1 



the Maxwell-Boltzmann distribution calls for many more atoms in these regions of 
higher energy than does the Bose-Einstein distribution in the degenerate regime. We 
therefore reduce the evaporation rate by the adjustable factor T of order 0.4. This T 
is the factor mentioned in figure |5-11| . Note that the rate of evaporation is reduced, 
but the energy of the escaping atoms is still approximately {rj + VjksT . 
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F.2 Simulation Details 



The input to the simulation is the initial temperature of the thermal gas and the 
initial condensate population (the proper thermodynamic parameters are T and /i; the 
chemical potential is directly linked to the peak condensate density, which is directly 
connected to the condensate population). The trap shape and depth {a,l3,9,et) are 
additional parameters. Finally, a one body loss rate with time constant ti is specified. 
Using this input, the program calculates the total occupation of the trap, Nt, and 
the total energy, Et, using equations |C.6| , |C.9| , and |2.14| , and the condensate energy 
Ec = fJ^Nc = npUoNc- It is the total energy and the total population that is evolved 
in time. At each time step the temperature is found by numerically solving equation 
C9 for the total energy of the thermal gas. We assume /i = for the thermal gas 
during the entire simulation. If we allowed fi to vary, then two equations would 
have to be consistently solved numerically, one for Nt and one for E^ (possible, but 
more difficult). Once the temperature is found, the occupation of the thermal gas 
is calculated. The condensate occupation is the total trap occupation minus the 
thermal occupation, Nc = Nt — Nt (we postulate equilibrium between the thermal 
gas and the condensate, a questionable assumption as discussed in section |2.3.3|) . 
The particle and energy loss rates due to dipolar relaxation from the thermal gas and 
condensate are calculated using equations |C.13| and |U.15| . Possible one-body losses 
are added. The particle and energy loss rates due to evaporation from the thermal 
gas are calculated for a classical gas of the same peak density as the Bose gas being 
simulated, rio = 2.612/A^(T). These rates are multiplied by the correction factor JF. 
The sum of all the particle and energy loss rates for the given time step is returned 
to the differential equation stepper routine. 

The simulation is stepped through time using the Romberg integration technique 



111 . Errors in Nt and Et are maintained below the 10 level 



F.3 Improvements 

One important improvement would be to determine the evaporation rate using the 
quantum Boltzmann equation, which is based on the correct quantum occupation 



function. This has been done very recently [112]. 

Further improvements would include treatment of the finite feeding rate of the 
condensate from the thermal gas. This bottleneck should be playing a role for the 



largest condensates in this thesis, as suggested in section p.3.3| . 
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Appendix G 

Trap Shape Uncertainties 



In order to calculate various properties of the thermal gas and condensate, it is nec- 
essary to combine measurements (of density and temperature, for example) with 
knowledge of the trap shape. In principle the trap shape can be calculated simply 
from the Biot-Savart law. In practice various complications arise. In this appendix 
we discuss these complications and estimate their effect on the trap shape. We con- 
clude that the trap parameters (radial gradient a and axial curvature /?) should have 
uncertainties of less than 20%. 

The trap shapes used for the experiments in this thesis are created by currents 
in 17 superconducting coils. Various compensation schemes are employed to produce 
steep gradients and large curvatures along the z axis, and shim coils are used to adjust 
the bias field. One source of uncertainty is knowledge of the precise geometry of the 
coils. Additional uncertainties exist because fields of unknown strength and direction 
arise from trapped fluxes Jl^ and magnetized materials in the cryostat. 

We first address the stray fields: how can these stray fields influence the trap 
shape? Let us consider the influence of stray fields on the three loffe-Pritchard trap 
parameters, a, /3, and 9. We label the three space coordinates x, y, z, and assume the 
trap exhibits cylindrical symmetry about the z axis. A uniform field perpendicular to 
the z axis will shift the location of the trap minimum in the plane, but will not affect 
a, f3, or 6. A uniform field parallel to the z axis will shift 6, but not affect a or f3. A 
gradient field with strength B = sx {s is the gradient strength) aligned perpendicular 
to the z axis will add linearly with a, but will not affect j3 oi 9. A field with strength 
B = s z aligned parallel to the z axis will move the trap minimum along the axis and 
shift 9, but will not affect j3. A field with strength B = s z aligned perpendicular to 
the z axis will tip the trap symmetry axis away from the z axis, but not affect any 
parameters directly. Finally, a field parallel to the z axis and with strength B = s x 
will bow the trap and directly perturb all three parameters. Defining the unitless 
parameter e = Jl + {hbs/oiY, we find that a ea, (3 (3/e, and 9 ^ 9/e. 



In the traps of immediate interest in this thesis, the samples are short (a few cm) 
and very narrow (less than 1 mm in diameter). The radial gradient fields are large 
(~ 10^ G/cm), and sources of stray fields are far (several cm) from the sample. Thus, 
we expect only very weak perturbations to the gradient parameter a and the curvature 
(3. The parameter 9 is easily perturbed, but we can measure it directly, as described in 
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section [4.2.4| . Trapped fluxes are expected to exist on the 1 G level, but the distance 
from the sources (superconducting wires) to the trap should cause these flelds to 
appear as simple gradients, which have been shown above to only slightly perturb (3 
and a. Magnetic materials in the cryostat could produce large flelds (~ 10 — 100 G), 
but they are very far away and the fleld should therefore appear uniform across the 
trap (i.e. not appear as gradients). We conclude that stray flelds should not perturb 
the trap parameters a and (3 signiflcantly. The parameter 6 may be perturbed, but 
it can be measured directly. 

Measurements of 6 can be complicated by the effects of gravity: for very cold 
samples gravity can pull the atom cloud away from the region of lowest magnetic fleld, 
thus making the measurement of 6 ambiguous. This distortion becomes important 
when the gravitational potential energy varies over the height, h, of the cloud by 
many times the thermal energy: h ^ kBT/mg = (0.84 mm//iK)T. For the present 
samples we are far from this limit. 

Since stray flelds do not strongly effect a or (3, the primary uncertainty in knowl- 
edge of the trap shape arises from uncertainty of the coil winding geometry. The fleld 
cancelation scheme used to create tight curvature /3 along the z axis amplifles the 
uncertainties since the net fleld is the small difference between multiple large flelds. 
Reasonable uncertainties in the winding geometry (elongations and translations of 
coils on the scale of 5 mm) give rise to no more than 10% changes to the calculated 
/3. 

A calculation of one of the trap shapes used in chapter |^ is shown in flgure p-1 



along with analytic flts to the bottom of the potential. The harmonic approximation 
of the axial curvature of the trap should be precise to the 10% level. 

We conclude that uncertainties in the trap shape due to both stray flelds and 
imprecise winding geometry should be below the 20% level. 
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Figure G-1: Calculation of potential energy contours in trap A, along the trap axis 
(p = 0) and offset from the axis by distance p = 0.854 mm. A harmonic potential 
with curvature /3 = 25 ± 3 /iK/cm^ is shown. This curvature is used in chapter ^ 
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Appendix H 



Symbols 



Here is a compilation of the symbols used in this thesis, and the pages on which they 
are introduced. 



symbol 


description 




a 


radial energy gradient in loffe-Pritchard trap, p|115 






axial energy curvature in loffe-Pritchard trap, pi 15 






elastic coUision rate, pl21| 




7 


phase space compression efficiency parameter, pl26 




e 


total energy of a trapped atom, pll4| 






trap depth (energy), p|ll4| 




e 


potential energy, p|ll4| 




V 


scaled trap depth, et/ksT, pll6 







bias field energy in loffe-Pritchard trap, pll5| 




A 


thermal de Broglie wavelength, h/ \^27rmkBT, p|114| 






chemical potential of sample, pp4| 






Bohr magneton, e^/2me = 9.274 x lO^^^ j/t , p41 




P 


radial coordinate relative to the trap center, p|115| 






total energy density of states, p|114| 




g{e) 


potential energy density of states (differential volume with potential energy e), pll4 


a 


collision cross section for identical particles, Svra^, pl21 





scaled bias field, 9/kBT, p|ll6| 




Up 


radial trap oscillation frequency, pp2| 






axial trap oscillation frequency, p|2^ 




Ajp 


total energy density of states prefactor for IP trap, p|ll4| 


Aip 


potential energy density of states prefactor for IP trap, p|114 


a 


scattering length for collisions between two IS" H atoms, 0.648 A, p|121 


5(0) 


indicates whether trap is predominantly linear or harmonic in radial direction, pl37 




dependence of heating ratio He/ Ht on trap shape, ] 


m 




average energy of atoms lost from thermal gas through dipolar decay, p|120|, p|139 


D 


phase space density when far from quantum degeneracy, noA^, p|126| 


E 


total energy of sample, pll7 




Ec 


total energy of condensate, pPS] 
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Et 


total energy of thermal gas, p|136| 




E2,c 


energy loss rate due to dipolar decay from the condensate, 






energy loss rate due to dipolar decay from the thermal gas, p|139| 




F 


condensate fraction, Nj{Nt + Nc) = //(I + /) , p29 




f 


condensate occupation ratio, Nc/Nt , p|29| 






dipolar decay event rate for a given channel, pll9| 




9 


dipolar decay rate for a cold gas of (i-state atoms, 1.2 x 10~^^ cm^/s 


pll9 




Bose function, E^o-^V^". 






rate at which process a heats the gas, p^ 




ks 


Boltzmann's constant, 1.3807 x IQ-^s J/^^K 




m 


mass of H atom, 1.0078 amu 




N 


population, number of atoms, pi 16 




Nc 


population of condensate, p|28| 




Nt 


population of thermal cloud, p|l35| 






particle loss rate due to (two-body) dipolar decay from thermal gas. 


pl38 


n 


sample density (atoms per unit volume), p|118|, p|134| 




no 


reference density of thermal gas, p|114| 
critical density for EEC, 2.612/A3(r), p26 




ric 




ncond 


condensate density, p^ 




rip 


peak condensate density, pp7| 




Qi 


unitless dipolar particle decay rate integral, pl38 




Q2 


unitless dipolar energy decay rate integral, pl39 




T 


sample temperature, p|114|, p|24| 




Uo 


mean field energy per unit density, Anfi^a/m, p27 




V{t) 


trap potential energy function, p|ll5| 




^cond 


volume of the condensate, assuming a Thomas-Fermi wavefunction, p35 


Veff 


effective volume of the sample, pi 17 




V 


mean speed of particles in a classical gas, ^JSkBT/irm, pl21 




z 


axial coordinate relative to the trap center, p|115| 




Z 


classical single particle partition function, pi 16 
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Appendix I 

Laboratory Units for Energy 



The chart on the following page may prove helpful when comparing the various ways 
energy is manifest in atom trapping experiments. The relevant constants are: 

constant value 

ks 1.381 X 10^^^ J//iK, Boltzmann constant 

fjLB 9.274 X 10-24 J/T, Bohr magneton 

h 6.626 X lO^'^^ J/Hz, Planck constant 

mg 1.6411 x 10"^^ J/cm, gravitational energy for H 
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= 1.4 MHz/G 
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